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ABSTRACT 


We  are  interested  in  the  time-asymptotic  behavior  of  solutions  to  viscous  conservation 
laws.  Through  the  pointwise  estimates  for  the  Green’s  function  of  the  linearized  system 
and  the  analysis  of  coupling  of  nonlinear  diffusion  waves,  we  obtain  explicit  expressions 
of  the  time- asymptotic  behavior  of  the  solutions.  This  yields  optimal  estimates  in  the 
integral  norms.  For  most  physical  models,  the  viscosity  matrix  is  not  positive  definite 
and  the  system  is  hyperbolic-parabolic,  and  not  uniformly  parabolic.  This  implies  that 
the  Green’s  function  may  contain  Dirac  ^-functions.  When  the  corresponding  inviscid 
system  is  non-strictly  hyperbolic,  the  time-asymptotic  state  contains  generalized  Burgers 
solutions.  These  are  illustrated  by  applying  our  general  theory  to  the  compressible  Navier- 
Stokes  equations  and  the  equations  of  magnetohydrodynamics. 


1991  Mathematics  Subject  Classification:  Primary  35K55,  76N10;  Secondary  35B40,  35A08,  35L65, 
76W05. 

Key  words  and  phrases:  quasilinear  hyperbolic-parabolic  systems,  large  time  behavior,  Green’s  function, 
nonlinear  and  linear  diffusion  waves,  conservation  laws,  compressible  Navier-Stokes  equations,  magneto¬ 
hydrodynamics. 
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1.  Introduction 


Consider  a  system  of  viscous  conservation  laws 
(L1)  ut  +  f(u)x  =  (B(u)ux)x, 

where  u  =  u(x,t)  is  the  density  of  physical  quantities,  f(u)  is  the  flux  and  B{u )  is  the 
viscosity  matrix.  Many  physical  models  are  of  this  form,  e.g.  the  compressible  Navier- 
Stokes  equations,  the  equations  of  magnetohydrodynamics  and  nonlinear  viscoelasticity 
models.  We  are  interested  in  the  explicit,  pointwise  large  time  behavior  of  the  solution 
of  the  initial  value  problem  of  the  system.  The  solution  approaches,  time-asymptotically, 
a  linear  superposition  of  linear  heat  kernels,  generalized  Burgers  solutions,  and  diffusion 
waves  of  algebraic  types. 

Several  features  of  the  system  contribute  to  the  rich  qualitative  behavior  of  the  solution. 
The  first  basic  feature  is  the  nonlinearity  of  the  flux  function  f{u).  Physical  models  are 
symmetrizable,  which  implies  that  the  inviscid  system 

(1.2)  ut  +  f{u)x=  0 

is  completely  hyperbolic  with  real  characteristics  Xi(u)  <  \2(u)  <  •  ■  •  <  Xn(u): 

f'{u)ri(u)  =  X  i(u)ri(u),  li(u)f'(u)  =  X  *(«)/*(«), 

(1-3) 

li(u)rj(u)=6ij,  i,j  =  1,2 

The  first  important  nonlinearity  is  the  behavior  of  the  characteristic  speed  A*(u)  in  the 
wave  direction  ri{u).  An  i- field  is  genuinely  nonlinear  (g.nl.)  if  A^(w)  is  strictly  monotone, 
linearly  degenerate  (l.dg.)  if  it  is  stationary,  [La], 


(g.nl.)  VA i{u)-ri(u)  ^0, 

(l.dg.)  VA i(u)  •  ri(u)  =  0. 

Waves  pertaining  to  a  (g.nl.)  field  either  compress  or  expand.  Compression  waves  lead  to 
shock  waves.  A  (l.dg.)  field  gives  rise  to  linear  hyperbolic  waves,  the  contact  discontinuities. 
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For  the  viscous  system  (1.1),  (g.nl.)  fields  yield  nonlinear  diffusion  waves,  solutions  of 
Burgers’  equation;  while  (l.dg.)  fields  yield  the  linear  heat  kernels  as  the  time-asymptotic 
state. 

The  second  nonlinearity  is  the  coupling  of  waves  of  different  characteristic  families.  This 
is  measured  by 

Uf"{rj,rj),  i^j. 

For  physical  models,  these  quantities  are  nonzero.  This  implies  that  the  linear  and  non¬ 
linear  diffusion  waves  interact  nonlinearly  and  produce  diffusion  waves  which  decay  alge¬ 
braically  in  x  and  t. 

Another  nonlinearity  is  the  dependence  of  the  multiplicity  of  the  characteristic  values 
A i(u),  i  =  1,  2,  . . .  ,  n,  on  the  state  u.  The  compressible  Euler  equations  are  always  strictly 
hyperbolic.  But  many  other  models,  including  the  equations  of  magnetohydrodynamics, 
full  elasticity  system  and  multiphase  flows  are  nonstrictly  hyperbolic.  For  instance,  Alfven, 
slow  and  fast  acoustic  speeds  may  equal  for  the  equations  of  magnetohydrodynamics,  see 
[LLP].  When  the  inviscid  model  is  nonstrictly  hyperbolic,  system  (1.1)  possesses  time- 
asymptotic  states  of  solutions  of  generalized  Burgers  equations.  These  are  systems  with 
single  characteristic  speed  and  quadratic  nonlinear  flux. 

The  second  basic  feature  of  the  viscous  conservation  laws  (1.1)  is  that  physical  viscosity 
matrix  B(u)  is  not  diagonalizable  with  the  inviscid  system,  and  not  even  positive  defi¬ 
nite.  This  is  dictated  by  physics.  For  instance,  the  continuity  equation  is  not  dissipative. 
Physical  models,  nevertheless,  are  dissipative  when  the  relevant  dissipation  parameters  are 
present.  The  general  mathematical  hypotheses  are  to  assume  that  there  exists  a  dissipa¬ 
tive  symmetrizer  and  that  no  characteristic  direction  of  the  flux  is  in  the  null  space  of  the 
viscosity  matrix.  With  such  physical  assumptions,  the  hyperbolic  nature  of  system  (1.1)  is 
reflected  in  the  propagation  of  the  discontinuities  in  the  initial  data  into  the  solution  at  all 
time.  This  is  reflected  in  the  present  study  as  Dirac  6-functions  contained  in  the  Green’s 
function  for  the  linearized  system  of  (1.1).  We  have  explicit  expressions  for  the  character¬ 
istic  information  on  the  6-functions  and  the  rate  of  exponential  decay.  This  non-smoothing 
property  complicates  our  analysis  on  large  time  behavior.  That  the  viscosity  matrix  is  not 
diagonalizable  with  the  Jacobi  matrix  of  the  flux  means  that  there  is  a  coupling  of  the 
heat  kernels  in  the  Green’s  function.  We  show  that  the  Green’s  function  equals  that  for  a 
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diagonalizable  system  plus  faster  decaying  terms. 

The  present  paper  is  concerned  with  the  large  time  behavior  of  solutions  of  (1.1)  that 
are  perturbations  of  a  constant  state.  Without  loss  of  generality,  we  take  the  constant 
state  to  be  the  zero  state, 

(1-4)  u(x,0)  =  u0(x),  u0(±  oo)  =  0. 

Our  approach  contains  two  new  major  elements,  the  accurate  pointwise  estimate  on  the 
Green’s  function  of  the  linearized  equations 

(!-5)  t*t  +  //(0K  =  B(0)tiaMBj 

and  the  analysis  on  the  nonlinear  coupling  of  diffusion  waves.  We  now  briefly  describe  the 
setup.  Since  the  system  is  dissipative,  the  solution  is  expected  to  decay  to  the  zero  state 
in  the  sup  norm.  Expand  the  solution  around  the  zero  state,  we  have  from  (1.1)  that 

n  1  1 

uu+  Ai(0)u*+  ^2  9W0)f(°)(wjTi(°)^fcffc(0)) 

j,k=  1  -lx 

n 

(1.6)  =  ^(OJBtOXu^O)),],  +  [0(1)(M3  +  MW)]*,  t  =  1, ...  ,n, 

3= 1 
n 

u(x,t)  =  ^Ui(x,t)ri(0). 

i=  1 

The  first  step  is  to  approximate  this  by  a  simpler  system,  which  yields  explicit  solutions. 
For  strictly  hyperbolic  case,  Ai(0)  <  A2(0)  <  •  •  •  <  An(0),  we  ignore  all  the  coupling  terms 
of  different  families  and  the  higher  order  terms  beyond  the  quadratic  terms  to  obtain 

(1.7)  eit  +  \i(0)6iX  +  [L|4 (0)/"(0) (rj(0),  n( 0))6f]x  =  li(0)B(0)ri(0)eixx,  t  =  1, ...  ,71. 

This  is  the  linear  heat  equation  or  Burgers  equation.  For  nonstrictly  hyperbolic  case,  the 
same  procedure  yields  a  system  of  equations  similar  to  (1.7).  These  equations  are  called 
generalized  Burgers  equations,  where  all  the  modes  corresponding  to  the  same  characteris¬ 
tic  speed  at  u  =  0  are  coupled.  In  any  case,  the  approximate  equations,  (1.7)  or  generalized 
Burgers  equations,  possess  self-similar  solutions.  These  are  the  linear  and  nonlinear  heat 
kernels: 


(1.8) 


a  (  j.\  j.— —  / Ai(0)£ 

6i{x^t)  t  2  Xi{ 
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for  some  functions  These  solutions  represent  the  time- asymptotic  state  for  the  approx¬ 
imate  equations.  They  are  also  the  main  term  in  the  time-asymptotic  state  for  the  solution 
of  (LI).  Equation  (1.8)  represents  a  one-parameter  family  of  solutions  with  the  parameter 

/oo 

6i{x,t)  dx. 

-oo 


Choose  Ci  so  that 


f. 


v{x ,  t)  dx  =  0, 


n 

v(x,t)  =  u{x,t )  -  ^YiOM)n(o). 
1 


For  simplicity,  we  use  the  strictly  hyperbolic  case  to  illustrate  what  we  need  to  estimate 
v  and  the  structure  of  the  paper.  Let  G(x,t)  be  the  Green’s  function  associated  with  the 
linearization  of  (1.6),  and  G*{ x,  t)  be  the  Green’s  function  associated  with  the  linearization 
of  (1.7).  Clearly  G*  is  a  diagonal  matrix  whose  ith  diagonal  element  G*t  is  a  heat  kernel 
along  the  Xt(0)  direction.  We  have  from  Duhamel’s  principle  that 


(1.9) 


/oo  r°° 

G*i(x  -  y,t)vi(y,0)dy  +  /  (G  -  G*)i(x  -  y,t)u(y,  0)  dy 

-oo  J  —  oo 

ft  roo  V  _  n 

+  G*i(x  -  y,t  —  t)  ^c^  +  ^c^  +  OWM2 

Jo  J- OO  L,-^t  7  =  1 

(y,r)dydr 


+  Y,0(l)8j(ek  +  vk) 


J  y 


ft  fOO 

+  (G  —  G*)i(x  —  y,t  —  t)  [0(1)(|#|2  +  \v\2)]y(y,r)dydT 

Jo  J —oo 

ft  fOO 

+  Gi(x  —y,t  —  r)  [O(l)(|0|3  +  |f|3 

Jo  J — OO 

+  (|0|  +  M)(l<y  +  kl))]y(v>T)  dvdT< 


i  —  1, . . . ,  ti, 

n 

v(x,t)  ='y]vi(x,t)ri(  0), 

i= 1 


where  (G  -  G*)*  and  G*  denote  the  ith  row  of  G  -  G*  and  G  respectively,  Cy ,  i,  j  =  1, 
...  ,  n,  are  constants,  and 


u  =  (ui, . . .  6  =  (0\,. . .  ,6ny> 
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A  more  precise  expression  of  v  in  the  general  case  (not  necessary  to  be  strictly  hyperbolic)  is 
derived  in  Section  8,  while  the  exact  formulation  of  our  main  theorems  about  the  estimates 
of  v  is  contained  in  Section  2. 

To  make  use  of  (1.9)  we  need  to  have  accurate  estimates  of  the  Green’s  function  G  and 
the  convolutions  of  the  heat  kernel  with  the  sources  J2  &j  and  the  nonlinear  term  v2. 
The  latter  is  the  study  of  nonlinear  coupling  of  diffusion  waves.  The  diffusion  waves  are 
either  of  the  linear  or  nonlinear  heat  kernel  types  (1.8)  or  of  algebraic  types 

^3/2 (^?  Aj(0))  =  [{x  —  Aj(0)(t  +  l))2  + 1  +  1]  3//4, 

^i(O))  =  [\x  —  Aj(0)(i  +  1)|3  +  (t  +  l)2]  1/2, 

(</)i'03/4)(x?  t]  A*(0))  =  (|x  —  A*(0)(£  +  1)|  +  t  +  1)  1//2 

x  [(ar  —  A*(0)(£  +  l))2  +t  +  l]”3/8. 

In  Section  3  we  study  the  dissipation  of  these  waves.  Our  analysis  is  motivated  by  and 
greatly  generalize  that  of  [Liu],  where  diffusion  waves  of  algebraic  types  were  identified. 

For  a  diagonalizable  system,  the  Green’s  function  consists  of  heat  kernels.  The  time- 
asymptotic  analysis  for  this  simplified  case  is  carried  out  in  Section  4. 

Physical  viscosity  matrix  is  non-positive  definite,  and,  as  a  result,  the  Green’s  function 
contains  Dirac  ^-functions.  Moreover,  the  system  is  nondiagonalizable  and  so  the  Green’s 
function  differs  from  that  for  a  diagonalizable  system  also  by  a  term  which  decays  at  a 
faster  rate.  These  phenomena  were  found  in  [Ze]  for  the  isentropic  compressible  Navier- 
Stokes  equations.  In  Section  5,  we  study  the  same  problem  but  using  a  simpler  method. 
This  new  method  can  be  extended  to  general  systems;  see  Section  6. 

In  our  study  of  dissipation  nature  of  the  linearization  of  (1.6),  or  equivalently  the  lin¬ 
ear  system  (1.5),  through  the  pointwise  estimates  of  its  Green’s  function,  we  make  uses 
of  Kato’s  perturbation  theory,  [Kt],  to  analyze  the  singularities  of  the  spectrum  of  the 
matrix  in  Fourier  variables.  In  consistence  with  physical  models,  we  assume  that  (1.1)  is 
symmetrizable.  That  is,  it  is  endowed  with  an  entropy  function.  This  removes  the  pos¬ 
sible  singularity  of  the  spectrum  at  the  origin  or  at  the  infinity.  We  then  show  that  the 
Green’s  function  G  for  (1.5)  contains  G* ,  the  heat  kernels  corresponding  to  the  associated 
diagonalizable  system.  Again,  motivated  by  physical  consideration,  we  assume  that  the 
null  space  of  the  viscosity  matrix  B  does  not  contain  the  inviscid  characteristic  directions. 
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This  guarantees  the  dissipativeness  of  the  system.  We  show  that  the  remainder  G  —  G* 
decays  at  a  faster  rate  than  that  for  G*.  The  remainder  G  —  G*  also  contains  exponen¬ 
tially  decaying  Dirac  6-functions  due  to  the  hyperbolic  nature  of  our  system.  In  Section 
6,  we  also  apply  this  general  theory  to  the  compressible  Navier- Stokes  equations  and  the 
equations  of  magnetohydrodynamics.  We  calculate  out  explicitly  the  heat  kernels  as  well 
as  the  Dirac  6-functions  contained  in  G  for  these  equations. 

For  the  isentropic  compressible  Navier-Stokes  equations,  it  has  been  shown  that  dis¬ 
continuities  in  the  initial  data  propagate  into  the  flow  along  the  particle  path,  [H].  This 
corresponds  to  the  6-function  in  the  Green’s  function.  For  the  full  compressible  Navier- 
Stokes  equations,  we  also  find  that  there  is  a  6-function  along  the  particle  path  when  both 
viscosity  and  heat  conductivity  coefficients  are  positive.  When  the  viscosity  coefficient  is 
zero,  this  6- function  splits  into  two  6-functions  about,  but  not  on,  the  particle  path.  For 
the  equations  of  magnetohydrodynamics,  there  may  exist  three  6-functions. 

Section  7  concerns  the  energy  estimate  due  to  Kawashima  [Ka] .  Since  the  presentation 
is  short,  we  present  it  here  for  completeness.  For  the  estimate  we  need  the  additional 
hypothesis  of  invariance  of  the  null  space  of  the  viscosity  matrix  with  respect  to  the  sym- 
metrizer;  see  Assumption  2.2  of  Section  2.  The  energy  estimate  is  needed  to  close  the 
stability  analysis  for  the  general  system,  which  is  carried  out  in  Section  8. 

Finally,  we  verify  the  aforementioned  dissipativeness  hypotheses  for  the  compressible 
Navier-Stokes  equations  and  the  equations  of  magnetohydrodynamics  in  Section  9. 

IP  optimal  convergence  rate  of  v  in  (1.9)  has  been  obtained  by  Zeng,  [Z],  which  general¬ 
ized  her  previous  work  on  p-system,  [Ze] .  Here  our  new  pointwise  approach  yields  explicit 
information  on  the  propagation  of  diffusion  waves.  The  analysis  is  effective  in  studying  the 
coupling  of  these  waves.  Many  physical  models  are  only  partially  dissipative.  This  is  so, 
for  instance,  for  the  compressible  Navier-Stokes  equations  with  positive  viscosity  and  zero 
heat  conductivity,  and  also  for  the  compressible  Euler  equations  with  thermal  nonequilib¬ 
rium.  For  such  a  system,  there  are  diffusive  and  non-dissipative  waves.  We  believe  that 
our  analysis  can  be  generalized  to  study  the  behavior  of  these  waves.  This  is,  however,  left 
to  the  future. 


2.  Viscous  Conservation  Laws 


Consider  the  Cauchy  problem  of  a  general  system  of  viscous  conservation  laws: 

(2>1)  ut  +  f(u)x  =  (B(u)ux)x,  -oo  <  x  <  oo,  t  >  0, 

(2-2)  u(a:,0)  =  uQ(x),  -oo  <  x  <  oo, 

where  u  =  u(x,  t)  is  an  n-vector,  f(u)  is  a  smooth  n-vector  valued  function,  B(u)  a  smooth 
n  x  n  matrix,  and  Uo(a;)  a  given  n-vector  valued  function.  We  assume  that  Uq  is  a  small 
perturbation  of  a  constant  state.  Without  loss  of  generality,  we  assume  that  the  constant 
state  is  zero.  To  assure  that  system  (2.1)  is  dissipative  we  make  the  following  three  basic 
assumptions: 

Assumption  2.1.  System  (2.1)  has  a  strictly  convex  entropy  U  ([Sm],  [K] ) . 

Assumption  2.1  means  that  there  exist  a  pair  of  smooth  functions  U(u)  and  F(u)  (called 
entropy  pair),  such  that  U  is  strictly  convex,  ( VC/)/'  =  VF,  and  ( V2U)B  is  symmetric 
and  semi-positive  definite,  where  /'  is  the  Jacobi  matrix  of  /,  and  V2U  is  the  Hessian  of 
U.  We  set  Aq(u)  =  V2U{u).  Clearly  A0  is  symmetric  and  positive  definite. 

Assumption  2.2.  For  small  u,  B(u)  ^  0.  There  exists  a  smooth  one-to-one  mapping 
u  —  fo(u),  /o(0)  =  0,  such  that  the  null  space  J\f  of  B(u)  =  F(/o(fi))/o(fi)  is  independent 
of  fi.  Moreover,  M1-  is  invariant  under  f^uf  A0(f0(u)),  and  B(u)  maps  Rn  into  J\fL ,  where 
is  the  orthogonal  completement  of  J\f. 

Assumption  2.3.  Any  eigenvector  of  /;(0)  is  not  in  the  null  space  of  B( 0). 

We  notice  that  several  interesting  physical  models,  e.g.  the  Navier-Stokes  equations  and 
the  equations  of  magnetohydrodynamics,  satisfy  Assumptions  2. 1-2.3  (c.f.  Section  9).  We 
also  notice  that  if  B(u)  is  nonsingular,  Assumptions  2.2  and  2.3  are  automatically  satisfied. 

Under  Assumption  2.1,  it  follows  that  A0f  is  symmetric  [FL].  Since  A0  is  symmet¬ 
ric  and  positive  definite,  there  exists  a  symmetric  positive  definite  matrix  Al  such  that 

ill  ^ 

Aq  (Aq  )  .  Thus  Aq  f  (A~ )  is  symmetric.  Consequently,  all  the  eigenvalues  of 

i  i_  ii 

*S  f'(A$ )  1  are  real,  and  A£  f'(A§  has  a  complete  set  of  eigenvectors.  Notice  that  if 

i  i  l 

ASfiAS)-1  has  an  eigenvalue  A  and  a  corresponding  eigenvector  r,  f'(A^)~1r  =  Ar, 
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then  f'(A§  )~lr  =  A (a\  )"1r,  i.e.,  A  is  an  eigenvalue  of  /'  with  {A§ )  V  as  a  corresponding 
eigenvector.  For  this  reason  f  has  real  eigenvalues  Xi(u)  <  ^(u)  <  ■■  •  <  An(u)  and  a 
complete  set  of  eigenvectors: 

f'(u)fi(u)  =  A  i(u)fi(u), 

(2.3)  Ii(u)f'(u)  =  A  *(u)/*(u), 

li(u)fj(u)  =  6ij ,  i,j,=  1,2, ...  ,n. 

In  other  words,  the  inviscid  system 

(2.4)  Ut  4-  f{u)x  ~  0 

corresponding  to  (2.1)  is  completely  hyperbolic. 

We  are  interested  in  the  perturbation  of  the  zero  state.  Hence  the  eigenvalues  and 
eigenvectors  of  /'( 0)  play  an  important  role.  We  use  special  notations  to  denote  them.  Let 
Ai  <  A2  <  *  •  *  <  As  be  the  distinct  eigenvalues  of  /'( 0)  with  multiplicity  mi,  m2,  ...  ,  ms, 
mi  _j_  m2  _l - \-ms—n.  Let  the  left  and  right  eigenvectors  associated  with  A*  be  kj  and 

T  ij ,  j  —  1 ,  ...  ,  771  j . 

/  (d)T’U  —  XiTij, 

Ujf  (0)  =  Xikj, 

(2.5) 

i,i'  =  l,...,s,  j  =  l,...,mi,  f  =  l,...,m*'. 

Introduce  the  notations 

/ l"  \ 

(2.6)  h  =  I  ;  J  )  7*i  —  (7**1,  .  •  *  ,  7 ^  ~  1,  ■  •  •  ,  5. 

\kmj 

Clearly, 

(2.7)  —  ImiXmn  hrj  ~  0 ^  7^  3* 

Decompose  the  solution  u  to  (2.1),  (2.2)  in  the  right  eigenvector  directions, 

s  rrti  s 

u(x,t)  = 

•=1  j= 1  i=l 


(2.8) 
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where 

(2.9)  Ui{x,  t)  (Uili  •  •  •  ,  Uirrii  )  (*^>  ~  ^w((C,  £),  2  =  1,  .  .  .  ,  S. 

Our  main  result  says  that  m  =  6t  +  Vi,  where  0t  is  a  self-similar  solution  to  the  Burgers 
equation,  the  heat  equation  ,  or  the  generalized  Burgers  equation,  while  the  remainder  V{ 
decays  more  rapidly  in  t. 

To  give  a  heuristic  explanation  of  the  construction  of  the  ansatz  Oi ,  we  first  consider 
the  equation  satisfied  by  v*.  With  (2.9)  and  (2.8),  multiply  (2.1)  by  lu  and  take  Taylor 
expansion  about  zero  to  obtain 

(2.10)  uit  +  \uix  +  \lif\0)(y^jrjuj,y^rjuj)x 

3  3 

=  hB( 0)  rjUjxx  +  (0(1)0^11^1  +  M3))*,  *  =  1,  •  •  • ,  s. 

3 

Each  of  the  above  equations  is  associated  with  an  eigenvalue  A*.  Neglecting  the  effects  of 
the  other  families  and  the  higher  order  terms,  we  have  the  following  approximate  equations: 

(2.11)  0it  +  \{6ix  +  ri6i)x  -  hB(0)ri0ixx,  i  =  1, . . . ,  s, 

where  Oi  is  an  m*- vector.  The  ansatz  Oi  of  Ui  is  defined  as  the  self-similar  solution  to  (2.11) 
carrying  the  same  mass  as  u*: 

/OO  /»CO 

Ot(x,t)dx~  /  liU0(x)dx  =  Si. 

■oo  J  —  OO 

We  call  i  =  1, . . .  ,s,  the  diffusion  waves  of  (2.1),  (2.2).  The  following  lemma  shows 
that  (2.11)  is  uniformly  parabolic,  though  the  full  system  (2.1)  may  be  only  hyperbolic- 
parabolic. 

Lemma  2.1.  [SK]  Under  Assumptions  2.1  and  2.3,  each  j?(0)r*  is  similar  to  an  mi  x  m* 
symmetric  positive  definite  matrix,  i  =  1, . . . ,  s. 

Proof.  Let  IU  and  rL,  i  =  1, . . . ,  s,  j  =  1, . . . ,  mi  be  other  sets  of  left  and  right  eigenvectors, 
respectively,  satisfying  (2.5).  Using  the  notations  in  (2.6)  it  is  clear  that 
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where  c*  and  d%  are  mt  x  m%  matrices.  Since  both  pairs  and  {£',r'}  satisfy  (2.7), 

we  have  di  =  c^1.  Thus 

l'B(0)r'  =  CiliB(Q)ri<Z\ 

i.e.,  /'£(0)r'  is  similar  to  kB( 0)rj.  The  lemma  is  proved  if  we  can  construct  kj  and  r^-, 
i  =  1, . . . ,  s,  j  =  1, . . . ,  rrii ,  such  that  each  /jl?(0)r*  is  symmetric  and  positive  definite. 

Under  Assumption  2.1,  both  AqB  and  Aq f  are  symmetric,  and  AqB  is  semi-positive 
definite,  where  A0  =  V2U.  To  simplify  notations  for  the  moment  we  set  Aq  =  A0(0), 
B  =  B( 0),  etc.  As  we  have  pointed  out  before,  {r^}  are  right  eigenvectors  of  /'  if  and 

1  X  1  1 

only  if  {Aq  r^}  are  right  eigenvectors  of  A§ff(A£)~l,  where  A£  is  a  symmetric  positive 
definite  matrix  such  that  A0  =  (A£)2.  Similarly,  {kj}  are  left  eigenvectors  of  /'  if  and 
only  if  {kjiAl)-1}  are  left  eigenvectors  of  A§f'{A§)~1.  We  can  choose  {ri:j}  such  that 
{A^rij}  is  an  orthonormal  set  since  Agf^Ag)-1  is  symmetric.  The  left  eigenvectors  of 

i  i 

Aq  /'(Aq)-1  are  then  chosen  as 

m4)_1  =  (4^)*, 

or 

hj  Oj 

i  =  1,. . .  ,s,  j  =  1, .. .  ,mi.  Clearly  {kj}  and  {r^}  satisfy  (2.5).  Let  R=(ru.. .  ,ra).  then 
each  UBri  is  an  x  m.t  principal  submatrix  of  RtA0BR.  Since  R1AqBR  is  symmetric 
and  semi-positive  definite,  so  is  UBri.  If  liBr%  were  not  positive  definite,  there  would  be  a 
w  e  Rm* ,  w  ^  0,  such  that 


wtUBriW  —  (riw)1  AoBriW  —  0. 


Thus 


A0BriW  =  0, 


B(riw)  =  0. 


Here  uw  Y  0  would  be  an  eigenvector  of  f  corresponding  to  A*.  This  contradicts  Assump- 

□ 


tion  2.3. 


QUASILINEAR  HYPERBOLIC-PARABOLIC  SYSTEMS 


11 


Consider  the  special  case  that  A*  is  simple,  m*  =  1.  Then  0*  is  a  scalar  and  (2.11) 
becomes 

(2-13)  On  +  \6ix  +  Cii{&i)x  —  HiOixx, 


where 

(2-14)  <ki  =  l/i/"(0)(ri,rj),  m=liB(0)ri. 

By  Lemma  2.1,  /i,  >  0.  Equation  (2.13)  is  the  heat  equation  if  cu  =  0.  The  corresponding 
diffusion  wave  is  then  the  heat  kernel: 


Oi(x,t)  = 


+  1) 


4^i(i+l) 


with  6i  given  by  (2.12).  If  cu  ^  0,  equation  (2.13)  is  the  Burgers  equation,  which  can  be 
solved  by  Hopf-Cole  transformation  [Ho],  [Co], 


C-ii 

On  4"  XiOix  —  (liOixx. 


The  diffusion  wave  is 


(2T6)  0i(Xyt) 


-l 


In  (2.15)  and  (2.16)  we  have  avoided  the  singularity  of  the  heat  kernel  at  t  =  0  by  replacing 
t  by  t+ 1.  Notice  from  differentiating  (2.3)  that  kfn{0){ru  n)  =  VA*(0)rv  Thus  pertaining 
to  a  linearly  degenerate  field  we  have  a  linear  diffusion  wave  (2.15),  and  to  a  genuinely 
nonlinear  field  we  have  a  nonlinear  diffusion  wave  (2.16),  [La]. 

When  Xi  is  not  simple,  mi  >  1,  system  (2.11)  is  called  the  generalized  Burgers  equations 
[Ch].  The  self-similar  solution  of  (2.11),  (2.12)  is  quite  similar  to  the  nonlinear  diffusion 
wave  (2.16).  Precisely,  we  have  the  following  lemma  due  to  Chern. 
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Lemma  2.2.  [Ch]  Under  Assumptions  2.1  and  2.3,  for  small  Si,  (2.11),  (2.12)  has  a 
unique  self-similar  solution  of  the  form  ,  with  =  Xi(-°°)  =  °* 

Furthermore,  \i  h as  & e  property 

2 

(2.17)  Xi(y)  =  e~,“iXi(y), 

where  m  is  the  maximum  eigenvalue  ofliB{0)rit  and  Xi  and  all  its  derivatives  are  uniformly 
bounded  by  C\St\,  C  >  0  is  a  constant. 

The  proof  of  Lemma  2.2  will  be  given  in  the  proof  of  Lemma  3.1. 

Again  to  avoid  the  singularity  at  t  =  0,  the  diffusion  wave  pertaining  to  the  z-family  is 
then  defined  as 

.  x  1  x  —  \i(t  1) . 

(2-18)  6i{-x’  ^  =  7f+TXi(  vT+ 1 

where  ^Xi(Xy^*)  is  the  self-similar  solution  in  Lemma  2.2. 

Introduce  the  following  notations, 

4>*{x,  t;  A)  =  (|s  -  A(t  +  1)|  + 1  +  1)-Q/2, 

(2.19)  ipa(x,t;  A)  =  [(a:  -  \(t  +  l))2  +  t  +  1]  “^2, 

A)  =  [|*  —  A(t  +  1)|3  +  (<  +  1)2]-1/2. 

We  call  them  generalized  diffusion  waves. 

The  main  result  of  this  paper  is  the  following  pointwise  estimate  of  the  solution  to  (2.1), 

(2.2). 

Theorem  2.3.  Suppose  that  system  (2.1)  satisfies  Assumptions  2. 1-2.3,  and  that  the 
initial  data  satisfy  the  following  smoothness  and  smallness  assumptions : 


(2.20) 

«o  €  H8( R), 

(2.21) 

v-(x)=  /  u0(y)dy  £  L1 

J  —  oo 

(K  ),  v+(x)  =  f  u0(y)dy  6  Ll{R+), 
J  X 

(2.22) 

«o(*)  =  0(1)(1  +  M)~8/4. 

«i(*)=o(i)(i+wr,/4, 

(2.23) 

«±(*)  =  0(1)(1  +  W)-1, 

II«o||h8  +  ll«-llt>(-oo,0)  +  1 

l|t'+||i‘(0,oo) 

+  sup{(l  +  |x|)5/4(|u0(a:)l  +  l«oWI)} 

x€R 

4-  sup  {(1  +  x)(|u_(-x)|  +  |v+(a;)|)}  =  6*  <  1. 

0<x<oo 


(2.24) 
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Then  for  all  x  £  R,  t  >  0,  the  solution  u  to  (2.1),  (2.2)  has  the  property 


u(x,t)  =  22ri0i(x,t)  +  v(x,t), 

i=  1 

(2.25)  v(x,t)  =  0(l)6*$(x,t),  vx(x,t)=0(l)6*(t  +  l)-1'2<Z(x,t), 

S 

*CM)  =  ^[^3/2(2,  t\ \i)  +  fafx,  t;  Xi)tp3/4(x,  t;  A^)], 

i=  1 


where  9it  i  —  1, . . . ,  s,  is  the  self-similar  solution  to  (2.11),  (2.12),  given  by  (2.15),  (2.16) 
or  (2.18);  A,,  r,  and  l,  ,  i  =  1, . . . ,  s,  are  eigenvalues  and  matrices  consisting  of  right  and 
left  eigenvectors  of  f'( 0),  given  by  (2.5)  and  (2.6);  ipa  and  cj>a  are  defined  in  (2.19). 

Remark  2.f.  From  (2.25),  (2.19)  we  see  that  for  k  =  0, 1, 


(2.26) 


dk 

dxk 


v{-,t) 


dk  ,  . 
dk  ,  s 


L1 


L 2 


L°° 


0(l)6*(t  +  l)~i~§, 


All  these  rates  are  optimal. 


Remark  2.5.  If  the  left  eigenspace  of  B  associated  with  the  zero  eigenvalue  is  independent 
of  u  (including  B  being  nonsingular),  to  obtain  the  pointwise  estimate  for  v(x,t)  in  (2.25), 
it  is  sufficient  to  require  u0  €  H4(R)  rather  than  u0  G  H8( R). 

From  (2.25)  we  see  that  the  remainder  v(x ,  t)  has  decay  rate  |;r|-5/4  as  x  — ►  oo,  which  is 
consistent  with  the  decay  rate  of  the  initial  data  in  (2.22),  and  can  be  improved  if  the  decay 
rate  of  the  initial  data  is  improved.  For  example,  if  uq  and  uf0  decay  exponentially,  so  do  v 
and  vx  in  x.  From  (2.25)  and  (2.19)  we  also  see  that  v  and  vx  have  decay  rates  (t  + 1)-3/4 
and  (t  + 1)-5/4,  respectively,  along  characteristic  directions  dx/dt  =  A;,  i  =  1, . . . ,  s.  These 
two  rates  are  optimal  and  cannot  be  improved  by  the  rapid  decay  of  the  initial  data.  Away 
from  the  characteristic  directions  v  and  vx  have  decay  rates  (t  +  l)-5/4  and  (t  +  1)~7/4 
respectively,  which  can  be  improved  up  to  (t  +  1)~3/2  and  ( t  +  1)“2.  More  precisely,  we 
have  the  following  theorem. 
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Theorem  2.6.  Suppose  that  system  (2.1)  satisfies  Assumptions  2. 1-2.3,  and  that  the 
initial  data  satisfy  the  following  smoothness  and  smallness  assumptions: 

«0  €  #9(R), 

u0(x)  =  0(1)(1  +  M)-3A  u'0(x)  =  0(  1)(1  +  \x\)-V\ 

/. x  r  oo 

uo(y)dy  e  L1^),  v+(x)  =  j  iio(y)dy  e  L1(R+), 

|l«o|l//9  +  11^—  oo,0)  +  ll^-f  llz/1(0,oo) 

+  sup{(l  +  |x|)3/2(M*)l  +  K(x)l)}  s  5*  «  1. 

Then  for  all  x  ER,  t  >  0,  the  solution  u  to  (2.1),  (2.2)  has  the  property: 

s 

u(x,t)  =  +  Vi(x,t)}, 

1=1 

(2.28)  «,(*,*)  =  0(l)«**i(*,t),  vix(x,t)  =  0(1)6"  (t+ 

$i(x,t)  =  xp3/2(x,t;\i)  +  Aj), 

j^i 

where  $i,  i  =  1, . . .  ,s,  is  the  self-similar  solution  to  (2.11),  (2.12),  given  by  (2.15),  (2.16) 
or  (2.18);  A*,  r*  and  li}  i  =  1, . . . ,  s,  are  eigenvalues  and  matrices  consisting  of  right  and 
left  eigenvectors  of  f'(  0),  given  by  (2.5)  and  (2.6);  ipa  and  (j)  are  defined  in  (2.19). 

Remark  2.7.  The  decay  rates  (t  +  l)"3/2  and  ( t  +  l)-2  for  vz  and  vix ,  respectively,  away 
from  the  characteristic  directions  are  optimal.  They  are  due  to  the  nonlinear  coupling  of 
different  families  and  can  not  be  improved  even  if  the  initial  data  are  smoother  and  decay 
faster. 

Remark  2.8.  Similar  to  Remark  2.5  to  Theorem  2.3,  if  the  left  eigenspace  of  B  associated 
with  the  zero  eigenvalue  is  independent  of  u ,  then  the  stringent  restriction  uq  E  H9( K)  in 
Theorem  2.6  is  relaxed  as  u0  E  H4(R)  to  obtain  the  pointwise  estimate  for  in  (2.28), 
and  uq  E  H6(M)  to  obtain  the  pointwise  estimate  for  vix. 

If  B  is  a  positive  definite  constant  matrix,  e.g.,  B  =  I,  we  further  have  the  following 
theorem. 

Theorem  2.9.  Suppose  that  B  is  a  constant  matrix,  there  exists  a  symmetric  positive 
definite  matrix  A0,  such  that  Aof'(0)  is  symmetric,  and  A0B  is  symmetric  positive  definite, 
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and  the  initial  data  satisfy  the  following  smallness  assumptions: 


(2.29)  u0(x)  =  0(1)(1  +  |x|)-3/2,  u'0  e  L°°( K), 

/x  poo 

u0(y)dy  €  L1^-),  v+(x)=  u0{y)dyeL1(i 

-OO  Jx 

(2.31)  ||MolU~  +  |b-IUi(-oo,0)  +  lk+ 1|  1,1(0,00) 

+  sup{(l  +  |x|)3/2|u0(x)|}  =  <5*  <  1. 
xeR 

Then  for  all  x  G  M,  t  >  0,  the  solution  u  to  (2.1),  (2.2)  has  the  property: 


u(x,t)  =  Y2ri[Si(x,t)  +Vi(a;,t)], 

i= 1 

Vi(x,  t)  =  o(i)«5*  w3/2(x,  t-,\i)  +  J2  &*.  t;  A,-)] 
Vix(x,t)  =  0(1)6*  + 


where  $it  i  =  1, . . . ,  s,  is  the  self-similar  solution  to  (2.11),  (2.12),  given  by  (2.15),  (2.16) 
or  (2.18);  \,  r*  and  li,  i  =  1, . . .  ,s,  are  eigenvalues  and  matrices  consisting  of  right  and 
left  eigenvectors  of  f'(  0),  given  by  (2.5)  and  (2.6);  V>Q  and  jj)  are  defined  in  (2.19). 

Throughout  this  paper  we  use  C  to  denote  a  universal  positive  constant. 


3.  Diffusion  Waves 


We  now  discuss  some  properties  of  the  diffusion  waves  defined  in  Section  2,  including 
their  convolution  with  the  heat  kernel.  We  also  discuss  the  convolution  of  the  generalized 
diffusion  waves  with  the  heat  kernel.  These  properties  will  play  an  important  role  when 
assessing  the  effect  of  the  nonlinear  error  terms  in  the  a  priori  estimates  in  Sections  4 
and  8. 

Lemma  3.1.  Letdit  i  =  1,  . . .  ,  s,  be  the  diffusion  waves  given  by  (2.15),  (2.16)  or  (2.18). 
Then  for  all  -oo  <  x  <  oo,  t  >  0,  we  have 

ei(x,t)=0(l)\6i\(t+l)-’e-vl 

6ix(x,t)  =  0(l)|«4|(t  +  ir'dl/il  + 

6a(x,t)  =  0(l)|tfi|(i  +  l)_1(?/»2  +  Ittl  +  l)e_yS 

(31)  [8tt  +  Xi9ix)(x,t)  =  0(l)|5i|(t  +  1  )~Hvi  +  IVil  +  l)e“yS 

eixx(x,t)  =  o(i)|«i|(t  +  M  +  i)e_y‘ . 

_ x  —  A  i(t  +  1) 

s/ijjfft  +  1) 

Proof.  By  straightforward  computation,  we  obtain  (3.1)  from  (2.15),  (2.16)  and  Lemma 
2.2.  We  now  prove  Lemma  2.2.  To  simplify  the  notation  we  suppress  the  subscript  i. 
Substituting  8{x,t)  =  j-t x(^)  in  (2-n).  we  have 

-\(yxY  +  \Q(x)'  =  dx", 

where  y  is  the  independent  variable  of  x>  and 

Q{9)=lf"(Q)(re,rd),  D  =  lB(0)r. 

Integrating  this  equation  and  using  the  assumption 

lim  yx(y)  = 

y-y-OO 

we  obtain 

(3.2)  DX'  =  -\yx+\Q(x)- 
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By  Lemma  2.1,  D  has  eigenvalues  pk  >  0,  k  =  1,  . . .  ,  m,  associated  with  right  and  left 
eigenvectors  pk  and  qk ,  satisfying 

QkPi  =  Skh  k,  l  =  1, . . .  ,m. 

m 

Set  p  —  max{/^i, . . . ,  pm}.  Substituting  (2.17)  in  (3.2),  expanding  \  XkPk,  and 

A;— 1 

noticing  that  Q  is  a  quadratic  form,  we  find 

(3-3)  Xk  =  “  7 -)XkJre~^Qk{x),  k  = 

z  pk  p 

where  Qk  =  ^ j^QkQ-  We  first  prove  that  (3.3)  (and  hence  (3.2))  has  a  unique  solution 
satisfying  the  initial  condition 

m 

(3.4)  x(0)  =  w  -  ^2  WkPk 

k= 1 

if  |w|  is  small.  With  the  abbreviation  ak  =  \{l/pk  -  1/p)  >  0,  (3.3),  together  with  (3.4), 
is  equivalent  to 

m 

.  x{y)  =  ^2(e-aky2wk  + 

=  T(x)(v)- 

Clearly  T  is  an  operator  from  Cb  (all  Rm-valued  bounded  continuous  functions  on  R  with 
the  sup  norm  ||  -  ||oo)  to  itself,  and 

ITOIloo  <c(H  +  llxllL)i 

II T(X)  -  T(x)||oo  <  c\\x  +  xlloollx  -  Xlloo 

for  some  constant  C  >  0.  Thus  for  small  e,  if  fy|  <  e,  T  is  a  contraction  mapping  on  the 
subspace  [\  £  Cb  |  ||x||oo  <  2Cs}.  Here  C  is  the  same  constant  as  in  (3.6).  The  fixed-point 
theorem  then  implies  that  (3.5)  has  a  unique  solution  x,  which  satisfies  \\x\\oo  <  2 Ce  by 
(3.6).  It  is  easy  to  check  by  (3.3)  and  (3.5)  that  all  the  derivatives  of  x  are  also  uniformly 
bounded  by  0(l)e. 

We  have  proved  that  for  small  w,  (3.2)  has  a  unique  solution  x  with  x(0)  =  w.  Notice 
that  the  mass  carried  by  the  self-similar  solution  9  is  time  invariant: 

Lx{y)d“' 


'/4,1Qkmv))dr})p, 
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We  consider  the  mapping 
(3.7) 


w 


x(y)dy- 


From  (3.2)  the  derivative  of  x{y)  —  x(P5 w)  whh  respect  to  w  for  w  —  0  is 


(Pi  *  •  •  Pin) 


e 


Thus  the  derivative  of  mapping  (3.7)  for  w  =  0  is  similar  to 


2v/7rdiag(v^1,...,v7Im), 

which  is  positive  definite.  The  inverse  function  theorem  then  implies  that  for  a  given  mass 
6  which  is  sufficiently  small,  there  exists  a  unique  self-similar  solution  of  (2.11),  (2.12). 
Moreover,  x  and  all  its  derivatives  are  uniformly  bounded  by  0(1) |6|  since  |u>|  =  0(l)|<5|. 

□ 


Next  we  consider  the  convolution  of  diffusion  waves.  Without  ambiguity  of  notation, 
we  denote  the  diffusion  waves  in  the  direction  A  by 

,2 


(3.8) 


x  a  (x-A(t+l))- 

ea{x,t\X,n)  =  (t  +  l)  >e  , 


where  a  >  0,  n  >  0,  and  A  are  constants.  Notice  that  by  completing  the  square,  we  have 

_  (a  -  y  -  A(t  -  s))2  _  (y  -  A'(s  +  l))2 

/j,{s  +  1) 


(3.9) 


K*  ~ s ) 

t  +  1 


v- 


(a  +  l)(s-(A-A')(t-*)) 


fj,(t  —  s)(s  +  1) 

(x  -  A(t  —  s)  -  X'(s  +  l))2 

+  1) 


t  T  1 


Lemma  3.2.  Let  a>0,  /3>0,  (i>0  and  X  be  constants.  Then  for  all  —  oo  <  x  <  oo, 
t  >  0,  we  have 

nOO  —  s))^ 

(t  -  s)_1  (t  +  1  -  sr‘e  Op{y,  s;  A,  p)  dyds 

-oo 

f  0(l)ftT(*,i;A,M),  tfP*  3 

\  O(l)0y(x,t;\,fi)\og(t  +  2),  if  0  =  3, 
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where  7  —  a  +  min(/3, 3)  -  1/ 

*£  roo 


(3.11)  [  [  + 

J t/2  J  —  00 

-{ 


■s)y 


QffiVt  si  X,/j,)dyds 


O(l)0^{x,t;\,n),  if  a  7^1 

O(X)0i{x,  t;  A,  n)  log(t  +  2),  2/0=1, 

where  7  =  min(a,  1)  +  /3  —  1. 

Proof.  Denote  the  left-hand  side  of  (3.10)  as  /.  By  (3.8)  and  (3.9)  we  have 

rt/2 

I=s/^i  /  (t  -  »)-*(«  + 1  -  a)-»(J  +  l)-*+*(f  +  lrie-^&tP  <fe 
Jo 

rt/2 

=  O(l)0a+2(z,t;  A,p)  /  (5  +  l)“2  +  5  rfs. 

Jo 

Evaluating  the  last  integral,  we  obtain  (3.10).  Equation  (3.11)  can  be  proved  in  a  similar 
way,  where  we  partition  the  interval  of  integration  [t/2,  t]  into  [t/2,  t- 1]  and  [t- 1,  t].  □ 

Denote  the  characteristic  function  of  a  set  V  as  char{£>}. 

Lemma  3.3.  Let  the  constants  a  >  0,  0  >  1,  p  >  0,  and  A  ^  A'.  Then  for  any  fixed 

e  >0,  K  >  |A  —  A'|,  and  all  —00  <  x  <  00 ,  t  >  0,  we  have 

(3.12) 

ft  foo  ,  .a  (i-y-A(t-s))2 

/  /  (^  _  5)  (t  +  l-5)  2e  M(t-S)  ^(y,  s;  A1 ,  fijdyds 

JO  J  —  00 

=  0(l)[61(x,t-,\,n  +  s)  +  81{x,t\\',ii  +  E)  +  \x-  X(t  +  l)|_2fi|a;  -  \'(t  +  l)!-2^ 

•  char{min(A,  X')(t  +  1)  +  Kjt+1  <x<  max(A,  A')(t  +  1)  -  KVt+l}] 

+  f  0(l)6^(x,t;\,n  +  e)log(J  +  1),  if /3  =  3 
l  0,  otherwise 

+  f  O(l)01(xi  t;  A',  fi  +  e)  log(t  +  1),  if  a  =  1 
l  0,  otherwise , 

wftere  7  =  min(a,  1)  +  min(/3, 3)  -  1. 

Proo/.  Denote  the  left-hand  side  of  (3.12)  as  I.  Notice  that  for  fixed  a  and  0,  I  = 
I(x ,  t;  A,  A',  /i).  By  change  of  variable,  we  have 

/(«,  u  A,  V, .)  =  |A  -  V|  (!  o,l,  . 


A'- A 
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Thus  it  is  sufficient  to  prove  the  lemma  for  A  =  0,  A'  =  1.  By  (3.8)  and  (3.9), 


(3.13)  I  =  VnL  [i(t-s)-Ht  +  l-s)-%{s  +  l)-f+i(tH-  l)-i 

Jo 


(x-(a+l))2 

e  ds. 


Case  1.  x  <  0.  Then 
I  = 


/*V2  n  i  x2 

0(l)(t  +  l)_1_f  /  (s  +  l)~2  +  2e  M(*+Dds 
jo 

pt 

+  0(l)(t  +  i)_5  /  (t-s)-i(t  +  i-«)-« 
Jt/2 


t+i 


=  O{l)6Jx,t;0,n)  + 


e  ds 

O(l)07(a:,t;O,^)  log(t  +  1),  if  P  =  3 
0,  otherwise, 

which  is  bounded  by  the  right-hand  side  of  (3.12). 

Case  2.  0  <  x  <  Ky/t  +  1.  If  t  >  (4i^)2,  then  by  (3.13), 

P t/2 

/  =  0(l)(i+l)“1_*  /  (s  +  l)_2+ids 
Jo 

(3.14)  +  0(l)(t  +  l)“*  /  (t-5)"i(t  +  l-5)-^dse“^ 

Jt/2 

=  0(i)(.  +  ir*  +  {o(1)(,  +  1)'il0‘(‘  +  1)’  “?  =  3 

L  0,  otherwise. 

It  is  easy  to  check  that  (3.14)  is  also  true  for  t  <  (4AT)2.  The  right-hand  side  of  (3.14)  is 
bounded  by  the  right-hand  side  of  (3.12)  since  0  <  x  <  K\fi . 

Case  3.  Ky/t  +  l  <  x  <  t  +  1  -  Ky/t  +  1.  Let  <5  >  0  be  small  and  <7  >  0  be  large.  By 
(3.13), 


L 


x/C—1 

•x+(t-x)(l—6) 


2 

(t  +  1  )“1_t(5  4.  i)-f+5e“(M+*)(*+i)  ds 


la  P  i  1  ,  1  _  (x-(a  +  l))2 

(t-x)~2  2x  2  +  2  (t  +  1)  2e  ds 


/  =  0(1) 

/•*+ l«- 

+  0(1)  / 

A/o-i 

+  0(1)  /  («-s)_5((+l-S)-;l(s+l)_^  +  5(t  +  l) 

Jx+(i-x)(l--<5) 

=  O(l)07(x,t;O,/i  +  e)  +  O(l)(l-3T)-i-?a:-‘i  +  ^  +  0(l)fty(*,t;  1,M +  0 


1  __  (t-j-1  -x)2 

“2g  (*i+e)(t+l) 


+ 


+ 


r  O(l)6>7(x,i;0,^  +  e)log(£  +  1),  if  P  =  3 
[  0,  otherwise 

(  O(l)01(x,  +  e)  log(t  +  1),  if  ol  =  1 

1  0,  otherwise. 
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which  is  bounded  by  the  right-hand  side  of  (3.12). 

Case  4 .  \x  -  (t  4-  1)|  <  K\Jt  +  1.  This  case  is  similar  to  Case  2. 

Case  5.  x  >  t  +  1  +  Ky/t  +  1.  Notice  that 

.  (a-(t  +  l))2  t 

0(l)e  *(*+u  c  if  0  <  s  <  i/2 
e  ,  if  t/2  <  s  <  t. 

This  case  is  similar  to  Case  1. 


□ 


Set 

ipa{x-,  t;  A)  —  [(a;  —  A(t  +  l))2  + 1  +  1] 

(3.15) 

</>a(x,tj  A)  =  (|x  —  X(t  +  1)1  +  t  +  l]~a/2. 

Lemma  3.4.  Te/  the  constants  1  ^  a  <  3,  fi,  fi  >  0  gtio?  X  ^  Xf .  Let  k  ^  0  be  an  integer. 
If  a  function  h(x,t)  satisfies 

h(x,t)  =  0(l)9a(x,t-i\',tJL'), 


dk 


(3.16) 


dxk^Xy^  ~  ^  >/0» 

/if  +  A  hx  —  —hxx  =  O(l)0a+2(s,  t;  X* ,[i')  +  (0(l)#a+i (x,  t;  A7 ,  /x/)); 


<9* 


^  ^a?*)  —  0(l)6a+k+2(x,  t\  X' ,  //)  +  (O(l)0a+k+i(x,t]  A',^'))^, 

then  for  any  fixed  e  >  0,  K  >\X-  A'|,  and  all  -oo  <  x  <  oo,  t  >  0, 

(3.17)  f*  r  (t  -  dk+1 

Jo  J — oo 


dyk+l 


h(y,  s)  dyds 


=  0(1) (<  +  1)  *  (x,  t-  A)  +  <?min(a,2) (x,  t;  A',  n*  +  e) 

+  |x  -  A(t  +  l)|~?|x  -  A'(f  +  1)|_5 

•  char{min(A,  A')(t  +  1)  +  KVt+1  <x<  max(A,  A  ')(t  +  1)  -  KsftTl}) 
where  /i*  —  max^,/). 


Proof  Again  it  is  sufficient  to  prove  the  lemma  for  A  =  0,  A'  =  1.  We  may  assume  (jl  <  g! . 
Otherwise  we  can  always  replace  //  in  (3.16)  by  /i.  Denote  the  left-hand  side  of  (3.17)  as 


(3.18) 


I  —  h  +  I 2  > 
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°  _i  (x-v)2  dk+1 

( t-s )  2e  h{y,s)dyds, 

1  (x-y)2 

(t-s)~*e  dyk+1h(y,  s)  dyds. 


We  assume  t  >  4.  The  case  that  i  <  4  is  a  consequence  of  Lemma  3.3.  By  (3.16)  and  (3.9), 


(3.20)  h  =  0(1) 


h(y,  s)  dyds 


fVl  k  1  1  ct  1  («-(a+l))a 

0(1)/  (t-s)_2~2(t  +  1)  2  (s  +  1)  2  +2e  +o(*+d  (is, 

Jo 


where  e  >  0  is  a  small  constant. 


Case  1.1.  x  <  0.  Then 


h  =  0(1)  J  (t-s)”5-i(t+l)-5(s  +  l)  ?  +  2e  «-'+«)<■+')  ds 

(3.21)  =  0(l)(t  +  l)_5_?_4e_('''+,)<t+1) 

=  0(l)(f  +  l)  *9a±i  ( x .  t;  0.  //  +  e). 

Case  1.2.  0  <  x  <  Cs/t+l  with  large  O.  By  (3.20), 

h  =0(1)  (t-s)-^{t+lT^{s+\)~^ds 

(3.22)  Jo 

=  0(l)(t  +  l)-^-»-i  =  0(l)(t  +  l)~$0«±i(x,t\0,ij!  +e). 


Cose  1.3.  x  >  Cy/t  +  1.  By  (3.20), 
r  Vi 

(3.23)  h  =  0(1)  /  (t-s)_  =  _5(t  +  i)-5(s  +  i) 

Jo 

=  0(l)(t  +  1)"*0«±I  +  26). 


a  ,  1  _*2(1-2/cU 
“T+2e  (/i'  +  e)(*+l)  rfs 


Notice  that  from  (3.8)  and  (3.15), 


(3.24) 


da(x,t;\,fi)  =  0(l)V'a(a:,f;  A) 
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for  any  a  >  0,  fi  >  0  and  A.  Thus  (3.21)-(3.23)  imply  that  I\  is  bounded  by  the  right-hand 
side  of  (3.17). 

Denote  the  solution  of 


(3.25) 


\i  ,  ^dk+1h 
Wt  =  4Wxx  +  Q ^5+r 
w(x,  t0)  =  0 


as  w(x,  t;  t0).  Duhamel’s  principle  and  (3.19)  imply  that  I2  =  w(x,  t;  y/i).  The  initial  value 
problem  (3.25)  can  be  rewritten  as 

,  ,  . — dkh.  n,  , dkh v 

(w  +  +  y/i^h 

. dkh w  „  _ dkh 


where 

(3.26) 
Hence 

(3.27) 


Qk 

h{x,t)  =  +  hx  —  ^ hxx)(x ,  t). 


k)L  dk  r 


dk  f°°  !  (x_, 

h  =  ~  Vi*K-Q^h(x,t)  + J  (t  - 

ft  f°°  1  (x_y)2 

+  /  /  (£  -  s)~ie~~rf*^h(y,s)dyds. 

Jy/t  J -oo 

Denote  the  first  integral  in  (3.27)  as  J21.  From  (3.16)  and  (3.9)  we  have 


21 


(3.28) 


/OO  { x_y )2 

(t  -  \/t)“i^iO(l)e_«‘+‘)(”-v^)h(j/,  v/t)  dy 

■oo 

=  o(i)(t+ir^-?  f 

=  o(i)(t  +  i)-t-^e-^a2 

—  +  1)“2<9o+i  +  2s), 


__  (x~v)2  (v-(\/t+l))2 

g  (M+e)(t->/t)  M/(v/i+1)  fly 


where  in  the  last  step  we  have  considered  x  <  0,  0  <  x  <  Cy/t  +  1  and  x  >  C\Jt  +  1 
respectively  for  a  large  C. 
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Denote  the  second  integral  in  (3.27)  as  122-  By  (3.26)  and  (3.16), 


n°°  1  k  (x-y)2 

(t~-s)~ 2  2  0(l)e  (M+e/2)(t-,T[6la+2(^,s;l,/x  ) 

-oo 


1  fi  roc  J  (x-y)2 

+  (<-s)  ^a+1(y,s;l,^')]^s+  /  /  (*-s)  2e 

^t/2  V— oo 

•  [0(l)6»a+A:+2(y,s;l,M')  +  (0(l)ea+fc+i(y,s;l,/i'))y]^s 

y>£  ^»00  1  (x-y)2 

=  0(l)(t  +  l)_i  /  /  (t  -  s)_i e“ o*+./»)(t-.)  0a+2(y,  s;  1, //)  rfyds 

J  \/i  J  —  oo 

k  ^*00  (x_y)2 

+  0(l)(t  +  l)-*  /  /  (t  -  s)_1e_ (M+«/a)(*-») eo+i(y, s;  1,  //) dyds 

*/  v7  •/ — oo 

=  /(1)  4-  I(2) 

—  J22  +  J22  * 


(3.29) 


We  now  estimate  the  first  term  on  the  right-hand  side.  By  (3.8)  and  (3.9), 

’Vi 

Case  2.1.  x  <  0.  Then 


r(D 

22 


k  1  a  1  (x-(s+l))2 

=  0(l)(t  +  l)“2-2  /  (s-hl)-2"2e  (^+E/2)(t+l )  ds. 

JVt 


(3.30) 


,-v  fc  a+i  C  i  _  (£±il 

1$  =0(l)(t  +  l)-5__r-  /  (i  +  l)-ie 

J  Vi 

=  0(1)  (t  +  l)“^a+i  (x,t;0 ,//  +e). 


+e)(t  +  i)  ds  e  (^/+e)(*+D 


Case  £.2.  0  <  x  <  C\/t  T  1  with  large  C.  Then 
4«  =  o(i)(t+i)-t-^  /  (t+ 

(3.31)  •'v't 


1  (x-(s+l))z 

2e  (m/+^)(*+i)  ds 


=  0(1)(£  +  1)  2  “41  -  0(l)(t  +  1)  20o+i  (x,t;0,/i'  -he). 


Case  £.3.  x  >  C\/t  +  1.  Then 


r  r2x/C—\  i  x2 

/W  =  0(1)(*  +  (s  +  l)“S-3e~«*'+‘><‘+i>  ds 

pt 

+ 

(3.32)  J2x/c- 

.  0  +  1  X - ■■■■-+- -  Q  +  l 

(t  +  1)  4  ,  ==e  (^'+o(t+D  +  x  2 


r 

J2x/C-1 

-0(l)(t  +  l)-» 


a  1  „  («-(«+!))■_ 

(C“2"2e  <M'+e)(t  +  l)  ds 


\/£  T  1 


=  0(l)(t  +  1)  0). 
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where  the  second  integral  is  zero  if  2 x/C  -  1  >t,  and  we  have  used  (3.24). 

Equations  (3.30)-(3.32)  imply  that  1$  is  bounded  by  the  right-hand  side  of  (3.17).  As 
for  the  second  term  1$  on  the  right-hand  side  of  (3.29),  we  apply  Lemma  3.3  to  obtain 
the  desired  result.  Together  with  (3.27),  (3.16)  and  (3.28),  I2  is  bounded  by  the  right-hand 
side  of  (3.17).  □ 


Next  we  consider  the  evolution  of  the  heat  kernel  with  generalizd  diffusion  waves. 


Lemma  3.5.  Let  a  >  0,  /?  >  0,  //  >  0  and  A  be  constants.  Then  for  all  — oo  <  x  <  oo, 
t  >  0,  we  have 


(3.33) 


(t-s)-a(t  +  1 


a  (x-y-X(t-s))2 

2  g 


(s  +  1)  *t/>3/2(y,s;  A)  dyds 


f  0(l)(t  +  l)-i^3/2(x,t;A)log(t  +  2),  if  a  =  1  or  0  =  3/2 
1  0(l)(t  +  l)~iip3/2(x,t;\),  otherwise, 


where  y  =  min  (a,  1)  +  min(/3, 3/2)  -  1. 


Proof.  Without  loss  of  generality  we  may  assume  A  =  0.  Denote  the  left-hand  side  of 
(3.33)  as  I.  Then 

I  =  h+h, 


where 


(3.34) 


h 


h 


°(1)  L 

0{1)  L  /v|<Vi+T 


(t  -  s)""1^  +  1  -  s)~%e~ 
(t  -  s)~l(t  +  1  -  s)~^e~ 
(t  -  s')~^ (t  +  1  -  s)~%e~ 


(x-y)2  _  0_  3 

M(t-S)  (s  +  1)  2  4  dyds 
^-■)^+|(y,  s;  0,/z)  dyds , 
n(t-S)  ( s  I)-  2  \y\~i  dyds. 


Applying  Lemma  3.2  to  /1?  we  see  that  it  is  bounded  by  the  right-hand  side  of  (3.33).  To 
estimate  I2,  we  consider  two  cases. 
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Case  1.  \x\  <  2 yjt  +  1.  Then 
•*/2 

/2  =  0(1) 


(3.35) 


/■t/z  „  /*oo 

/  (t-s)_1(t  +  l-s)~?(s  +  l)“  /  dj/tte 

JO  J  \/s+l 

+  0(1)/  (t  —  s)_1(i  +  1  —  s)_^(s  +  l)_2_l  /  e- dt/ds 

Jt/ 2  J  |!/|>V*+1 

/•*/2  a  i 

0(1)(1  +  /  («  +  l)  2  «ds 

Jo 

+  0(l)(t  +  1)”^”*  [  (t-  s)~^(t  +  l  -  $)~i  ds 

J  i/2 


_  f  0(l)(t  +  1)"2-|  log(t  +  2),  if  a  =  1  or  0  =  3/2 

|  0(1)(£  +  1)_ 2-| ,  otherwise, 

which  is  bounded  by  the  right-hand  side  of  (3.33). 


Case  2.  |x|  >  2 \Jt  +  1.  We  consider  x  >  2 \Jt  +  1.  The  case  x  <  —2 y/t  +  1  is  the 
same.  Considering  the  integrals  with  respect  to  y  on  (— oo,  —y/s  +  1  ] ,  [Vs  -FT,  x/2]  and 
[x/2, +oo)  respectively  in  (3.34),  we  obtain 

I2=0(  1)  [  (t  -  s)-1(£  +  1  -  s)_^e”4^*— >  {s  +  1)“ 2-i 

Jo 

+  0(1)  f  (t  -  s)"*(t  +  1  -  s)_^(s  +  1)"^"'  ds 

Jo 

(3.36)  =0{l)(t  +  l)-1-^-^  [  '  (s  +  l  +  -i 

Jo 

rt  2 

+  0(1)(<  +  1)“  2 -4  /  (t  -  s)_1(t  +  1  -  s)-f  e“4^Xt-s>  ds 

J  t/2 

+  0(l)(f+l)"*aT*. 

The  first  term  is  bounded  by  the  right-hand  side  of  (3.33),  cf.  the  last  step  in  Case  1.  The 
third  term  is  0(1)(*  +  l)”7/2^3/20M;O).  To  estimate  the  second  term,  set  rj  =  x/Vt-s 
and  it  becomes 


0(l)(t  +  1) 


*/ 


2r]a' 


»  (i r ]2  +  x2 

VtJ 2 


\a/2 


U—  , 

dr] 


—  0(l)(t  +  1) 

—  0(l)(t  + 1)- 


0_l_OL  _xl_ 

2  4  2  £  4 fit 


r]a  1e  dy 

=  0(1)(£  +  1)  2^3/2(^7^0)- 


□ 
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Lemma  3.6.  Let  the  constants  a  >  0,  (3  >  0,  (i  >  0,  and  A  ^  A'.  Then  for  any  fixed 
K  >  2|A  —  A'|  and  all  —  oo  <x<oc,t>0,  we  have 

ft  I*00  a  (a;  — ;/  —  A(t  —  s))2  a 

/  /  (*“«) (t  +  l-  s)~^e  (s  +  l)_2^3/2(j/,s;A,)dj/d5 

o  0  J —oo 

=  0(l)(i  + 1)  2  [^3/2(2,  f;  A)  +  ip3/2(x, t;  A')] 

+  0(l)|as  -  A (t  +  1)|“5  mi"W’§)-3  j®  _  A'(l  +  l)|“i  min(“'1M 

(3.37) 

•  char{rain(A,  A  '){t  +  1)  +  Ky/t  +  l  <x<  max(A,  A  '){t  +  1)  -  KVt+ 1} 

(  0(1)(i  +  l)“=  log(t  +  l)[^3/2(a:,t;A)  +  ^3/2(®,t;A,)])  if  a  =  1 
+  \  0(l)(i  +  I)-?  log(t  +  l)^'3/2(^s  Z;  A),  if  a  ^  1  and  f3  =  3/2 
[  0,  otherwise , 

where  7  =  min(a,  1)  +  min(/?,  3/2)  —  1. 

Proof.  Without  loss  of  generality  we  assume  A  =  0,  A'  =  1.  Denote  the  left-hand  side  of 
(3.37)  as  I.  Then 

/  =  /j  +  J2j 


m  I'J 


|j/-(s  +  l)|<\/5+I 


{t  -  s)  1  (t  T  1  —  s)  “e  mX(*-S)  (s  +  1)  2  4  dyds 


pt  p 00 

0(1)/  /  (*_S)-I(t  +  1-S) 

J  0  J —00 

0(1)  f  f 

Jo  J\y- (s+1)|>a/s+T 


_  (sc  —  I/)2 

e  ep+3/2{y,s;l,fj,)dyds, 


(t  —  s)  (t-hl-s)  2e  /*(*—>  (s  +  1)' 
■  \y  ~  (s  +  l)|“^  dyds. 


Applying  Lemma  3.3  to  Ji,  we  see  that  it  is  bounded  by  the  right-hand  side  of  (3.37).  To 
estimate  I2 ,  we  consider  five  cases. 

Case  1.  x  <  —2y/t  +  1.  In  this  case 


h  =0(1) 


pt  px/2 

/  /  (t 

J  0  J  —  00 


s)  ^■(t  +  l-s)  “e  **<*-■)  (s  +  1)  2  (s  H-  1  —  — )  2  dyds 


+omJ  I. 


t  rs+i~^+i 


0  Jx/2 


(t  -  s)  x(t  +  1  -  s)  2  e  4m(*-S)  (s  +  1)  2  (5  -(- 1  _ 
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t  /*oo 


+  0(1) 


a 


(t-  s)-1(t  +  l-  s)~%e  ^(5  4-1)  *(y-(s  +  1))  2  dyds 


+1++^+! 


0(1)  [  {t-s)-i{t  +  l-s)-%(s  +  l)-%\x\-§ds 

Jo 

+  0(1)  [  (t-  s)_1(t  +  l  -  5)-fe_*^^(s  +  l)_^"i  ds. 
Jo 


The  right-hand  side  is  exactly  what  we  had  in  (3.36),  which  is  bounded  by  the  right-hand 
side  of  (3.33),  hence  of  (3.37). 


Case2.  |x|  <  Ky/t  +  1.  From  (3.38), 

rt/2  Q 

h  =  0(1)  /  (i  -  s)-1(i+  1  -  s)-Y(s  +  l)-2-3 
Jo 

ft 


ds 


+  0(1)/  (t-s)  2(t  +  l-s)  a(s  +  l)  2  4  ds, 

Jt/2 

which  is  what  we  had  in  (3.35)  and  hence  bounded  by  the  right-hand  side  of  (3.37). 
CaseS .  Ky/t  +  1  <  x  <  t  +  1  -  Ky/t  +  1.  From  (3.38), 


(3.39) 

Here 


=  0(1)  [  [  +0(1)  [  f  =hi+ fa- 

Jo  J y<s+l  —  a/s+1  Jo  J + 

;—y/x 


121 


(3.40) 


fX  V*  __  a  _  (*-(s  +  l))2  __  0  _  1 

=  0(1)  /  (t  -  s)  x(t  +  1  -  s)  2e  (5  +  1)  2  4  ds 

Jo 

+  0(1)  fX  (t-s)+t+l-s)=(s  +  l)=“*<fe 

JX—y/x 

+  0(1)  f  /  +1 

J X  +  y/x  Jy<—. 5 


a  _ 


(t  -  s)  1  (t  +  1  -  s)  2  e  '*<*-•>  (s  +  1) 


/•t  /*s  +  l  —  \/  s  +  l 

(s  +  1  —  x)~2  dyds  +  0(1)  /  /  (t  —  5)  J(t  +  1  —  s) 


-  .<d~2 


■  e  ( ^0-4  (s  +  1)  2(s  +  l-y)  2  dyds 

=  /211  +  7 212  +  7^213  +  7214- 
Estimate  each  term  respectively.  We  have 

1211=0(1)  J  (t-  |r1-fe“*+s  +  l)=-i 


[X  'f*  1  a  1  (s  +  l-s)2  1 

+  0(1)  /  (t  -  x)_2_T(£  -  s)  2e  ^-s>  x  2  4  ds 

Jx/2-l 


=  0(l)(t  +  l)  1  *e  4+* 


y»t/2  ^  3 

/  (s  +  l)_  =  _Jds  +  0(l)(t-x)_i-*a:'“2-i 

Jo 
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which  is  bounded  by  the  right-hand  side  of  (3.37). 


PX  +  y/x 

hl2  =  0(1)  /  (t  -  X  -  y/x)~*~%  (x  \fx  +  1)_  2 -! 

J  X—a/x 


la  0  1 


=  0(l)(i  +  1  —  x)  2  2x  2  4  , 

which  is  also  bounded  by  the  right-hand  side  of  (3.37).  For  e  >  0  small, 


rt 

^213  =0(1)  /  (*  -  s)~^(t  +  1  -  $)~%x~%(s  +  1  -  x)~%  ds 

J  x-\-y/x 

=  0(1)  f  (t  -  +  1  -  x)~%  ds 

J  X  +  \/x 


+  0(1)  f  (t  —  s)  2(t  +  l  —  s)  2  x  2  (t  —  x)  2  ds 

J  t—e(t—x) 


+  0{l)x  *(t-x)  2  /  +  /  (t  -  s)-z(t  +  1  -s)-;5  ds. 

\Jt~e(t-x)  Jt-eJ 


The  last  term  is 


o(i).-f (. -»)-!+( 0(111  ;(1  - iia*1 

[  0(l)x  2  (t  -  x)  2  log (£-x),  if  a  =  1. 

Notice  ( t  —  x)-1/2  =  0(l)x-1/4.  Then  the  above  expression  can  be  rewritten  as 

+  f  0(l)r?  Iogt[ip3/2(x,t;0)  +  ip3/2(x,t;l)],  if  a  =  1 
1  0,  if  a  ^  1. 

Thus  the  right-hand  side  of  (3.41)  is  bounded  by  the  right-hand  side  of  (3.37). 

J2i4=0(l)  /  (t  -  s)-1(t  +  1  -  s)~^e  w-*)  (s  +  l)_2“4ds 

J  x  +  i/x 


„  ,  x  _0_1.  .  I  _  a  f*  £^t  X\  1  _  (s  +  l-x)2 

=  0(1)X  2  4(t  —  X)  2  2  /  (t  -  S)  2g  4M(i-S) 

./  s-l-yT 

+  0(l)x“^_i  f  (s  +  1  —  x)~2(t  +1  —  s) 
Jt~e(t~x) 

r°°  2 

=  0(l)x_2  “i (t  —  x)~i~ 2  /  dr) 

Jo 

-\-0(l)x~2~i(t-x)-2  f  (t  +  l-s)_2c 

Jt—e(t—x) 


- s)-*ds 


0(l)x" 2  i  (t  -  x)“i"*  min(a’1}. 
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Thus  we  have  estimate  I21 •  In  (3.39),  I22  can  be  estimated  in  a  similar  way: 

\—y/x  -S+i±* 


122  — 


1  a  (g-(s  +  l)U  /  x  0 

(t-s)  1(t+l-s)  2  e  V(t-s)  (s  + 1)  2 


5+l  +  yS+I 


pX~y/X  p 

0{1)l  I 

■(y-(s  +  l))-idyds  +  0(l)  [  [  (t  -  s)~\t  +  1  -  s) 

Jo  Jy>*±¥JE 


■e  (s  4-  1)  2  (x  —  (s  T  1))  2  dyds  + 12\2 

+  0(1)  f  {t  -  s)-1(t  +  1  -  s)_te"(Mt-^ (s  +  l)-^-*  ds 

J  X  +  y/x 

=  I22I  +  I222  +  h\2  +  -^214? 


where 


I221  — 


/*X  1  a  (x-(3  +  l))2  x  _  0  _  1 

0(1)/  (t-s)_1(t  +  l-s)~2e  (5  +  1)  2  4  ds, 

Jo 


which  is  /211  in  (3.40); 

»X~y/x 


J-222 


=  0(1)/  (t-s)+t  +  l-s)“?(s  +  l)  *(x-(s  +  l))  *  ds 

Jo 


=  0(l)(t  +  1  -  x)-i-t  '  {s  +  1)- 
=  0(l)(i  +  1  -  x)-i-?x~i  min(/3’4)- s 


£  _3  , 

2  a:  2  ds  + 


Jx/2 


2  (z  -  (s  +  1))  2  ds 


which  is  bounded  by  the  right-hand  side  of  (3.37);  I212  is  the  same  as  in  (3.40),  and  J214 
is  in  (3.42). 


Case  4.  \x  -  {t  +  1)|  <  Ky/t  +  1.  If  (5  ^  3/2,  this  case  is  the  same  as  Case  2.  If  (3  =  3/2, 
we  divide  I2  into  two  parts.  One  is  the  integral  with  respect  to  s  on  [0,  t/2],  and  the  other 
is  on  [t/2,t].  Clearly  the  second  part  is  bounded  by  the  right-hand  side  of  (3.37).  We 
consider  the  first  part, 


rt/2  /  /»  rh(s+ l+x)\ 

Jo  \  J y<s+l  —  y/s+l  J $■+■!  + y/s  +  l  J 


0(1) 


(t  -  s)  1{t  +  1  -  s)  2  e 


(x-(3  +  l))- 
4M{t-s) 


(s  +  l)~i\y  -  {s  +  l)\  2  dyds 


+ 


ou)  rj , 

Jq  Jy>k 


(x-y) 


(t  -  +  1  -  s)~^e  “<*-»>  (s  +  1)  «(x-(s  +  l))  2  dyds 


y>h(s  4-l+x) 
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ft/2 

0(1)  /  (t  —  s)  1(t  +  1  —  s)~  2e~64^:  (5  +  i)-i  (fs 
Jo 

ft/2 

+  0(1)/  ( t-s )  2(i  +  l-s)-9(s  +  l)-?(z-(s  +  l))-t 

J  0 


=o(i)(t  +  i)-f-J, 


which  is  bounded  by  the  right-hand  side  of  (3.37), 


Case  5.  x  >  t  +  1  +  Ky/t  +  1.  In  this 


h  =0(1)  f(f  +  /-|(S+1+*)N\  (f  _  a)-i(i  +  !  _  s)-f 

J°  \Jy<S+l-y/s+l  Js+  l  +  v/I+I  J 


(x-(s  +  l)): 
g  4/x(t  —  s) 


(s  +  1)  2  |l/  —  (5  +  1)1  ldyds  +  0(X)(  f  1 

JO  + 


•(*  +  !-«)  ?e  (s  +  l)_^(x-(s  +  l))-idyds 

=  0(1)  f  (t  -  «)"*(*  + 1  -  s)~Te_<Vu-1.))  (s  +  l)~f-4  ds 

V  0 

+  £*(1)  f  (*  _  5)  2  (i  +  1  -  s)-"2  (s  +  1)“  2  (x  _  (5  _|_  l))-f  ds 

J  0 

=  0(l)(t  +  [t/2{s  +  ds 

Jo 

+  0(l)(t  +  1)_  2  -4  f  (t  -  S)-1-fe-< 

Jt/2 

+  0(l)(f  +  1)"*  (*  -  (t  +  l))-i  f  l0g(i  +  2)5  if  a  =  1 

l  1,  otherwise. 


The  first  two  terms  are  bounded  by 


0(l)(t  +  1)-*  03/2(x,  t;  1,  8m)  =  0(l)(i  +  l)-^3/2(x,  t;  1). 


Thus  /2  is  bounded  by  the  right-hand  side  of  (3.37). 


Lemma  3.7.  Let  a  >  0,  0  <  (3  <  2,  (j,  >  0,  and  X  be  constants.  Then  for  all  ~oo  <  x  <  oo 
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t  >0,  we  have 


(3.43) 


n°°  (f  -  s)-1^  +  1  -  s)“7e-<*  ee-*)’”  (s+1)  ?<t>i(y,s;  A)^3/4 (y,  s;  A)  dyds 

-OO 


0(l)(t  +  l)-^  log(t  +  2)<j>i(x,t\  \)^s/i{x,t;  A),  if  a  =  1  or  1  <0<  3/2 
0(1)  (t  +  1)-^  <l>i  (x,  1;  A)V>3/4(*i  <;  A),  otherwise 

(  0(l)(t  +  l)_^log(t  +  2)^3/2(x,«;A),  i/ a  =  1  or  (3  =  3/2 
+  \  0(l)(t  +  l)-^V’3/2(a:,<;  A),  otherwise, 


where  71  =  min(a,  1)  4-  (min(/3, 1)  +  min(0, 3/2))/2  -  1,  72  —  min(a,  1)  +  min(/3, 3/2)  1. 

Proo/.  Without  loss  of  generality  we  assume  A  =  0.  Denote  the  left-hand  side  of  (3.43)  as 
I.  Then 


rt 

=  /  (t-s)-1(t+l-s)“?(s  +  l)'! 

Jo 

e-§7^T  01(y?  S;  0)^3/4(y,  5;  0)  dy 

=  A  +  I2  +  ^3- 


7  +/ 

J|yj<\/s+l  « v/s+l<|yl<s+1 


/  ) 


ds 


By  (3.15), 

r  (s  +  i)-*,  li/l  <  V*  + 1 

</>i(j/,s;0)V'3/4(y,s;0)  =  0(1) <  |»|-*(«  +  i)-*,  vm<|y|<,  +  i 

l  |j/|-+  |y|  >  s  + 1. 

The  integral  I\,  corresponding  to  |y|  <  y/s  +  1,  is  estimated  in  the  same  way  as  in 

Lemma  3.5.  For  |x|  <  Vt  +  1,  h  +  h  is  estimated  in  a  similar  way  as  in  (3.35).  For 

y/t  +  T  <  x  <  t  +  1,  considering  y  <  x/2  and  y  >  x/2,  we  have 

I2  +  I3  =  0{  1)  [  (t-s)-1(t  +  l-s)_fe”I?^T(s  +  1)_'”’cis 
Jo 

+  0(1)  [  (t-s)-2(t  +  l-s)-»(s  +  l)-UCTx-^+<Tds 

Jo 

where  0  <  a  <  1/2.  The  first  term  has  been  estimated  in  (3.36).  In  the  second  term,  we 
let  <7  =  min(l/2, 1  -  0/2)  for  0  <  s  <  tj 2,  and  a  =  1/2  for  t/2  <  s  <  t.  Thus  the  second 


QUASILINEAR  HYPERBOLIC-PARABOLIC  SYSTEMS 


33 


term  is  bounded  by 


(3.44) 

0(l)(t  +  l)-i  *-i  +  | 


0(l)(t  +  l)-i  min<“'1)-f  log(t  +  1)3.-^ 
0{l)(t  +  log(t  4-  , 


0, 


Notice  that 


if  a  =  1 

if  1  <  fi  <  2 

otherwise. 


(t  +  1)  *x  4  =  O(l)0i {x,t\0)ip3/4(x:t-,0), 

and 

a;2-f  =  0(l)(t  +  l)i“S  for  1  <  fi  <  2. 


We  conclude  that  (3.44)  is  bounded  by  the  right-hand  side  of  (3.43).  For  x  >  £  +  1,  /2  +  J3 
is  estimated  in  a  similar  way  as  in  (3.36).  Similarly  we  can  prove  the  case  x  <  —  y/t  +  1. 


□ 


Lemma  3.8.  Le£  the  constants  a  >  0,  0  <  fi  <  2,  (i  >  0,  and  A  ^  A'.  Then  for  any  fixed 

K  >  2 1 A  —  A'|  and  all  — oo  <  x  <  oo,  t  >  0,  we  have 
(3.45) 

(s  +  1)”  0! (*/,  5;  A')^3 /4(y,  s;  A')  dyds 

«0(l)(t  +  l)-^ 

'  <^i(^,  t;  A)^3/4(£c,^;  A)  -F  ^!(or,  A'),  i/a  #1,  fit  [1,  |] 

•  <  \og(t  +  2)</>!  (or,  t;  \)ipz/4(x,t;  A)  +  <j>i(x,t\  X'tyz/^x,  ^  A'),  i/a^l,  fie  [1,  |] 

w  M^  +  2)[01(a:,<;A)^3/4(^,i;A)  +  ^i(x,t;AO^3/4(^^;A/)],  »/a  =  l 


a 


a  (x-y-A(t-s))- 


(t  —  s)  1(i+l—  s) 


f  ^3/2 (®» t;  A)  +  ^3/2 (x,  £;  A'),  i/a^l,/?^~ 

+  0(l)(t  +  l)-:?  j  log(t  +  2)Vty2(M;A)+ife/a(*,t;A'),  /3  =  § 

i  l°g(i  +  2)[^3/2(x,t;A)  +  V3/2OM;  A')],  i/a=l 
+  0(l)|a;  -  A(t  +  1)|_5  min(/5,2)- J+e^  _  +  ]_) | - i  min(a,l)-i 


•  char{min(A,  A')(t  +  1)  +  KVt+1  <x<  max(A,  A')(f  +  1)  -  KVt  +  1}, 


where  71  =  min(a,  1)  +  (min (fi,  1)  +  min (/?,  3/2))/2  -  1,  72  =  min(a,  1)  +  min(/3, 3/2)  -  1, 
and  e  >  0  arbitrarily  small. 


Proof.  Making  use  of 

0i(2/,s;1)V’3/4(2/,s;1)  =  (s  +  l)_<T|y  -  (s  +  l)r?+<T,  0  <  a  < 

2 

the  proof  of  this  lemma  is  totally  parallel  to  the  proof  of  Lemma  3.6.  □ 
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Lemma  3.9.  Let  /i  >  0 ,  a  and  X  be  constants.  Then  for  all  -oo  <  x  <  oo,  t  >  0,  we 
have 

(3.46)  [  e~('t~sS)/l*'il>3/2{x  -  <r(t  -  s),s;\)ds  =  0(l)^3/2(aJ}  t\  A). 

Jo 

Proof  By  (3.15),  the  left-hand  side  of  (3.46)  is 

I  =  f  e-(t~s)/tl[{x  -  \(t  +  1)  +  (A  -  a){t  -  s)f  +  s  +  1]~3/4  ds. 

Jo 

If  |x-A(t  +  l)|  <  \/tTl,  then 

I  =  f  e~{t~s)/lJ0(l){s  +  l)-3/4  ds  =  0(l)(t  +  1)_3/4  =  0(1)^372(0:,  t;  A). 

Jo 

If  \x  -  A(i  +  1)1  >  y/t  +  1,  we  consider  two  cases: 

Case  (i).  (x  -  X (t  +  1))(A  -  <r)  >  0.  In  this  case 

I  =  f  e~{t~s)/fiO(l)\x  -  X(t  +  1)|_3/2  ds  =  0{l)\x  -  X (t  +  1)|~3/2  =  0(l)ip3/2{x,t;  A). 

Jo 


Case  (ii).  (x  -  A (t  +  1))(A  -  a)  <  0.  In  this  case 


I  = 


f  e~s/tiO{  1) 

Jo 


x  —  X(t  +  1) 


+  5 


+  1 


-3/4 


ds. 


X  —  a 

To  bound  the  integrand,  consider  whether  |(x  -  A (t  +  1))/(A  —  o)  +  5I  ^  2\(x  ~  ^  + 
l))/(A-a)|.  Thus 


-  0(1)  f 

Jo 


e  “ 


x  —  A  (t  +  1) 


A  —  a 

3 


1  ^ 

ds  4-  0(1) 

Jo 


x  — A(t  +  1) 

e  2/.(a— <*)  ds 


=  0(l)|x  -  A (t  +  1)|“3  =  0{l)ip3/2(x,t\  A). 


□ 


Similarly,  we  can  prove 


Lemma  3.10.  Let  \i  >  0,  a  and  X  be  constants.  Then  for  all  -00  <  x  <  00,  t  >  0, 


/ 


-  a{t  -s),s;  A)t/>3/4(x  -  a(t  -s),s;  A)  ds 

=  0( t;  \)ip3/4{x,  t;  A). 


(3.47) 


4.  Systems  with  Diagonalizable  Linearizations 


In  this  section  we  prove  Theorem  2.9  for  the  systems  whose  linearizations  are  diagonal¬ 
izable.  The  nondiagonalizable  case  will  be  discussed  at  the  end  of  Section  8. 

The  Cauchy  problem  we  are  now  considering  is 

Ut  +  f{u)x  =  Buxx ,  —  oo  <  x  <  oo,  t  >  0, 

(4.1) 

u(x,0)  =  Uo(x),  -oo  <  X  <  00. 

Here  B  is  a  constant  matrix,  and  there  exists  a  symmetric  positive  definite  matrix  Ao, 
such  that  Ao/'(0)  is  symmetric,  and  A0B  is  symmetric  positive  definite.  Let  A*,  i  —  1, 

s 

. . .  ,  5,  be  the  distinct  eigenvalues  of  /'( 0)  with  multiplicity  m*  =  n,  and  U  and  r* 

i— 1 

be  the  matrices  consisting  of  left  and  right  eigenvectors  kj  and  r^,  respectively,  of  /'( 0): 

(4  2)  /'(°)»«  =  Vy, 

T'i'j '  ~  fijj '  i 

1 ,  .  .  .  ,  5,  j  1 ,  .  .  .  ,  TTli ,  j  —  1,  .  .  .  ,  Tfli' . 

Then  we  further  assume  that  B  satisfies 


(4.3)  kBrj  =  0  for  i  ^  j,  z,  j  =  1,  . . .  ,  s. 

For  the  initial  data  we  assume  that 

«o(z)  =  0(1)(1  +  |a;|)-3/2,  ^  e  L°°( K), 

/X  rtOO 

uo(y)  dy  e  L1(R“),  v+(x)  =  u0(y)dy  €  i1(K+), 

-oo  L 

llWolU°°  +  ||v_||Li(_oo,0)  +  ||v+||li(0,oo)  +  sup  {(1  +  |t|)3/2|w0(^)|}  =  6*  <  1. 

— oo<x<oo 

Decompose  the  solution  w  to  (4.1)  in  the  right  eigenvector  directions, 

s 

(4.5)  u(x,t)  =  ’^2riui(x,t),  Ui{x,t)  =  ku{x,t),  i  =  l,...,s. 

i= 1 
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Let  Oi(x,t),  i  =  1,  •  •  •  ,  s,  be  the  self-similar  solution  to 


(4.6) 


On  4-  A iOix  +  (0)(r*^u ri@i)x  —  liBriOi; 


ip3/2{x,t;\i)  +  '^2i>{x,t-,Xj) 


/oo  r°° 

8i(x,  t)  dx  =  /  k u0(x)  dx  =  Si, 

-oo  J- oo 

given  by  (2.15),  (2.16)  or  (2.18).  Our  purpose  here  is  to  show  that  for  all  -oo  <  x  <  oo, 
t  >0,  2  =  1,  ...  ,5, 

Vi(x,t)  =  Ui(x,t)  -  Oi(x,t)  =  0(1)63 

(4.7) 

Vix(x,t)  =  0(1)6*  (t  +  l)-5^4, 

where 

t/?a(^,t;A)  =  [(x- A(f +  l))2  +  t  +  l]  a/2, 

^  ^  ^(x,t;  A)  =  [\x  -  A(t  -1-  1)|3  +  (t  +  l)2]  1/2- 

Under  our  assumptions,  equation  (2.10)  satisfied  by  Ui  takes  the  form 

u it  +  A iUlx  +  hif"(0)  (XrJuJ’  Y^rjUj)  =  hBnuixx  -  §i, 


j=  1  J  =  1 


where 


1  /  s  ^  5  ^ 

§i  =  <</(«)*  -  AiUte  -  2^/"(°)(^rJuJ’SrA 
=  (0(l)lw|3)a:  =  0(l)|u|2|u*|. 

Thus  with  (4.6)  Vi  satisfies 

(4.9) 
where 

s 

(4.10)  gi  =  -'<X)/"(o)Mj*rA 


7  =  1  J=1 


Vit  T  At'Uix  —  liBTiVixx  4“  Qix 
V ixt  4"  AjU^xx  —  I'iBTiVixxx  4“  $ix 


j=l 


+ 


0(1)  X  N(M  + N)  +  °(i)  X)  m2  +  X  !^l3’ 

m  J=1  >=1 


+ 


o(i)  d^-iivjb*i  +  imm  +  ieti2i^i + ww) + °(uXim*i- 

,./c= i  ^ 


QUASILINEAR  HYPERBOLIC-PARABOLIC  SYSTEMS 


37 


By  Lemma  2.1  we  may  assume 


(4.11) 


kBn  =  diag(/t, 


with  mj  >  0,  j  =  1,  . . .  ,  rrii,  since  for  each  1  <  i  <  s  we  can  perform  the  same  hnear 
transform  to  vt  and  9,  in  (4.9)  and  (4.10)  to  diagonalize  UBrt.  Let 


(4.12) 


Vi  =  vimj)‘. 


Duhamel’s  principle  and  (4.9)  imply  that  for  1  <j<  mi, 


(4.13) 


Vij(x,  t)  =  I  +  11  +  III  +  IV  +  V  +  VI, 


with 

(4.14) 

I 


-£ 


t  pOO 


(x-y-X it)* 

e  4"«'  vij(y,0)dy, 


“-U. 


{x-y-Xjjt-T))2  r 

g  4  Mij(t-T) 


li]'52f"{°)(rk6k,rkdk)(y, 
k^i 


>t  poo 


III  = 

IV 


r  f  , _ 

JO  J—  oo  y/4irfiij(t  -  r) 

-u:  1  - 


{x-y-Xj  (t-T))‘ 


Uj  f"(0)(rkdk,rkvk)(y, r] 


t  poo 

0  J- oo  y/4nfjLij(t  -  t) 


k=\  J  y 

e  ^TrJ”  \0(l)J2mM  +  M(y,r) 
k^k'  ->  y 


dydr , 
y 

dydr , 
dydr , 


/‘oo 


V 


V7 


-a 

=  r  r  , _ 

«/o  i-oo  y' 47T/Liij (t  -  r) 


_  (j-y-Aj(t-r)) 


2  r 


oo  y/4wfJ,ij(t  -  T) 

t  r°°  1  ( x-y-Xiit-T ))2  r 

u - 

oo 


IO(l)^|»*|2(y, 

k= 1 

o(i)Ei«*i*(tf. 


fc=l 


dydr , 


dydr, 
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with 

(4.16) 


/=r  a 

J-oc  J 4?r 


//  = 


///  = 


V^TT^jt 

-ur 


(x-y-Xjt)* 

e  4uyl  vijy(y,0)dy, 


1 


(x-y-Aj  (t-r))^ 


t/2  /*oo 


n 


y/4vHij(t-T) 

1 


k^i 


dydr, 


J  yy 


_  (x-y-Aj  (t-r))' 


—  OO  \/47r /Xij  (£  -  r) 

S 

+  0(l)5>fc|2  +  N3)|  [y, t)  dydr 


0(1)  E  IflfclKI  +  0(i)  E  I^H^' 

k^fik' 


fc.fc'  = 1 


fc= 1 

■t  /'OO 


-i  yy 


u 


(x-y-A,(t-T))2  r 
4Mi7(t-x) 


i/2  7-oo  ~T) 


0(1)  ^  (l^fcllvfe'yl  +  l^fcyll^fc' 


+  |0fc|2  A'vl  +  |v*:l|wfc'wl)  +  0(1)  E 


k,kf= 1 

(2/,r)dydr. 


We  now  derive  a  priori  estimate  for  v*.  Set 

(4.17)  M(t)  =  sup  maxi  t>t(*,  t)[^3/2(',  t; A*)  +  ^#,t; A*)]  1 
o<r<ti<i<^l 

+  ||^ix(-,^)IUoc(T  +  1)5/4|* 


Then  for  — oo  <  x  <  oo,  t  >  0,  i  =  1,  .  •  •  ,  s 
(4.18) 


\vi(x,t)\  <  M(t)[il>3/2(x,t;\i)  +  ^T4>{x,t-,\k)}, 

k^i 

ViX{x,t)\  <  +  l)-5/4. 


We  estimate  each  term  in  (4.13). 

From  (4.5),  (4.4),  (3.1)  and  (4.6)  we  see  that 

^(£,0)  =  kuo(x)  -  6i(x, 0)  =  0(1)6*  (1  +  \x\)  3^2, 

Vix(x,  0)  =  0(1)6*, 

/oo 

|i>i(x,0)|dx  =  0(1)5*,  /  Vi(x,0)dx  =  0, 

-oo  J-°° 

roo  fx 

/  |j3i(x)|dx  =  0(1)5*  with  rii{x)=  Vi{y,0)dy. 
J  —  OO  0-00 


To  estimate  /  we  consider  two  cases. 
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Case  1.  \x  —  \it\  <  (t  +  l)2/3.  Let  the  jth  component  of  r]i  be  rjij.  By  (4.14), 

/°°  1  ( x-y-Xjt )2  rOO 

.»  4"y‘  ^{y) dy  =  °(1)t_1  L  !^(J/)I 


(4.20) 


sfinHjt 

=  0(1)6*  (t  +  1)-1  =0(l)8,i>(x,t-,\i). 


Case  2.  \x— A;l|  >  (t+1)2/3.  We  consider  a:— A*  /  >  (t+1)2/3.  The  case  a;— A;f  <  —  (t+1)2/3 
is  similar.  By  (4.14)  and  (4.19), 


/ 


(x  \it)/2  (x  — A,-t)2 

1  =  I  0(l)t~ie  ,6"‘T  |ny(y,0)|dj/ 


(4.21) 


+ 


/*°°  j  (tt-y-Ajt)5 

/  0(l)l_5e  4e«*  tf*(l+j,)-S  dy 


,()/2 


_  ,  .  1  t1/3  (*  —  Ajt)2  O  ^ 

-  0(l)t“2e-^  ^—,5*  +  O(l)0*(z  -  A*t)"2  -  O(l)0>(a;,*;  A,). 


Equations  (4.20)  and  (4.21)  imply  that 
(4-22)  I  =  Q(  l)0*#r,t;A*). 


To  estimate  //  we  set 


(4-23)  hk(x,t)  =  lijf''{0)(rk6k,rk6k)(x,t),  1  <  k  <  s. 

we  recall  (3.8), 

(4.24)  0a OM;  A,  A4)  =  (t  +  l)”f  e-"  m(U+d)}  . 

By  (3.1)  it  is  clear  that 

(4-25)  hk(x,t)  =  0(l)6*2e2(x,t;\k,2ixk), 

and 

hkt  T  A khkx  ftijhjc  xx 

=  2iijf"(0)(rkek,rk(ekt  +  \kekx))  +  (0(  I)\ek\\ekx\)x 

(4.26) 

=  o(\)(\ek\\ekt  +  xkekx\)  +  (0(1)1^11^1), 

=  0(l)6*%(x,  t;  Xk,  2nk  +  e)  +  (0(l)6*263(x,  f,  Xk,2nk  +  e))„ 
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where  (xk  >  0  is  the  maximum  eigenvalue  of  lkBrk,  and  e  >  0  is  an  arbitrarily  fixed  small 
constant.  We  now  apply  Lemma  3.4  to  II  in  (4.14)  with  a  =  2,  k  =  0,  A  =  A*,  and  X'  =  Xk: 


II  = - - - V  f  [°°  (: t 

JO  J-J 


t)  2  e 


t))2 

hf.y(y,T)  dydr 


(4.27) 


=  y^O(l)5*2[^3/2(x,i;  Aj)  +e2{x,t;\k,fj.*) 

+  \x-\i(t  +  l)\-1\x-\k(U-l)\-i 

•  char{min(Aj,  \k)(t  +  1)  +  K\ft  +  1  <  x  <  max(Aj,  \k){t  +  1)  —  Ky/t  +  l}] 
=  0(1)A'2  [rh/2(x,  t;  A,)  +  53  x{>(x,  t;  A*)] , 

k^i 


where  jx*,  K  >  0  are  constants. 

To  estimate  III  similarly  we  set 


hk(x,t)  =  lijf"(0){rk9k,rkvk)(x,t),  1  <  k  <  s. 


We  have  by  (3.1),  (4.18)  and  (4.8), 

hk(x,t)  =  0(l)\9k\\vk\  =  0(1)8*  M(t)95/2(x,t,\k,4fik), 


and  also  by  (4.9), 

hkt  4"  Xkhkx  —  ftijh'kxx 

=  lijf"(0)(rk(9kt  +  Xk6kx),rkvk )  +  lijf"{0)(rk6k,rk(vkt  +  \kvkx)) 

+  (0(l)\9kx\\vk\  4-  O(l)|0*lhfc*|)* 

—  0(1)8*  M(t)94(x,  t\  Afc,  4 fik  +  s)  4*  (0(1)8*  M(t)  9$ (x,  1 ;  Xk ,  4 fik  +  £‘))x 
+  kjf  (0)(rk9klvk(lkBrkvkxx  d-  gkx))i 

where  e  >  0  is  small.  By  (4.10)  the  last  term  is 

lijf'/(0)(rk9k,rk(lkBrkvkx  +  gk))x  ~  kjf"(0)(rk9kx,rk(lkBrkvkx  +  gk)) 
=  (0(1)8* (8*  +  M(t)  +  M2(t))93(x,t]  A*,  4/x*))x 
+  0(1)8*  (8*  +  M(t)  +  M2(t))94(x ,  t;  Afc,  4jxfc  +  e). 
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hkt  +  A khkx  -  iHjhte*  =  0{1)6*(8*  +  M(t)  +  M2(t))64(x,t ;  A*,4/*fc  +  e) 

+  (0(1)6*  (6*  +  M(i)  +  M2(*))03(;M;  A*,4/x*  +  e))B. 

Now  for  those  terms  of  III  in  (4.14)  with  k  ^  i  again  we  apply  Lemma  3.4,  while  for  the 
term  k  =  i  we  apply  Lemma  3.2.  We  obtain 


JJ/  =  - 


1  /*  Z*00  1  (x-y-A,(t~T))2 

.  /  /  (t  —  T)  2e  4/iy(t-r) 

yjAnflij  Jo  J-OO 

■  [0(1)6*  M(T)6s/2(y,T-,  Xi,4m)  +  'Y^hk{y,T)]ydydT 

k^i 


=  0(1)6*  M(t) 


(t~r)~ *c  05/2(y, T; A*, 4^) dydr 


+  0(l)tf*(«*  +  M(t)  +  M2(t))  [rp3/2(x,  t;  X{)  +  £  ft*,  t;  A,)] 

k^i 

=  0(1)6*  (6*  +  M(t)  +  M2(t))  fa/2(x,  t;  A,)  +  £  ft*,  t;  A*)] , 

k^i 

where  in  the  last  step  we  have  used  (3.24)  with  a  =  3/2. 

To  estimate  IV  and  VI  in  (4.14)  we  use  Lemmas  3.2  and  3.3.  By  (3.1),  (4.18),  (4.8), 
(4.24),  (3.10)-(3.12),  for  small  e, 

(4.29) 


IV  +  VI  = 


<**>£/ 


_  (x-y-Aj(t-T-))2 

(t-T)~le  (<*«+■><«-■>  £*(£*+ Af(r)) 


Y  r;  Afc,  4/^fc  +  e)  dydr 


=  0(1)6* (6*  +  M(t))<  d2-e(x,t;Xi,fi*)  +  YJ[^(x,t]Xk,ti*) 

'  k^i 

+  I®  -  Aj(t  +  l)!"*1!®  -  Afc(t  +  1)1  2 

•  char  { min  (A*,  Afc)(t  +  1)  +  if  vT+T  <  a;  <  max(A;,  A&)(t  +  1)  -  Ky/t  +  l}]  | 
=  0(1)6*  (6*  +  M(t))  [^3/2(®»  Ai)  +  $(&,  t;  A*)] , 

k^i 

where  /a*,  K  >  0  are  constants. 

To  estimate  V  in  (4.14)  we  use  Lemmas  3.5  and  3.6.  Notice  that  by  (4.8) 


$(x,  t ;  A)  =  0(l)^3/2(®,  L  A). 
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By  (4.18),  (4.8),  (3.33)  and  (3.37),  for  small  e, 


(4.30) 
V  = 


°*u 


{t  -  t)~ 


M2(r)(r  +  1)  ii'Y^Tl>3/2(y,T-,\k)dydT 


fc  =  l 


=  0(l)M2(t)l  (<  +  l)“?+£^3/2  (Mi A  i)  +  [(t  +  1)  h3/2  (Mi  Afc) 

t  k^i 

+  |x  —  A,(t  +  1)1  l\x  —  \k{t  +  1)1  2 

•  char{min(Ai, \k){t  +  1)  +  KVt  +  1  <  x  <  max(Aj, Afc)(i  +  1)  -  K\/t+  l}] 
=  0(l)M2(t)  [^3/2(s,  t;  Aj)  +  Y,  J(x’ f'- A*)]  - 

k^i 


where  K  >  0  is  a  constant. 

Equations  (4.13),  (4.22),  (4.27)-(4.30)  imply  that  for  all  -oo  <  x  <  oo,  t  >  0,  i  =  1, 
(4.31)  k(z,t)|  <C(.5’'  +  5*M(t)  +  M2(t))[^3/2k,i;Ai)  +  S^k.*;k)]- 

k^i 


Next  we  estimate  v*x.  By  (4.19),  /  in  (4.16)  is 

/°°  3  ( x-y-Xji )2  3 

r*c  (4^+e)t  %-(y)dt/  =  0(i)(t  +  i)  ^ 

-oo 


for  t  >  1,  which  is  obviously  also  true  for  t  <  1.  For  II  in  (4.16),  let  hk  be  defined  by 
(4.23).  Then  besides  (4.25)  and  (4.26)  we  have  by  (3.1) 


hkx(x,t)  =  0(l)\6kx\\6k\  =0(l)6*263(x,t-:\k,  2Mfe+£), 

(hfcj  +  ftijh'kxx'jx 

=  Y  [2Zy/"(O)(r*0*,rfc(0*  +  A*M)  +  O(l)(|0fc||«te|  +  l^x|2)] 
=  (O(1)<S*204(:M;  A*,  2/jfc+ <•))*. 


Apply  Lemma  3.4  to  II.  Similar  to  (4.27)  we  obtain 


11=  - 


=£  f  [ 

KVa 


4V^J 


t  /*0O  1  (a:  —  y— Aj(t  —  r  ) ) 2 

(t-r)_Je  4“o«-T>  hkyv{y,T)dydr 

OO 


—  0(l)<5*2(i  +  1)  2  ['03/2(**'i Ai)  4"  ^  ^ Afc)] 

i 

=  0(l)(5*2(t  +  l)-*, 


(4.33) 
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using  (4.8).  By  (3.1)  and  (4.18),  rewrite  III  in  (4.16)  as 

't/2  pOO 


III  =  0(1) 


+ 


/*E/Z  poo 

/  /  {t 

JO  J  —  co 


r)  2e 


(x-y-X^t-T)) 
(477- 


+,,<-T)  +r2)05/2(2/,r;Alt,4W) 


k=  1 


+  M2(r)(l  +  r)  ^3/2(2/,  t;  Afc)]  dydr 


/■(  /»oo 

0(1)/  /  (t  -  r)_1e_  '<‘^<-.5  £[(«*M(r)  +  «*a)flr/a(»,r;  Afc,4W  +  e) 

Jt/2J-oo  ^ 

+  M2(r)(l  4-  r)~ 5 ?/j3/2 (j/,  r;  A*)]  rfj/dr 

=  0(l)(6*2  +  6*M{t)  +  M2(t))(t  +  l)-i  /  /  (t  -  r)-ie-  (««„+i)(. 

Jo  J —co 


(X-y-Xi{t-T)y 


r)V 


'  y^[^5/2(g/>T;  Afc,4^fc  +  £)  +  (1+t)  ^3/2(2/,  t;  Afc)]  dydr. 

fc=i 

Apply  Lemmas  3.2,  3.3,  3.5  and  3.6  to  the  right-hand  side.  We  obtain 
(4.34)  777  =  0(1)(6*2  +  6*M(t)  +  M2  (*))(*  +  l)-l. 

Equations  (4.15),  (4.32)-(4.34)  together  imply  that  for  -00  <  x  <  00,  t  >  0,  i  =  1,  . . .  ,  5, 
(4-35)  \vix(x,  *)|  <  C{8*  +  8*M(t)  +  M2(t))(t  +  l)“i . 


By  the  definition  of  Af(t)  in  (4.17),  inequalities  (4.31)  and  (4.35)  imply  that 
M(t)  <  C(8*  +  8*M(t)  +  M2(t)). 


If  M(t)  and  <5*  are  small,  we  have  M(t)  <  C8*.  By  continuity,  M(t )  <  C8*  if  8*  is  small. 
With  (4.18)  we  thus  have  proved  (4.7). 


5.  Green’s  Function  for  a  2  x  2  Linear  System 


In  this  and  next  sections  we  study  the  Green’s  function  for  the  linearization  of  the 
viscous  conservation  laws  (2.1).  In  this  section  we  explore  a  typical  2x2  system,  the  p- 
system.  This  simplest  model  provides  us  a  clear  picture  how  the  behavior  of  the  eigenvalues 
of  — /'( 0)  +£B{ 0)  for  small  £  and  for  large  £  plays  an  important  role  in  Green’s  function. 
Most  estimates  in  this  section  apply  to  general  n  x  n  systems. 

The  nonlinear  p-system  describes  the  compressible,  isentropic,  viscous  1-D  flow  in  La- 
grangean  coordinates.  The  linearization  of  the  p-system  is 

f  wu  ~  W2x  =  0 

\  W2t  +  p'lVix  =  pW2xx,  <  x  <  °°5  t  >  0, 

where  p'  <  0  and  p  >  0  are  constants.  Physically  wx  and  w2  represent  the  perturbations  of 
specific  volume  and  velocity  about  a  constant  state,  respectively,  while  p'  is  the  derivative 
of  pressure  and  p>  the  viscosity  at  the  constant  state.  Set 


Equation  (5.1)  can  be  written  as 


(5.3)  ut  +  Aux  =  Buxx. 

The  Green’s  function  G{x,t)  of  (5.1)  or  (5.3)  is  the  solution  matrix  to  (5.3)  satisfying  the 
initial  condition 

(5.4)  G(x,  0)  =  8{x)I, 

where  <5(x)  is  the  Dirac  6-function  and  /  is  the  2  x  2  identity  matrix. 

Our  purpose  is  to  obtain  a  pointwise  estimate  for  G.  Here  the  difficulty  lies  in  the 
fact  that  B  is  singular  and  hence  system  (5.3)  is  not  diagonalizable.  To  cope  with  the 
difficulty,  we  approximate  (5.3)  by  a  diagonalizable  system,  whose  Green  s  function  is  a 
combination  of  heat  kernels,  and  then  estimate  the  remainder  by  studying  the  Fourier  and 
inverse  Fourier  transforms.  The  result  in  this  section  is  due  to  Zeng  [Ze].  Here,  however, 
we  give  a  simpler  proof,  which  can  be  extended  tonxn  systems  in  the  next  section. 
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(5.6)  E{z)  =  —A  +  zB,  z  =  i£. 

With  (5.2),  it  is  straightforward  to  compute  the  spectral  representation 

(5.7)  E(z)  =  A -(z)P-(z)  +  A  +(z)P+(z), 

where  A±(z)  are  eigenvalues  of  E(z ),  and  P±(z)  are  corresponding  eigenprojections,  satis¬ 
fying 

(5.8)  P-+P+  =  /,  P±=P±,  P-P+  =P+P~  =0. 

In  fact,  let  l±(z)  and  r±(z)  be  the  left  and  right  eigenvectors  satisfying 


l±(z)E(z)  =  A  ±(z)l±(z),  E(z)r±(z)  =  A  ±{z)r±(z), 

l+r+  —  l_r  —  =  1,  /+r_  =  Z_r+  =  0. 

Then  we  can  set  P±(z)  =  r±(z)l±(z ).  The  result  is 


A±w  =  Y  =fpW, 


From  (5.7)  and  (5.8)  it  is  clear  that 
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Thus  the  solution  to  (5.5)  is 

(5.10)  G(£,t)  =  +  eiiX+w)tP+(iO  = 

Perform  the  inverse  Fourier  transform 

1  f°° 

F_1  (<?)  =  7T  Pr- v-  /  tyxi 

Z'K  J -oo 

to  G.  We  obtain 

1  fOO 

(5.11)  G{x,t)  =  ±  [e^-WPyiO  +  e^^PyiOh^dt;. 

^7r  J- oo 

From  (5.9),  we  have 

*^±(*0  =  -|e2  T  i£\J-p'  ~  (y  )  • 

Since  (i  >  0  and  p'  <  0,  we  have 

Lemma  5.1.  For  real  £  ^  0,  Re{f£A±(z£)}  <  0. 

Again  from  (5.9)  we  notice  that  A±(z)  and  P±{z)  are  branches  of  algebraic  functions  of 
a  complex  variable  2,  having  common  branch  points  yf^i  of  order  one.  At  each 

branch  point  A±(z)  are  continuous  while  P±(z)  have  a  pole  there.  We  also  notice  that 
when  z  describes  a  small  circle  around  z+  or  z_,  {A+(z),  A_  (2)}  and  {P+(z),  P-(z)}  both 
will  undergo  the  transposition  of  themselves  after  analytic  continuation.  This  implies  that 
G(£,t)  in  (5.10)  is  single-valued  in  £.  It  is  easy  to  show  that  d(£,t)  is  bounded  near  the 
isolated  singularities  =  z± : 

\G(^t)\  <  <  e|€|(|A|+|€||B|)t 

_►  as  ^  ^  z±m 


Hence  we  have  proved  the  following  lemma: 

Lemma  5.2.  in  (5.10)  is  an  entire  function  of  f. 

we  now  approximate  system  (5.3)  by 


(5.12) 


u>i  +  Aux  —  P 
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where 


b-(§  })• 


System  (5.12)  is  a  diagonalizable  system,  and  its  Green’s  function  G*{x,t)  can  be  found 
by  diagonalization.  Here,  however,  we  go  through  the  analysis  by  Fourier  transform  so  as 
to  compare  with  the  estimates  for  G.  From  (5.12),  G*  satisfies 


G*  =  %£&&)& 

6*«.o)  =  J, 


where 


E*(z)  =  -A  +  zB*. 


With  (5.2)  and  (5.13),  we  have  the  spectral  representation 


E*(z)  =  \*_(z)Pl  +  \:(z)P* 


where  the  eigenvalues  and  eigenprojections  PI  of  E*  are 


-  in — 1 

2  +2v/=F 

Z  1 
"•“2  2 


Similar  to  (5.10),  the  solution  to  (5.14)  is 


=  e^A1^P*  +  c*A+(*)*p< 


The  inverse  Fourier  transform  gives  us 
(5.17) 


]  i+tv-p' 


1  1 

2 

^  I 

2  2 


-g  2/xt 


1  1 

2  2V^P 

I 

2  2 


The  Green’s  function  G  is  basically  determined  by  the  behavior  of  A ±(z)  and  P±(; 
when  z  is  small  and  when  z  is  large. 
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Lemma  5.3.  The  eigenvalues  A ±(z)  and  eigenprojections  P±(z)  of  E{z),  given  by  (5.9), 
are  holomorphic  at  the  origin,  with  the  property 


(5.18) 


A±(0)  T 

p±(  0)  =  P±, 


where  A^(«)  and  PJ  are  eigenvalues  and  eigenprojections  ofE*(z),  given  by  (5.15).  More¬ 
over,  A±(0)  and  A'±(0)  are  real  with  A'±(0)  >  0. 

The  proof  of  Lemma  5.3  is  a  simple  Taylor  expansion  of  (5.9). 

Lemma  5.4.  For  large  z,  p_(z)  and  ezX+(z)tP+(z),  or  the  other  way  around,  have 

the  expressions 


(5.19  a) 


Q(l)g(c-i2:24-Co2:+0(l))f 


and 

(5.196)  +  e~at[z~1CM  +  0(z~2)(l  +te°*tz  ^)] 

respectively,  where  c_i  and  a  are  positive  constants,  c0  is  a  real  constant,  and  Cm  a 
constant  matrix. 


Proof.  For  large  z,  p(z)  in  (5.9)  has  the  expression 


(5.20) 


p{z) 


0(z~% 


where  the  sign  on  the  right-hand  side  depends  on  how  we  define  the  square  root  of  a 
complex  number.  Suppose  that  we  have  the  sign.  Then  by  (5.9)  and  (5.20), 


\+{z) 


—  +  0(z~3), 

pz 

1  \ 


ju.z(l-f  0(z' 


1  -r  1 

2  +  2(l+0(z-2)) 


(\  ±  |(1  +  0{z-2))  T^(l  +  0{z  2))  \ 
(  ±£(l  +  0(*-2))  |=F§(1  +  0{z-2))) 


(5.21) 
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P_(*)  =  0(1), 


0(z~2). 


Prom  (5.21)  and  (5.22)  clearly  ezX-^tP_(z)  has  the  first  expression  in  (5.19),  where  c_, 
p  >  0,  Co  =  0,  and 


(5.23)  e*^P+(z)  =  e^  (  J  °Q) (e°^ -1  )P+(z)  +  (z~'CM  +  0(z~>)) 


where  CM  is  the  second  matrix  in  (5.22).  Notice  that  \ez  -  1|  <  \z\e^.  The  right-hand 
side  of  (5.23)  gives  the  second  expression  in  (5.19)  with  a  =  -p'  /p  >  0. 

If  we  have  the  sign  in  (5.20),  interchanging  the  subscripts  and  for  A  and 
P  gives  us  the  desired  result.  □ 

Equation  (5.18)  suggests  that  G*  is  the  leading  term  in  G,  while  (5.19)  suggests  that  G 
contains  a  ^-function  which  decays  exponentially  in  t.  Set 


R(x,t)  =  G(x,t)  -  S(x)  ^ 


Lemma  5.5.  For  all  -oo  <  x  <  oo,  t  >  0,  we  have 


(5.25)  |i?(a:,t)|  <  C(t  +  1)  1  e  i,+>ct<0),)  +  e~ l*+xJt(0),>  1  +  ct~^e~t/c  +  Ce~t/c\x\ 


Proof.  From  (5.24),  (5.11)  and  (5.16), 


R{x,t)  =  T  J  |e<«i-«)tp_(i^)  +  e‘£A+(<£)tp+(^)_e«AI(«{)ip. 


ei(Kmtp*  -  eT  (  *  q)  }ete?  d£ 


—  Rl  +  i?2  +  R3 
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where  for  e  >  0  small  and  N  >  0  large, 

(5.26) 

R.=—  f  le^-Wp-W  +  e^+WP+te)  -eiix~mtPl  -  eiiX+{i()tP+}eix(  d^, 

2w  y_e 

eixidt, 


=  i-(  [  N  +  rWe^-^P-iiO  +  e^^P+iiO-e^ 

\ J  — oo  J P 


R2 

R3  =  ~  — 


pLl  I  l  0  \  |  ^ 


0  0 


= _ L  (y  E  +  y°°j  {ei?v-(i?)tP!  +  ei?Ai(if)tP;}e“?  d£ 

+  i-  (y_E  +  j”')  {c«*-ra>‘P_(*0  +  « 

To  estimate  i?i,  we  consider  each  mode  separately: 

(5.27a)  #i  —  #i,-  +  #i,+j 

where 

(5.276)  flllT  =  ±f  {c^«i‘PT«)  -  de 

Consider  Ri  _.  From  (5.18), 

(5.28)  i?!  _  =  A  f  e»«>-(0X-«2A'_(0)t[eou3t)p_(jC)  _  P*]c**«de 

v  J  27t  y_£ 

=  A  [e  ei5A-(o)«-?2A'_(o)t+ix«  [o«3t)e°K',()  +  0(0]  dO 

2?r  y_e 

Since  the  integrand  is  holomorphic  at  the  origin  by  Lemma  5.3,  we  move  the  path  of 
integration  to 

(5.29)  o’-  =  &(-£,£,  c), 


where 

cr(a, 6, c)  ={f  |  Ref  =  a,  Imf  is  from  0  to  c} 

n  { f  I  Imf  =  c,  Ref  is  from  a  to  b } 
fl  {  f  |  Re  f  =  6,  Im  f  is  from  c  to  0} 
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a(  a,b,c ) 


Im  £  =  c 


R 


Figure  5.1 

as  shown  in  Figure  5.1,  and  c  will  be  given  later.  Let  f  =  (  +  irj,  (  and  r)  are  real.  Then 

(5.30)  Re{i£A_(0)<  -  £2A'_(0)t  +  ix£}  =  -V(x  +  A_(0 )t)  -  (2Y_(0)t  +  rf\'_(0)t. 

Here  we  have  noticed  that  A_(0)  and  A'_(0)  are  real  by  Lemma  5.3.  Also  notice  that 
A'_  (0)  >  0. 

For  IFfflT  ^  let  c  =  in  (5.29).  Then  by  (5.28)  and  (5.30)  we  have 


2A'_(0R 

.(Q)t 
4A'_(0)t 


.  O(|C|3t+^-(0)‘|3) 


0(  |ci3t  +  jjL+M°)*l3 

c 


+  o(lC  l  + 


dc 


|x+A_(0)t| 


(5.31)  +  -  f  2i'<0>1  e-’'|x+A-(0)‘|-s2A'_(°)(+r)2A'_(0)t|-Q(£3j)eO(£3t)  _|_0(e)l  dr] 

JO 


* CL 


e  (x+A-  (o)t)2 

8A'_(0)f  2S  A-(UJ£ 

3  ^2  \  / 


ICI  + 


|x  +  A_(0)£| 


dC 


+  Cc-^A-(°)*[e2teCeJt  +  l]es 

(x+A_(0)i)2 


<  C(t  +  l)_1e_~c;  +  Ce~t,c. 

For  -  le’  let  c  =  §  sign(a;  +  A_(0)t)  in  (5.29).  Again  by  (5.28)  and  (5.30), 

\Ri  _|  <cj  e-SI^-(0)(|-C2A'_(°)(+iE2A'_(o)t^3tec£3t  +  gj  d(- 


-77|x+A_  (0)i|  — e2  A'_  (0)<+772  A'_  (0 )t  _Ce3t 


£  te 1  +£  \  dr) 


<  C 


+cf 

Je  e-i«*X'_(0)f-C»X'_<0)t  [e3 uc*t  +  e]  dC  +  Ce-t/c 


<  Ce~t/c. 
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Thus  we  have  proved  that  for  all  — oo  <  x  <oo,  £  >  0, 

,  (x  +  A_(0)t)2  ,n 

(5.32)  <  C(t  +  l)-1e  ^  +  Ce  t/c. 


Similarly, 

-  (x  +  A+(0)t)2  ./r, 

(5.33)  \Ri,+\  <  C(t  + 1)-^  ^  T  Ce  . 


Inequalities  (5.32)  and  (5.33),  together  with  (5.27a),  imply  that  \Ri\  is  bounded  by  the 
right-hand  side  of  (5.25). 

Next  we  consider  R2.  From  (5.26)  and  Lemma  5.4, 

—  N 


|B,I  - 


,(— c-i£  +icoC+G(l))t 


+  e~at  [-CM  +  0(r2)(l  +  te0(t£_1))]  } eix(  <%  , 

K 

where  c_  1  and  a  are  positive  constants,  c0  is  a  real  constant,  and  CM  a  constant  matrix. 
Thus 


\R2\<C^J  +J^  )e_^t/c  di  +  Ce  at 

— (/->/>- 

<Cr!e-*'°  +  C«-|2/”!i3j£ 


<  Ct~ie~t/c  -f  Ce~t/c\x\, 


which  is  bounded  by  the  right-hand  side  of  (5.25). 

To  estimate  R3,  notice  that  by  (5.18)  Re{z£A±(*£)}  =  ~^±(0)£2  with  ^±(0)  >  0  for  real 
£.  Also  notice  that  Re{i£A±(i£)}  <  —1/C  for  —  N  <  £  <  —£  and  e  <  $  <  N  by  Lemma  5.1 
and  by  the  fact  that  A±  are  continuous  functions.  However,  as  we  have  pointed  out  before 
P±(i£)  have  a  pole  at  the  algebraic  singularities  =  \z±  =  ■  Let  I±  denote  the 

small  intervals  [£±  -e,£±  4*  e]*  From  (5.26), 

(5.34)  |i73|<c(/E  +  /  ){e-v-(0«2t  +  e-A+(0)?2t}dC  +  Ce-(/c 

+c\{L+L ) 
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By  (5.10)  and  Lemma  5.2,  the  integrand  of  the  last  integral  in  (5.34)  is  analytic.  Hence 
we  can  move  the  path  of  integration  from  I±  to 

CT±  =  o'tefc  -  £,  £±  +  e,  £ sign(x)), 
see  Figure  5.1.  The  last  term  in  (5.34)  is  then  bounded  by 

+/  j{eReWA-(<f)}i|P_(^)| +eRe{<fA+W»t|P+(^)|}e-:c,)d^|. 

Take  Puiseux  expansion  for  A ±(i£)  at  £_  and  notice  that  A±  are  continuous  there.  The 
integral  on  cr_  is  bounded  by 

C  f  {eRe«-X-<if-)>t+°(£l/2)‘|P_(^)|  +  eRe«-A+(if-)}(+°(£l/2)t|P+(^)|}  d\t\ 

J  (7  — 

<  Ce~^c. 


Similarly  the  intergral  on  cr+  is  also  bounded  by  Ce  t/c ,  it  is  clear  that  the  first  term 
in  (5.34)  is  bounded  by  Ct~ ^2e~^c .  Thus  \R3\  is  bounded  by  the  right-hand  side  of 
(5.25).  □ 


Lemma  5.6.  Let  K  >  0  be 


(5.35)  \R(x,t) |  <  Ct~*(  e~&  H-e~^+‘aT^  +e~— 


For  -oo  <x  <  oo,  t>  0,  if  \x\/t  >  K,  then 

)• 


Proof.  From  (5.24), 


(5.36) 


R(x,t)  =  R4-G*{x,t), 


where 


fl4  =  G(*,0-cVtf(*)(j  J)  =±-J^G(U)-e^  (j 


The  integrand  on  the  right-hand  side  is  an  entire  function  by  Lemma  5.2.  Hence  we  can 
move  the  path  of  integration.  Set  (  =  f  +  irj  with  real  f  and  rj.  let  v  >  0  be  a  constant 
such  that  vc- \  <  where  c_j  is  defined  in  Lemma  5.4.  Let 
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see  Figure  5.1.  For  \x\/t  >K,K  large,  from  (5.10)  and  Lemma  5.4  we  have 
(5.37)  \R4\  =  lim  [  {0(l)e<-c-’«2+ic°«+o(1))t 

27T  N^oo  Ja 

+  e~at \-CM  +  0(r2)(l  +  te°(tr '>)]  }eix «  d£, 

K 

where  c_i,  a  >  0,  co  is  real,  and  Cm  an  constant  matrix.  The  integrand  on  the  line 
segments  £  =  d=iV  is  bounded  by 

c^-c-tNH+c-^/t+omzi+t)  +  e-°t[i_  +  _L_(i  +  te«/")]  j, 

which  goes  to  zero  as  IV  — ►  oo  for  fixed  x  and  t.  Therefore,  in  (5.37)  the  integral  on  the 
line  segments  (  =  ±N  goes  to  zero  as  N  — >  oo.  Equation  (5.37)  becomes 


f?4  = 


in 

2tt  J_  ocl 


+  e' 


^(l)e(-c_i(C4-i^/U2+icn(C+^^/*)+0(1))i 

1  _  -  /  1 


■at 


f£  —  vx/t 


Cm  +  O 


C2  4- 


(1  +  te°W) 


eix(C+ivx/t) 


Notice  that  uc-i  <  h  and  K  is  large.  We  have 


\R*\  <  C 


/: 


■c_iC2t+c_it/2^-coi/x+0(l)f-i/^r^  _j_  Qe 


-at-v^r 


£  r00  e^_ 

J- oo  -  vx/t 


d( 


s: 


■at 


+  CJ 

J  —  cx 

<  C'.-".'.'*'1'*1  J 


-c_iC2i 


d£  +  (7e 


—  ckZ- 


-oo 

1 


+  T 


i/x/t 


i 


d( 


5ix£ 


o  [X 
+  Cc-t/c“^ 


t  \(2  +  {vx/t)2  '  i£  C2  +  {vx/t)2 

<  (l  +  JfU  (f)2w)  +  Ce~t,C~ ** 

1  __  T2 

<  Ct  2  e  c-t . 


Thus  | /?4 1  is  bounded  by  the  right-hand  side  of  (5.35).  From  (5.17)  it  is  clear  that  G  (x,t) 
is  also  bounded  by  the  right-hand  side  of  (5.35).  Equation  (5.36)  then  implies  (5.35).  □ 

we  now  have  the  main  theorem  of  this  section. 
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Theorem  5.7.  [Ze]  For  all  — oo  <  x  <  oo,  t  >  0,  the  Green’s  function  G  of  the  p-system 
(5.1)  has  the  property 

1  1  /  (x+t^-p')2  (x  —  t  sZ-p'f2  \ 

(5.38)  G(x,t)  =  G*(x,  t)  +  0(l)(t  +  l)-*r  *  (  e"  c«  ;  +  e~V  ct  ;  ) 


+  e  ^(a:) 


o  or 


where  G*  is  the  Green’s  function  of  the  uniform  parabolic  system  (5.12),  given  by  (5.17), 
and  6  is  the  Dirac  6 -function. 

Proof.  Fix  a  large  K  in  Lemma  5.6.  If  \x\/t  >  K ,  by  (5.35)  we  have  for  a  constant  C\  >  0, 

(5.39) 

,  /  (x+ty/-p')2  (x-t\/^p')2  \ 


R(x,t)\<Ct  2  e  cP+e  cif  + 


<Ct~*\e 


(x+fyr^7)2 _ ^ 

,  2C1(l  +  v/rp7//fV2t  2Cl 


(x+ty^v)2  (x-tv^)2  ^  (k-V^p7)' 

e  2clt  +  e  2C11  ]e-“  s CT 


i  i  /  (-+*v^v)2  (x-t/v)' 

<  C(i  +  l)’ir  i  f  e-1 — +  e"1 — cr^- 


If  \x\/t  <  K,  by  (5.25)  we  have 

(5.40) 

/  (a=  +  *VCF)2  r*v/V)2  \ 

|H(x,t)|  <C(f  +  l)-‘  e  ci*  +  e~  ci‘  )  +Ct~ie~t/Cl  +Ce~t/Cl\x\ 


<  C(t  +  i)-ir  i  (  e-v  c*  ;  +  e-k  ct  ; 


t  (*+ty^)2 

Cl  2C1(  K  + 


since  AT(0)  =  ±y/IIpr  by  (5.9)  and 


e“ c/01  <  e 


Equation  (5.24),  together  with  (5.39)  and  (5.40),  gives  us  (5.38).  □ 

We  have  similar  result  for  the  derivatives  of  G. 

Theorem  5.8.  For  all  —  oo  <  x  <  oo,  t  >  0,  the  Green’s  function  G  of  the  p-system  (5.1) 
has  the  property 

$  1  £4-1  {  (x  +  t  y7—  p/ )  ( x  —  ty/  ~p'\2  \ 

(5.41)  jrjGfat)  =  r  G*(x,t)  +  0(l)(i  +  lrit-V  (  e-1  c  j  +  e-L--c.  1  ) 


+  e^tJ2S(l-i)(x)Qj, 
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where  l  >0  is  an  integer,  G*  given  by  (5.1 7),  and  Qj,  j  =  0,  . . .  ,  l,  are  2x2  polynomial 
matrices  in  t  with  degrees  not  more  than  j.  Especially , 

*-(!  S)- 

Proof.  We  have  from  (5.10)  and  properties  of  the  Fourier  transform, 

(5.42) 

Similar  to  Lemma  5.4,  one  of  zlezX~^tP-(z)  and  zlezX+^tP^(z)  has  the  expression 

O(l)^e(c-lz2+C0Z+o(1))* 
for  large  2,  while  the  other  has  the  expression 

J2  z‘~jQj  +  e~at  V~lQi+ 1  +  0(z~2)(  1  +  •  •  •  +  tl+1)  +  0{z-2tl+2)e0(tz~">] , 

3=0 

where  c_i  and  a  are  positive  constants,  Co  is  a  real  constant,  and  Qj,  j  =0,  ...  ,  l  +  1, 
are  2x2  polynomial  matrices  in  t  with  degrees  not  more  than  j.  Especially,  Qq  is  given 
explicitly  as  in  the  theorem.  Set 

(5.43)  RW(x,t)  =  J^G(x,t)  -  £iG*(x,t)  -  e£  ^S^(x)Qj. 

3=0 

As  in  the  proof  of  Lemma  5.5,  we  write  R ^  by  (5.10)  and  (5.16)  as 

R(i)  =  —  [  J(j£)'  [e«A-«)tp_(^)  +  ei«A+(iO <p+(^)  _  e*£Al(t ()tp* 

in  j_  oo  t 

(5-44)  -  e«A;«)‘p*]  -  tir  Y&r'Qi 

j=o  ' 

=  r 5° + 4° + 4'\ 

where  for  £  >  0  small  and  N  >  0  large, 

p(0  =  —  f  (i£)'{e**-(i?)t.P_(tf)  +  e*«A+(»Ofp+(^)  _  e^-M^p* 

2tt  i_£ 
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R*)  =  ^U-oo  +L  ){(*OieieX-(<e)‘P-«)  +  «)Iei^+W)tP+«) 

Rs)  =  -^(f  + 1  ^j(iOl{eiiX-(i0tPl  +eitx+WtP;}eixSdt 

+  ~ '{f  N  +Je  ) +  e*^+(i5)fp+(^)}e“? d4 

f_ 

3=0 

The  estimating  for  and  is  almost  the  same  as  for  R2  and  R3  in  the  proof  of  Lemma 
5.5,  where  we  notice  that 

/oo 

ile~^t/c  di  =  0{l)t~^  e~t/c . 

The  result  is 


(5.45)  +  \RW\  <  Ct-l-¥e-Pc  +  Ce~t!c\x\. 

To  estimate  R[l\  we  have  —  R[r>_  +  R[1\,  where 

R{%  =  L  f  (it)1  {e*x*(ili)tPT(iO  -e‘€*i(‘€)«p*}ctee  df. 


Compared  with  the  definition  of  R1jT  in  (5.27b),  we  see  that  the  integrand  for  iq  ^  has  an 
additional  factor  Following  the  same  procedure  for  R\y-,  we  have  for  <  k£j 


lM?-l  <cf 


(=c  +  A_  (0)t)2  , 

Si'  (0)t  2^  A-iU)( 


(+1  ,  I*  +  A_(0)t|z+1 


ici  + 


^+i 


d(  +  Ce"*/c 


i  +  2  O  +  *_(0)t)2 

<  C(t  +  1)  2" e - ci - +  Ce~t/c, 


cf.  (5.31).  For  JSM  >  l£j  has  the  same  bound  Ce  t/c  as  Similarly  we 

can  estimate  R®+-  Together  with  (5.44)  and  (5.45),  we  thus  have 


(5.46)  |i?«|  <C(t  +  l)-^[e 


(z  +  *_(0)t)‘ 
Ct 


+  e 


(x  +  A+(0)t)' 


+  Ct 


2  e 


-*/°  +  Ce-‘/c|*|, 
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which  is  an  analogue  of  (5.25).  Similar  to  (5.35),  we  can  also  prove  that  for  \x\/t  >K,K 
large, 

;  +  ,  /  t2  (x+tx/^)2  (^-tVZP7)2\ 

(5.47)  |H(i)|  <  Ct~~  ( +  e  ™  +e  1. 

Equation  (5.43),  together  with  (5.46)  and  (5.47),  implies  (5.41).  □ 


6.  Green’s  Functions  for  n  x  n  Systems  with  Applications 


We  now  extend  the  result  of  last  section  to  general  n  x  n  hyperbolic- parabolic  linear 
systems.  At  the  end  of  this  section  we  give  applications  to  the  Navier-Stokes  equations 
and  the  equations  of  magnetohydrodynamics. 


6.1.  Green’s  Functions. 

We  consider  the  general  system 

(^•1)  ut  +  Aux  =  Buxx ,  — oo  <  x  <  oo,  t  >  0, 

where  u  =  u(x,t)  is  an  n-vector,  while  A  and  B  are  n  x  n  constant  matrices.  Equation 

(6.1)  can  be  viewed  as  the  linearization  of  (2.1).  Our  basic  assumptions  for  (6.1)  are  the 
following  linear  version  of  Assumptions  2.1  and  2.3. 

Assumption  6.1.  There  exists  a  symmetric  positive  definite  matrix  A0,  such  that  A0A 
is  symmetric,  and  A$B  is  symmetric  semi-positive  definite. 

Assumption  6.2.  Any  eigenvector  of  A  is  not  in  the  null  space  of  B. 

The  purpose  of  this  subsection  is  to  obtain  an  estimate  similar  to  (5.38)  for  the  Green’s 
function  G(x,t)  of  system  (6.1),  which  is  the  solution  matrix  to  (6.1)  satisfying  the  initial 
condition  G(:r,0)  =  8{x)L 

As  we  have  pointed  out  in  Section  2,  Assumption  6.1  implies  that  all  the  eigenvalues  of 
A  are  real,  and  that  A  has  a  complete  set  of  eigenvectors.  We  use  the  same  notations  here 
as  in  Section  2.  Denote  all  the  distinct  eigenvalues  of  A  as  Ai,  A2, . . .  ,  As,  with  multiplicity 
mi,  m2,  . . .  ,  ms,  mi  +  7712  +  •  •  *  +  ms  —  n.  Denote  the  left  and  right  eigenvectors  associated 
with  A i  as  lij  and  respectively,  j  —  1,  ...  ,  m^,  satisfying 

Arij  =  XiT{jy 

lij  A  =  Xikj, 

(6.2) 

hjTi'j'  —  8u>8jj>, 

M7  =  1  j  =  l,...,m*,  f  =  1 
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We  also  use  the  notations 


/  \ 

l'=  (j-1  /’ 

ri  ~  (ril>  •  •  ■  ?  rim i 

(6.3) 

Pi  =  TiU, 

II 

i—1 

y> 

(ll\ 

L=  A 

R=  (n,...,rs). 

\U 

Here  P{  is  the  eigenprojection  of  A  associated  with  A*,  which  is  unique.  Clearly  LR  =  I 
by  (6.2). 

We  now  approximate  (6.1)  by 
(6.4)  Wt  +  Aux  —  B  ’U'xxi 


where 

B*  =RDL,  D  =  diag(Di, . . . ,  D$)> 

(6'5)  Di  =  liBne  RmiXm‘,  i  = 

By  Lemma  2.1,  we  see  that  under  Assumptions  6.1  and  6.2  each  Dt  is  similar  to  a  diagonal 
matrix  with  positive  diagonal  elements.  Hence  there  exist  T,  £  Rm’ xm’ ,  i  =  1,  . . .  ,  s,  such 
that 


Mil 


T~xDiTi  = 


—  Mi 


i  —  1,  .  .  .  ,  S, 


where  %  >  0,  j  =  1,  . . .  ,  ro*,  are  eigenvalues  of  D*.  Denote  the  identity  in 
and  set 


T  =  diag(Ti, . . .  ,TS), 

M  =  diag(Mi,..-,Ms)» 

A  =  diag(Ai/i, . . . ,  A SIS). 


Multiply  (6.4)  by  T~XL  on  the  left.  We  thus  diagonalize  it  as 

(T~1Lu)t  +  A(T_1Lw)x  =  M(T~lLu)xx. 

The  Green’s  function  of  this  system  is  straightforward,  so  is  that  of  (6.4).  We  notice  that 
the  row  vectors  of  are  left  eigenvectors  of  A,  and  the  column  vectors  of  r &  are  right 

eigenvectors,  associated  with  A*.  Denote  them  as  l*j  and  r*j  respectively,  i  =  1,  . . .  ,  s, 
j  =  1,  . . .  ,  rm.  Clearly  and  r*j  satisfy  (6.2).  We  thus  have 
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Lemma  6.1.  The  green’s  function  of  (6.4)  is 


_S  _  Tit  1  (*-A,t)2 

G*(x,t)  =  YY—  -e~  r*/*7- 


where  fori  =  1,  . . .  ,  5,  j  =  1,  . . .  ,  mi,  r*j  and  l*j  are  the  eigenvalues ,  right  eigenvectors 
and  left  eigenvectors  of  A  satisfying  (6.2)  and 


with  the  Hij  >  0. 


Notice  that  Assumption  6.1  also  guarantees  that  all  the  eigenvalues  of  B  are  nonnegative, 
and  that  B  has  a  complete  set  of  eigenvectors.  We  are  now  ready  to  give  the  main  theorem 
of  this  section. 


Theorem  6.2.  Suppose  that  Assumptions  6.1  and  6.2  are  satisfied,  and  that  B  has  the  zero 
eigenvalue  of  multiplicity  m,  0  <  m  <  n.  Let  loj  and  r$j  be  the  left  and  right  eigenvectors, 
respectively,  of  B  associated  with  the  zero  eigenvalue,  satisfying  =  6jk,  j,  k  =  l, 

. ..  ,  m.  Then  for  —00  <  x  <  00,  t  >  0,  the  Green’s  function  G  of  system  (6.1)  has  the 
property 

s  (i-a  t)2  m' 

(6.7)  G(x,t)  =  +  0{l)(t  +  ~  MVk, 

3= 1  k= 1 

where  G*  given  by  (6.6)  is  the  Green’s  function  of  (6.4);  \j,  j  =  1,  ...,  s,  are  the  distinct 
eigenvalues  of  A;  C  >  0  is  a  constant;  vn!  <  m;  for  k  =  1,  . . .  ,  m! ,  >  0  and  j3k  are 

constants,  Vk  are  constant  projections,  satisfying 

VkPk>  =  hkk'Vk,  k,k'  = 

m!  m 

fc=l  j=l 

m' 

Y.PkPk  =  Q0AQ0, 

k= 1 

and  8  is  the  Dirac  8 -function. 

Remark  6.3.  If  B  is  nonsingular,  m!  =  0  and  the  second  summation  in  (6.7)  disappears. 
We  may  say  that  the  ^-functions  in  (6.7)  are  due  to  the  hyperbolic  ingredient  of  the  system. 
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Remark  6.4.  Suppose  that  B  is  singular.  Prom  (6.7)  we  see  that  dx/dt  =  f3k,  k  =  1,  . . .  , 
m',  are  the  directions  in  which  an  initial  jump  will  propagate.  From  (6.8)  the  are  the 
eigenvalues  of  QoAQq  when  restricted  to  the  range  of  Qq.  We  can  also  compute  ak 
and  Vk ,  k  =  1,  ...  ,  m',  explicitly  using  a  reduction  process.  In  fact,  let  1/1 ,  . . .  ,  vr  be 
all  the  nonzero  (and  hence  positive)  eigenvalues  of  B  with  corresponding  eigenprojections 
Qu  ...  ,  Qr  (cf.  Pi  in  (6.3)).  Set 

T  1 

(6.9)  S  =  —Qh 

j= 1  3 

which  is  the  value  at  zero  of  the  reduced  resolvent  of  B  with  respect  to  the  zero  eigenvalue. 
Associated  with  each  distinct  eigenvalue  /3k,  Q0AQ0  has  the  eigenprojection  Qk-  Here  if 
/ 3k  happens  to  be  zero,  Qk  is  taken  as  the  subprojection  QkQo •  Corresponding  to  this 
f3k  there  may  be  several  ak->  which  are  all  the  nonzero  distinct  eigenvalues  of  QkASAQk- 
Each  ctk  then  associates  with  an  eigenprojection  of  QkASAQk ,  which  is  the  corresponding 
Vk-  At  last  m'  is  the  total  number  of  the  ak  so  constructed  for  all  the  distinct  pk.  In  the 
special  case  that  all  the  j3k  are  simple,  we  have 

(6.10)  m!  =  m,  Vk  determined  by  (6.8),  ak  =  tr {ASAVk),  k  —  1,. . .  ,m. 

Remark  6.5.  Prom  Remark  6.4  we  see  that  if  B  is  replaced  by  (xB  in  (6.1)  with  M  >  0, 
then  ak  =  0(m_1)  while  (3k  and  Vk  are  independent  of  /i,  k  =  1,  . . .  ,  m! .  This  can  also 
be  seen  from  G(x,t;/0  =  JC?(a://x, t//x ;1),  which  comes  from  (6.13)  below.  Therefore,  in 
(6.7)  the  ^-functions  disappear  when  fi  — >  0.  We  may  say  that  the  ^-functions  are  due  to 
the  parabolic  nature  in  addition  to  the  hyperbolic  ingredient  (Remark  6.3)  of  the  system. 

The  proof  of  Theorem  6.2  consists  of  a  series  of  lemmas.  First  we  notice  that  under 
Assumption  6.1  we  may  decompose  Aq  as  A$  =  (Aq)2,  where  A q  is  a  symmetric  and 
positive  definite  matrix.  Introduce  a  new  variable  w  —  Ag  u.  Then  w  satisfies  a  symmetric 
system 

wt  +  AlA(Al)~xwx  =  Aq  B(A£  )~lwxx 

since  A0A,  A0B  and  (dj )  ~ 1  are  all  symmetric.  Therefore  it  is  sufficient  to  prove  Theo¬ 
rem  6.2  for  symmetric  systems.  For  the  rest  of  this  subsection  we  assume  that  A  and  B 
are  symmetric  and  Aq  =  I . 
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Perform  the  Fourier  transform  with  respect  to  x  for  G.  From  (6.1)  we  have 


(Gt=  iCEUOG 

1  6(C,o)  =  J, 

where 

(6.11) 

E(z')  =  — A  +  zB . 

The  solution  is 

(6.12) 

G{£,t)  = 

Hence 

(6.13) 

G(x,t)  =  F“1{ei5£(^)t 

From  the  case  of  the  p-system  we  see  that  G  is  determined  by  the  behavior  of  the 
eigenvalues  and  eigenprojections  of  E.  We  now  establish  a  series  of  lemmas  analogous  to 
Lemmas  5. 1-5.4.  The  following  lemma  is  due  to  Shizuta  and  Kawashima. 

Lemma  6.6  [SK].  Under  Assumption  6.1 ,  Assumption  6.2  is  equivalent  to  any  of  the 
following  statements: 

(i)  There  exists  a  K  E  Mnxn,  such  that  KAq  is  skew  symmetric  and  ^(KAqA  + 
AqAK*)  +  AqB  is  positive  definite. 

(ii)  Let  A (z)  be  any  eigenvalue  of  E{z)  defined  by  (6.11).  Then  Re{^A(^)}  <  0  for 
real  £  ^  0. 

(Hi)  There  exists  a  constant  c  >  0,  such  that  Re{if  A(i£)}  <  -c£2/(  1  +  £2)  for  real 
where  A(z)  is  any  eigenvalue  of  E(z). 

Proof,  we  have  assumed  that  A  and  B  are  symmetric  and  Aq  =  I .  For  the  general  case, 
we  can  prove  the  lemma  by  making  use  of  the  decomposition  A0  =  (AJ)2. 

We  first  prove  that  Assumption  6.2  implies  (i).  Define  an  inner  product  {■,  •}  in  ]RnXn 
as 


{X,Y}  =  tv(XYt) 
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for  any  X,  Y  e  RnXn.  Define  a  linear  transform  TA:  RnXn  ->  RnXn  as 

rA{X)  —  ax  -  xa,  vx  e  Rnxn. 

Denote  the  null  space  and  range  of  VA  as  Si  and  S2  respectively.  We  have 

(6.14)  dimSi  +  dimS2  =  dimRnXn. 

For  any  X  €  Si,  AY  —  Y A  €  52, 

{X,  AY  -  YA}  =  tr(Xy*j4  -  XAYf)  =  tr (XYtA  -  AXY1) 

=  tr  (XYtA)  -  tr  {AXY1)  =  0, 

where  we  have  used  the  facts  that  A  is  symmetric,  ^4X  =  XA  since  X  e  Si,  and  tr(ST)  = 
tr(TS).  Thus  Si  _L  S2.  Together  with  (6.14),  we  conclude  that  S2  -  Sj~. 

By  (6.2)  and  (6.3),  we  have  the  spectral  representation 

s 

(6.15)  A  =  Y/XJPj, 

3= 1 

It  is  easy  to  see  that  each  Pj  is  symmetric  since  A  is  symmetric.  If  Y  e  RnXn  commutes 
with  A,  Y  must  commute  with  Pj.  In  fact,  multiply  YA  —  AY  from  the  left  by  P3  and 
from  the  right  by  Pk ,  use  (6.15)  and  notice  that  PjPk  —  $jkPj-  We  arrive  at  A kP jYPk  = 

s 

A jPjYPk,  or  PjYPk  =  0  for  j  ^  k.  This  implies  that  PjY  =  PjY  Y  Pk  =  PjYPj. 

J  J  k=  1 

Similarly  we  obtain  YPj  —  PjYPj  and  therefore  YPj  = PjY . 

Define  a  linear  transform  HA:  Rnxn  — ►  Rnxn  as 

s 

(6.16)  n  A(X)  =  J2P>XPP  VXeR"Xn 

i=i 

For  any  X  €  Rnxn,  Y  €  Su 

S 

{X  -  nA(X),Y}  =  tr (XY*  -  n^(X)y‘)  =  tr(Xy')  -  tr (J^PjXPjY1) 

j= 1 

$  S 

=  tr  (XY*)  -  tr(£  YPjX'Pj)  =  tr  (XY*)  -  tr(£  PjYX'Pj) 

j= 1 

=  tr(xyf)  -  tr(^b  YXtPj)  =  tr(Xy*)  -  tr(yX()  =  0. 
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Hence  X  -  nA(X)  e  Sj~  =  S2.  Especially,  there  exists  a  K  e  Mnxn,  such  that 

B-Ua{B)=AK-KA. 

Set  Kx  =  \{K  +  K1),  K2  =  \(K-  K *).  Then 

B  —  HA(B)  =  A(Ki  +  K2)  —  ( K\  +  K2)A. 

The  left-hand  side  is  symmetric,  so  is  the  right-hand  side,  which  implies  that  AKX~KYA  = 
0.  Thus 

K2A-AK2  +  B  =  IIa(B). 

From  (6.16),  using  the  notations  in  (6.3)  and  (6.5),  we  have  E [A{B)  =  RDL  =  B* .  As¬ 
sumptions  6.1  and  6.2  imply  that  all  the  eigenvalues  of  B*  are  positive,  as  we  have  pointed 
out  before,  hence  II^B)  is  symmetric  and  positive  definite.  Set  K  =  2K2  and  obtain  (i). 

We  next  prove  that  (i)  implies  (iii).  For  real  £,  let  A(i£)  be  any  eigenvalue  of  E(i £)  with 
eigenvector  r(z£): 

(6. 17)  (—A  +  i£B)r  =  Ar. 

Multiplying  (6.17)  from  the  left  by  i^f1  and  taking  real  part,  we  have 

(6.18)  Re{^A}|r|2-h?VBr  =  0, 

where  we  have  used  the  fact  that  f* Sr  is  real  if  S  G  Rnxn  is  symmetric.  Multiplying  (6.17) 
from  the  left  by  £2r*AT,  K  given  in  (i),  and  taking  real  part,  we  have 

(6.19)  -  Re{i£A }r*Kr  +  ^2f*(KA  +  AKl)r  =  Re^^KBr) 

<  K|3|r||tf||Br|  <  ^2|r|2  +  ±-?\Kf\Brf, 

where  £  >  0  and  we  have  used  the  fact  that  ff5r  is  pure  imaginary  is  S  e  MnXn  is  skew 
symmetric.  Multiplying  (6.18)  by  (1  +  £2)  and  (6.19)  by  a  >  0,  addition  of  the  results 
gives  us 

(6.20)  Re{i£A}(l  +  £2)  (|r|2  —  +  £2[f*Br  +  f*^(KA  +  AK*)r]  +£iftBr 
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Clearly  from  (i)  and  Assumption  6.1  that  B  is  semi-positive  definite  we  have 


r'Br  +  \f*{KA  +  AKl)r  >  cx  |r|2, 

Zt 


\Br\2  <  Ci^Br, 


-t 


i+e 


f  Kr 


<  c3a|r|2, 


where  c}  >  0 ,j  =  1,  2,  3,  are  constants.  Applying  these  inequalities  to  (6.20)  and  choosing 
e  =  cua<  min{l,  £},  we  have 


Re{i|A}(l  +  £2) 


**  -i*Kr 


i+e 


Ci  a 


Re{if\}  <  - 


cia/2  £2 
1  +  c3a  1  +  £2 


Hence  (iii)  is  proved  with  c  =  |cia/(l  +  c3a). 

Obviously  (iii)  implies  (ii).  At  last  we  prove  (ii)  implies  Assumption  6.2.  If  Assump¬ 
tion  6.2  were  false,  there  would  be  an  r  €  In,  such  that  r  ^  0,  Br  =  0,  and  Ar  =  A r  for 
some  real  A.  Thus 


(—A  +  zB)r  =  -Ar, 


or  there  would  be  a  A (z)  =  -A.  Then  Re{zf  A(i£)}  =  0  for  anY  real  f •  This  contradicts 

(n).  D 

We  now  state  some  basic  facts  from  complex  analysis  and  perturbation  theory.  The 
reader  is  referred  to  Kato’s  book  [Kt],  Chapter  2,  and  the  references  therein. 

Consider  E(z)  given  by  (6.11)  as  a  matrix- valued  function  of  the  complex  variable 
z.  The  number  s'  of  the  distinct  eigenvalues  of  E(z)  is  constant  if  2  is  not  one  of  the 
exceptional  points,  of  which  there  are  only  a  finite  number  in  the  whole  complex  plane.  In 
each  simply  connected  domain  containing  no  exceptional  points,  the  eigenvalues  of  E(z) 
can  be  expressed  as  s'  holomorphic  functions 


Ai(z),  ^2(2), . . . ,  As'  (z), 

s' 

which  have  constant  multiplicities  mi,  m2,  . . .  ,  7?v  respectively,  X)  =  n*  The 

3= 1 

1  <  j  <  s',  are  branches  of  one  or  several  algebraic  functions,  which  are  everywhere 
continuous  since  the  coefficient  of  the  highest  power  in  the  characteristic  equation  is  1. 
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These  algebraic  functions  will  also  be  denoted  by  A j(z).  An  exceptional  point  z0  is  either 
a  branch  point  of  some  of  the  A j(z)  or  a  regular  point  for  all  of  them;  in  the  latter  case 
the  values  of  some  of  the  different  A j(z)  coincide  at  z  =  z0.  Hence  there  is  always  splitting 
of  eigenvalues  at  (and  only  at)  an  exceptional  point. 

In  general,  E(z)  of  the  form  (6.11)  is  not  diagonalizable,  and  has  the  spectral  represen¬ 
tation  for  any  z  which  is  not  an  exceptional  point: 


The  spectral  representation  of  a  matrix  is  unique.  Here  Pj  and  Dj,  called  the  eigenprojec- 
tion  and  eigennilpotent  for  A j  respectively,  must  be  given  by 


Dj(z)  =  (E(z)-\j(z)I)Pj(z), 


where  Tj(z)  is  a  positive-oriented  circle  enclosing  A j(z)  but  excluding  all  other  A k(z). 

The  eigenprojections  Pj(z )  and  the  eigennilpotents  Dj(z)  for  A j(z)  are  also  holomorphic 
in  each  simply  connected  domain  containing  no  exceptional  points,  being  branches  of  one 
or  several  algebraic  functions  (again  denoted  by  Pj(z)  and  Dj(z)). 

The  algebraic  functions  Pj(z)  and  A j(z)  have  common  branch  points  of  the  same 
order.  Let  zq  be  an  exceptional  point.  When  z  describes  a  small  circle  around  zq , 
{Ai  (^),  A2(^), . . . ,  As/(^)}  and  {P\(z),P2(z), . . .  ,Pa>(z)}  will  undergo  the  same  permuta¬ 
tion  among  themselves  after  analytic  continuation.  We  emphasize  here  that  z$  may  or 
may  not  be  a  branch  point,  A j(z)  is  always  continuous  there,  while  Pj(z )  in  general  has  a 
pole  at  z0  even  when  z0  is  not  a  branch  point.  If  z0  indeed  is  a  branch  point,  we  have  the 
following  lemma  due  to  Butler. 
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Lemma  6.7.  ('[Bu]  or  [Kt]j  If  z  =  z0  is  a  branch  point  of  \j(z)  (and  therefore  also  of 
Pj(z))  of  order  p  -  1  >  1,  then  Pj{z)  has  a  pole  there;  that  is,  the  Laurent  expansion  of 
Pj(z)  in  powers  of  (z  -  z0)1/p  necessarily  contains  negative  powers. 

Proof.  Without  loss  of  generality  we  assume  that  Zo  =  0.  When  z  is  subjected  to  a 
revolution  around  z  =  0,  P,(z)  is  changed  into,  say,  Pj+i{z).  Suppose  that  the  lemma 
were  not  true.  Let 

Pfc(z)  =  P i0)  +  z1/pp{k1)  +  •  •  •  ,  k  =  j,j  +  1 

be  the  Laurent  expansions  of  Pj  and  Pj+i-  We  must  have  Pj+i  =  e,2lrl^pPj  \  l  =  0,  1, 
...  ,  since  P,(z)  becomes  P,+i(z)  after  the  revolution  of  z.  Especially,  =  Pj0)-  On 
the  other  hand,  the  relation  Pj(z)Pj+1(z)  =  0  in  (6.22)  for  z  ^  0  implies  Pf  ] Pf^x  =  0, 
and  the  idempotent  character  of  P,(z)  implies  {Pf'!)  =  Pj  ^  ■  Hence  P-  =  (Hj  )  — 

p(°) =  o,  and  lim  |P,(z)|  =  0.  We  notice  that  any  eigenprojection  is  not  zero.  In 
fact,  the  rank  of  P,(z)  is  rhj.  Thus  (6.22)  implies  that  |Pj(z)|  >  1.  Therefor  we  reach  a 
contradiction,  and  Lemma  6.7  is  proved.  Ol 

We  now  explore  the  \}(z),  P^(z)  and  Dj{z)  under  Assumptions  6.1  and  6.2. 

Lemma  6.8.  Assume  Assumption  6.1.  For  j  =  1 ,  ,  s',  Dj(z)  =  0,  while  Aj(z)  and 

Pj{z)  are  holomorphic  at  the  origin.  At  the  infinity,  the  \j{z)  are  single-valued,  having  at 
most  a  pole  of  order  1,  while  the  Pj  (z)  are  holomorphic  there . 

Proof.  Under  Assumption  6.1,  both  A  and  B  are  symmetric.  Thus  E{z)  defined  by  (6.11) 
is  real  symmetric  for  real  2.  We  then  have  an  orthonormal  set  of  eigenvectors  for  E(z). 
For  any  real  z  which  is  not  an  exceptional  point,  let  the  orthonormal  set  be  rjk{z ),  j  —  1, 

. . .  ,  s',  k  =  1,  . . .  ,  rhj,  where  the  rjk  are  eigenvectors  associated  with  A  j.  We  set 

mj 

(6.23)  j  =  l,...,s', 

k=i 

which  clearly  satisfy  (6.21)  and  (6.22).  The  uniqueness  of  spectral  representation  then 
implies  that  the  P3  and  Dj  in  (6.23)  are  the  ones  under  consideration.  Notice  that  we 
only  have  a  finite  number  of  exceptional  points.  After  analytic  continuation,  Dj(z)  —  0 
on  the  whole  complex  plane.  Meanwhile,  from  (6.23)  \Pj(z)\  =  1  for  real  z  which  is  not 
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an  exceptional  point.  By  Lemma  6.4  the  analytic  functions  A j(z)  and  Pj(z)  can  not  have 
real  branch  points.  Especially,  they  are  single- valued  at  the  origin.  Therefore,  the  A j(z) 
are  holomorphic  there  since  they  are  continuous,  and  the  Pj  (z)  are  also  holomorphic  there 
since  they  are  bounded.  At  the  infinity,  we  write 

E(z)  =  z{--A  +  B)  =  zE(z). 
z 

Then  E(z)  has  eigenvalues  ^A j(z)  and  eigenprojections  Pj(z),  j  —  1,  ...  ,  s' ,  by  (6.21). 
As  we  have  proved,  ~Xj(z)  and  Pj(z)  are  holomorphic  at  the  infinity  0).  Hence  we 
have  completed  the  proof  of  Lemma  6.8.  □ 

By  Lemma  6.8,  the  spectral  representation  (6.21),  (6.22)  becomes 
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where  B*  is  defined  in  (6.5),  Pj  in  (6.3),  and  the  summation  in  (6.28)  is  for  all  the  k  such 
that  Afc(0)  =  -A j,  1  <  k  <  s'. 

Proof.  By  Lemma  6.8  the  A j(z)  and  Pj (z)  are  holomorphic  at  the  origin.  As  z  -+  0  in 
(6.24),  we  obtain  (6.26)  by  (6.11). 

From  (6.5)  and  (6.3)  we  have 

(6.29)  B*  =  rjDjlj  =  PjBPj. 

j=l  j— 1 

From  (6.2)  and  (6.3)  we  have  the  spectral  representation 

5 

(6.30)  A  =  YJ^ipi- 

3=1 

Equations  (6.26),  (6.30)  and  the  uniqueness  of  spectral  representation  imply  that  we  can 
classify  the  Xj (z)  and  Pj{z)  as  the  A jk{z)  and  Pjk{z),  such  that 
—  Ajfc(0)  =  Aj,  j  =  l,...,s,  k  =  l 

nj 

(6.31)  0)  =  Pj,  j  =  l,...,s, 

k= 1 

m~\ - \-ns  =  s'. 

rij 

Multiplying  E(z)  by  £  Pjk(z)  from  the  left  and  from  the  right  and  using  (6.11)  and 

k= 1 

(6.24),  we  have  for  j  =  1,  ...  ,  s, 

(f2pjk(z))(-A+ zb)(j2p^))  =  (EAtw)  (ELwAw)- 

Vfc=i  /  \i=i  '  vfc= 1  7  l  7 

Having  the  Taylor  expansion  for  small  2,  comparing  the  coefficients  for  2  on  both  sides, 

and  using  (6.31)  and  (6.24),  we  have 
+  P]DP] + 

\k=1  7  vz= 1  7 

=  (E^fc(°))(-^A)  +  E^(0)^i(°)  +  (-Aipi)(E^(0))- 
Vi  7  1=1  xi=i  7 

Since  APj  =  Pj  A  =  XjPj,  we  arrive  at 

PjBPj  =  jh\'jk{0)PJk{0), 

k=l 

which  is  (6.28).  Summation  for  j ,  1  <  j  <  s,  together  with  (6.29)  and  (6.31),  gives  (6.27). 

□ 
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Lemma  6.11.  Let  G  (x,t)  be  the  Green’s  function  of  ( 6.4 )•  Under  Assumptions  6.1  and 
6.2,  the  Fourier  transform  G*  of  G*  has  the  expression 


G*(£,t)  =  0) 

3= 1 


where  the  Aj(0)  are  real  and  the  A' (0)  are  positive. 
Proof  Clearly  from  (6.4)  G*(£,t)  =  where 


E*(z)  =  -A  +  zB*  =  £(A,(0)  +  zA'-(0))Pj(0) 


by  (6.26)  and  (6.27).  Notice  that  Pj(0)i\(0)  —  6jhPj{ 0)  by  (6.24).  Expression  (6.32) 
is  then  straightforward.  FYom  (6.26)  and  (6.27)  we  see  that  the  A^O)  and  A'(0)  are 
eigenvalues  of  -A  and  B*  respectively.  The  Aj(0)  are  real  since  A  is  symmetric.  The  A' (0) 
are  positive  by  Assumptions  6.1  and  6.2,  see  the  explanation  following  (6.5).  □ 

Lemma  6.12.  Under  Assumptions  6.1  and  6.2,  for  large  z  we  have 
(6.33) 


ezA  = 


e-{0Jz+«J)tVj  +  c-(ftz+ay)t  Z-1  cMj  +  0(z-2)(  1  + 1  +  t2e°(^-1)) 

\  s'  A  ^  _ f 


3X  }  I  1  <  j  <  m  , 

{  0(l)e^.-^2+c^o,+o(i))t)  m'  <j<  s', 

where  m!  <  m,  m  is  the  multiplicity  of  the  zero  eigenvalue  of  B;  for  k  —  1,  . . .  ,  m!  and 
j  —  m  -j-l,  ...,  s' ,  ak,  0k,  Cj- 1  and  cJ?0  are  real  constants  with  a*., Cj_i  >  0;  for  k  =  1, 

. . .  ,  m  ,  Vk  are  constant  matrices  while  Cm,)c  are  polynomial  matrices  in  t  with  degrees 
not  more  than  1;  the  fik  and  Vk  satisfy  (6.8).  Here  we  have  arranged  A j(z),  j  —  1,  ,  s', 

in  appropriate  order. 

Proof.  By  Lemma  6.8,  as  2  — >  00, 

\j(z)  ~  Cj'—\Z  +  Cjfi  ~f-  Cjti - f-  Cjt2~z  +  •  •  •  , 


pj(z)  =  Pj,o  +  Pj,i-  + 


By  Assumption  6.1,  A j(z)  and  P3  (z)  are  real  for  real  2  since  E(z)  is  real  symmetric. 
Therefore,  all  the  coefficients  in  (6.34)  are  real.  Comparing  the  leading  coefficients  on 
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both  sides  of  (6.24),  using  (6.34),  we  have 


(6.35) 


B  = 

3= 1 

Pj,()Pk,0  =  $jkPj, 0?  ji  ^  1,  .  .  .  ,  <S  , 


E^.o  =  /. 
j=i 


Hence  Cj  _i,  j  =  1,  - . .  ,  s',  are  eigenvalues  of  B,  which  are  nonnegative  by  Assumption  6.1. 
After  rearranging  the  \j{z),  we  have  ckf-i  =  0,  1  <  k  <  m',  and  c^-i  >  0,  mf  + 1  <  j  <  s  , 
with  m'  <m.  Set 


(6.36)  A  =  -Cfc,o5  afc  =  -cfc)i,  Vk  =  Pk,  o,  = 

We  prove  that  a*  >  0.  In  fact,  for  large  real  £  we  have  from  (6.34)  and  Lemma  6.6  (iii), 

ReRA*(i£)}  =  -a*  +  0(d)  fc  = 

where  c  >  0  is  a  constant.  This  implies  that  ak  >  0. 

For  m'  <  j  <  s',  the  second  expression  in  (6.33)  is  straightforward  by  (6.34).  For 
1  <  k  <  mf,  by  (6.34)  and  (6.36), 

ez~XkMtPk(Z)  =  e~{fikX+ak~Ck-^0{^))\Vk  +  Pk,i~  +  O(L)) 

z  z 

=  e-^z+ak)t  {rk  +  ~{ck>2  tVk  +  Pfc.i)  +  0(4)  (!  +  t  +  «2e°(t*'l))}> 

z  z 

which  gives  us  the  first  expression  in  (6.33)  with  Cm, A:  —  ck,2tPk  +  Pk, l- 

Clearly  the  Vk  satisfy  the  first  two  equations  in  (6.8)  by  (6.36),  (6.35)  and  the  uniqueness 
of  spectral  representation.  Notice  that  the  ~A j(z)  and  Pj(z )  are  the  eigenvalues  and 
eigenprojections,  respectively,  of 

(6.37)  -E{z)  =  B+-(-A). 

z  z 

Regard  w  =  l/z  as  a  new  variable  and  compare  (6.37)  with  (6.11).  We  then  apply  (6.28) 
to  (6.37)  to  obtain 

m!  mf 

Qo{-A)Qo  =  Ck,oPk,Q  =  J^(  — Pk)Pk 

k= 1  k= 1 
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using  (6.34)  and  (6.36).  The  third  equation  in  (6.8)  is  proved.  □ 

We  now  have  Lemmas  6.6,  6.8,  6.9,  6.11  and  6.12  analogous  to  Lemmas  5.1-5.4  in  the 
last  section.  Set 

m! 

(6.38)  R{x,t)  =  G(x,t)  -  G*(x,t)  -  -  fat) Vk. 

k=  1 

By  (6.25)  and  (6.32), 


R(x,t)  =  —  /  \Y,e*Xj{i0tPj(iti)  ~  ^e^W^JW^P^O) 


m!  >. 

Y^e~{Pki(+ak)tVk  \eixt  d£. 

fc— 1  > 


With  tiny  modifications  in  the  proofs  for  Lemmas  5.5  and  5.6,  we  obtain 

Lemma  6.13.  Under  Assumptions  6.1  and  6.2,  for  all  -oo  <  x  <  oo,  t  >  0,  we  have 


(6.39)  \R(x,t)\<C(t  +  l)  1^e  +  Ct  1/2e~t/c  +  Ce~t/c  ^  \x  -  j3kt\. 


Lemma  6.14.  Let  K  >  0  be  large.  Under  Assumptions  6.1  and  6.2,  for  all  -oo  <  x  <  oo, 
t  >  0,  if  \x\/t  >  K,  then 

(6.40)  \R(x,t)\  <Cri(e~£  +  £V 

V  j= i 

As  Theorem  5.7  follows  from  Lemmas  5.5  and  5.6,  (6.7)  follows  from  (6.38)-(6.40). 
Theorem  6.2  is  proved.  We  further  have  the  following  analogue  of  Theorem  5.8,  which 
extends  Theorem  6.2  to  higher  derivatives. 

Theorem  6.15.  Under  the  same  assumptions  and  same  notations  as  in  Theorem  6.2,  the 
Green’s  function  G  of  system  (6.1)  has  the  property 

(6'41)  &G{X’ t]  =  +  OUX*  +  !)-'■ ± 

3  =  1 

m!  l 

+  e~akt  H  (X  ~  MQkj, 

k= 1  j=0 


74 


TAI-PING  LIU  AND  YANNI  ZENG 


where  l  >0  is  an  integer,  — oo  <  x  <  oo,  t  >  0,  and  Qkj,  k  =  1,  ,  m! ,  j  =  0,  ...  t  lt 

are  n  x  n  polynomial  matrices  in  t  with  degrees  not  more  than  j.  Especially,  Qk o  =  Vk, 
k  =  1,  ...  ,  m! ,  with  Vk  being  the  same  as  in  Theorem  6.2. 

We  finish  this  subsection  with  a  justification  of  the  computation  of  m' ,  ak  and  Vk , 
which  was  described  in  Remark  6.4.  The  reader  is  referred  to  Kato  s  book  [Kt] 

We  have  taken  m'  as  the  number  of  \\ j(z)  which  splits  from  zero  at  z  =  oo  in 
Lemma  6.12.  Then  some  of  the  pairs  (0k,  &k)  in  (6.36)  may  be  identical.  Suppose  that  in 
(6.7)  we  sum  up  those  terms  corresponding  to  identical  pairs  and  let  m!  be  the  number  of 
the  distinct  ones.  Now  we  use  a  new  notation  m*  to  denote  the  number  of  |A j(z)  splitting 
from  zero  at  z  =  oo.  Set 

m*  i 

(6.42)  P(z)  =  Yi  A(z)  =  P{0)  +  -PW  +  "  •  • 

fc=I 

Then 

(6.43)  P(0)  =  Qo,  ^(1)  =  SAQo  +  QoAS, 

where  S  is  given  by  (6.9).  The  first  equality  is  clear  by  (6.34),  (6.36)  and  (6.8).  To  prove 
the  second  equality,  by  (6.11),  (6.24)  and  (6.42)  we  have 

ra* 

{-A  +  zB)P(z)  =  Y  h(z)Pk{z)- 

k= 1 

Comparing  the  constant  terms  on  both  sides,  using  (6.34)  and  noticing  ck-i  =  0,  k  =  1, 
...  ,  m* ,  we  arrive  at 

m* 

-AP<®  +  BP(1)  =YCk’°Pk’°- 

k= 1 

Clearly  SPk,0  =  SQ0Pk,o  =  0  by  (6.36),  (6.8)  and  (6.9).  Thus 

(6.44)  -SAPW  +  SBPm  =  0. 

Similarly  we  have 


(6.45) 


-P^AS  +  P{1)BS  =  0. 
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Addition  of  (6.44)  and  (6.45)  yields 

(6.46)  SBPW  +  P^BS  =  SAPW  +  pW  AS. 

Notice  that  under  the  notations  in  Remark  6.4,  we  have  the  spectral  representation  B  = 
JC  ujQj  and  thus  SB  =  BS  =  I  —  Q0,  Also  notice  that  p(°)p(i)  +  pi1) p(o)  —  p(i)  by 

comparing  the  coefficients  of  z"1  in  (P(z))2  =  P(z).  Equation  (6.46)  then  gives  us  the 
second  equality  in  (6.43)  by  using  the  first  one. 

Consider  the  matrix- valued  function 

(6.47)  H{z~l)  =  P(z)E(z)P(z)  =  -ff0  +  z^H,  +  0(z~2) 

for  large  z,  where  by  (6.42),  (6.11),  (6.43)  and  (6.9), 

Ho  =  QoAQo, 

(6.48) 

Hi  =  -QoAP(1)  -  PmAQ0  +  pWfiptt)  =  -Q0ASAQo  -  H0AS  -  SAH0. 

m* 

Clearly  H{z~l)  =  ^  A k(z)Pk(z)  by  (6.42)  and  (6.24).  Thus  all  the  nonzero  eigenvalues 
of  H(z~')  are 

M*)  =  cjt,0  +  cMz_1  +  •  ■  •  ,  k=l,...,m*, 

by  (6.34)  and  the  meaning  of  m*.  By  (6.36)  we  may  assume  -01  =  clfi  =  •  •  •  =  Cpfi  ±  Cj  0, 

p  <  j  <  m* .  We  can  apply  (6.28)  to  H{z~1)  even  it  has  a  higher  order  term  in  (6.47). 
Under  the  notations  in  Remark  6.4,  the  result  is 

p 

QiH\  Qi  =  ^2  ck,iPk,o  =  ^2  ~akPk, 

k= 1  k 

where  we  have  used  (6.34)  and  (6.36),  and  summed  up  those  terms  with  the  same  ck,i 
to  make  the  ak  distinct.  In  case  =  0,  we  multiply  the  equation  by  Q0  from  the  left 
and  from  the  right.  Notice  that  the  right-hand  side  is  unchanged  by  (6.8).  Then  we  use 
the  same  notaton  Qx  to  denote  the  subprojection  QXQ0  =  Q0QX.  Thus  whatever  /3X  is, 
Qi  —  QiQo  —  QoQi-  By  (6.48),  the  above  equation  becomes 

J2<*kVk  =  -Q1Q0H1Q0Q1  =  Q1ASAQ1. 


(6.49) 
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The  ak  and  Vk  associated  with  (3\  are  so  uniquely  determined.  In  the  same  way  we  obtain 
all  the  ak  and  Vk  for  the  distinct  fa,  and  hence  the  total  number  m!  of  the  ak. 

In  the  special  case  that  all  the  Pk  are  simple,  we  have  m  distinct  fa  since  the  range  of 
Qo  is  m-dimensional.  Thus  m'  =  m,  and  the  Vk  are  uniquely  determined  by  (6.8).  Notice 
that  Qk  =  Vk.  Equation  (6.49)  becomes  a{P\  =  V\ASAV\ .  The  rank  of  V\  is  1.  Thus 
oq  =  =  tr (AS AVi ),  using  the  fact  that  tr  (ST)  =  tr(TS). 


6.2.  Applications. 

We  now  consider  two  important  examples  in  continuum  mechanics,  the  Navier- Stokes 
equations  and  the  equations  of  magnetohydrodynamics.  These  equations  will  be  studied 
in  more  details  in  Section  9.  Here  we  simply  apply  the  result  of  last  subsection  to  obtain 
the  Green’s  functions  of  their  linearizations. 

First  consider  the  linearization  of  the  Navier-Stokes  equations  (9.1)  around  a  constant 
state  (i>*,0,e*)  with  v* ,  e*  >  0: 

(6.50)  wt  +  Awx  =  Bwxx, 

where 

w  =  (v  —  v*,  u,  e  +  ~u2  —  e*)*, 

(6-51)  /0  -1  0\  /  0  0 

A  =  \  pv  0  pe],  ^=(  0  ^ 

Vo  p  o)  \^ev  o 

and  p  >  0,  0e  >  0,  k  >  0,  p  >  0,  pv,  pe  and  6V  are  constants.  In  consistency  of  (9.3)  and 
(9.5),  we  assume 

6 

(6.52)  Pv=Pv-  PPe  <0,  pv=Pv-  Yje  <  0,  Pe{0v  ~  pOe)  <  0- 
The  matrix  A  has  eigenvalues 


(6.53) 


A2  =  0,  A3  = 
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The  right  eigenvectors  r*  and  left  eigenvectors  lit  i  =  1,  2,  3,  for  A  are 

ri  =  (1, -Aj,  — p)*,  i  =  l,3, 

f2  =  (Pe,0,  -?„)*, 

^6'54^  it  =  ^-(p«,Aj,pe),  *  =  1,3, 

2 p„ 

i2  =  --=-(1,0,  -), 

Pv  P 

satisfying  liVj  =  6ij ,  j  —  1,  2,  3.  Assumption  6.1  is  implied  by  Assumption  2.1  for  the 
nonlinear  system,  which  will  be  varified  in  Section  9.  Assumption  6.2  is  satisfied  by  (6.51) 
and  (6.54).  Applying  Theorem  6.2  we  obtain  the  Green’s  function  G  of  (6.50),  (6.51), 

(6.55)  G(x,t)  =  yj-U===e“l  w’  +  0(l)(t  +  l)_*r  *  ^e_<*  +  H, 


i- 1 


where  pi  =  kBri ,  and  H  —  ^  e  Qfct^(x  —  j3kt)Vk<  Direct  computation  gives 

fc=i 


(6.56) 


Pi  =  7T— —(UPv+KPeiOv  -l>0e))  >0,  %  =  1,3, 


P2  = 


2v*pv 

KOefiv 

v*pv 


>0. 


To  compute  ify  we  consider  two  cases. 

Case  (i).  p  >  0.  In  this  case  ra'  =  m  =  1.  The  right  eigenvector  r0  and  left  eigenvector  l0 
of  B  associated  with  the  zero  eigenvalue  are 

r0  =  (1,0, -0v/6eY,  Jo  =  (1,0,0), 

satisfying  l0r0  —  1.  From  Theorem  6.2  and  Remark  6.4, 


{3  —  IqAtq  —  0,  V  —  ro/o5 

s  =  j(0, 1, 0)*(0, 1, 0)  +  ^(0, 0,  iy(6v/0e, 0, 1), 

a  =  tr  {ASAV)  =  l0ASAr0  =  pv. 

Thus  in  (6.55), 

/I  0  0\ 

(6.57)  H  =  e^p^S(x)  0  0  0. 

V-Sf  0  0 ) 
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Case  (ii).  fi  =  0.  In  this  case  m  =  2.  The  right  eigenvectors  r0j  and  left  eigenvectors  l0j, 
j  =  1,  2,  of  B  associated  with  the  zero  eigenvalue  are 


roi  =  (lA-0v/0e)\  r02  =  (0,1,0)*, 

loi  =  (1,0,0),  to  =  (0,1,0), 

satisfying  lojrok 

(6.58) 
i.e. 

(6.59)  Pi  =  -y/-pv ,  P2  =  V~~Pv 


=  bjki  j,  fc  =  1,2.  The  Pk  are  the  eigenvalues  of 

"o1)’ 


Both  /?i  and  ft2  are  simple.  Thus  m!  =  2.  We  have  from  (6.58) 


QqAQq  =  (r0i,r02)  ^  ?  q  ^  (/02)  =^1^1 

where 

We  also  have 


(6.60) 


Thus  in  (6.55), 


5=  -^(O,O,l)‘(0»/0e.O.l), 

Kl7g 

04  =  ti(ASAVk)  =  (^ioi  - 

=  ^(p-ev/ee)  =  a>0,  k  =  1,2. 


—ln)ASA(r  m-hrta) 
2pk 


,2  (  1  -1/A  0\ 

F=ie-QtY<S(x-At)  -A  1  0 

2  \-ev/ee  ev/(8A)  0  J 


with  a  and  A.  k  =  l,2,  given  by  (6.60)  and  (6.59). 

Next  consider  the  linearization  of  the  equations  of  magnetohydrodynamics  (9.13)  around 
a  constant  state  v  —  v*  >  0,  u\  —  u2  —  n3  =  0,  B2  =  -BJ,  B3  —  B%,  e  =  e*  >  0.  We  thus 
have  system  (6.50), 


wt  +  Awx  =  Bw  XX  5 
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where 


W  =  (v-  v*,  «1 ,  «2,  «3,  i;B2  -  v'B^vBa  -  u*BJ,  E-E'f, 
E*=e*  +  (B*2  +  B3*>  7(2/10), 


Cl 

C3B2 

A  =  C3B3 


0  0  0 

c2B|  c2S|  pe 
-c3  0  0 

0  c3  0 

0  0  0 


P  —C3B2V*  -c3B3v' 

0  0  0  0 

£  0  0  0 

0  £  0  0 

004  0 

r  0  0  0 

f  0  0  0  0 

0  0  0  <14.62 


1  ( 1 

C2  =  —  — 

MO  \V* 


P  =  p  +  (622  4-  BJ2)/(2p0), 


£fl.  +  (BS»+W=£ 


2p0  cr/ig; 


<24  =  — 


Mo  ^Mov’ 


with  p,  0C,  po  and  «  being  positive  constants,  */,  p  and  1/cr  nonnegative  constants,  and  pw, 
Me 5  -^1  and  constants.  Again  we  assume  (6.52)  for  those  constants  related  to  thermo¬ 
dynamics.  The  matrix  A  has  eigenvalues  in  nondecreasing  order  as 

(6.62)  cj ,  ca,  Cg ,  0,  Cg ,  ca,  Cf , 

where  c^,  ca,  cs  are  the  fast,  the  Alfven,  and  the  slow  wave  speeds,  respectively,  given  by 

.2  _  Bf 


cf,s  =  2  [  ±  \d4  +  4Pvb¥/(hov*)  , 

d2  =  —pv  +  ( B l2  +  Bf  +  Bf)j{^\ 
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with  pv  defined  in  (6.52).  In  the  following  cases  some  of  the  eigenvalues  coincide: 

Case  A.  When  B*  =  0,  cs  =  ca  =  0,  thus  A  =  0  is  an  eigenvalue  of  multiplicity  5. 

Case  B.  When  B^  ^  0  and  =  B%  —  0,  we  have  c2  =  ma x(-pv,cl),  and  c2s  = 
min (~pv,cl).  If  -pv  ^  cl,  either  cj  =  c2a  or  c2s  =  c2a,  hence  A  =  ±ca  are  eigenvalues 
of  multiplicity  2.  If  —  pv  =  cl,  the  multiplicity  of  A  —  ±ca  is  3. 

Set 

r(c)  =  (1,  -c,  cB'BJ/p, cB\BHp,  - B?B*2/p , -B?B*3/p, 

-p-  (B2*2  +  BfXpov’c2  +  B{2)/(2p0p))\ 

(6.64)  1(c)  =  (2  +  2BJ2(BJ2  +  B3*2)/p2)_1(l  -PPe/c2,-l/c,B1*B27M,B1*B37M, 

peB*2/(p0c2)  -  Bl/p,peBl/(p0c2)  -  Bl/p,  -pe/c2), 
p  =  p(c)  =  p0 v*c2  -  B[2. 

We  now  write  down  the  eigenvectors  of  A  satisfying  (6.2)  for  each  case. 

Case  A.  B{  =  0.  The  matrix  A  has  eigenvalues 

(6.65)  Ai  =  —Cf,  A2  =  0,  A3  =  Cf, 

with  mi  =  m3  =  1,  m2  —  5.  The  right  eigenvectors  are 

ni  =  r(-Cf),  r31=r(cf), 

r2 1  =  (1,0,0,0,0,0,-Ci/pe)', 

(6.66a)  r22  =  (0, 0, 1, 0, 0, 0, 0)(, 

r23  =  (0,0, 0,1, 0,0,0)*, 

r2. 4  =  (0,0,0,0,l,0,-c2B2/Pe)‘, 

r25  =  (0,0,0,0,0,l,-C2B37pe)t. 

The  left  eigenvectors  are 

hi  =  K~cf)i  hi  =  Kcf)i 

/  b*2  4-  B*2\  -1 

hi  =  (  ~~pv  "1 - - - )  (Ppei  0?  0?  c2^2  )  c^B3  ,Pe)i 

V  Pov*  ) 
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(6.66b)  l22  =  (  0,0, 1,0, 0,0,0), 

«23  =  (0,0, 0,1, 0,0,0), 

Jm  =  (0,0, 0,0, 1,0,0), 
hs  =  (0,0, 0,0, 0,1,0). 

Case  Bl.  B1  ^  0,  B2  =  S3  =  0,  and  —pv  5^  c^.  The  matrix  A  has  eigenvalues 

(6.67)  Ai  =  —\/—pv,  A2  =  —  ca,  A3  =0,  A4  =  ca,  A5  =  C~Pv, 

with  mi  =  m3  =  ms  =  1,  m2  =  =  2.  Here  Ai  and  A5  are  the  fast  waves  or  the  slow 

waves.  The  right  eigenvectors  are 

rn  =  (i>  V—pH,  0, 0, 0, 0,  ~pY, 
r51  =  (l,--/^,  0,0, 0,0,  -p)\ 
r3i  =  (1,0, 0,0, 0,0,  -Pv/Pe)1, 

(6-68«)  »-2i  =  (0,0,8ign(BI)/V5Sr,0, 1,0,0)*, 

r22  =  (0,0,0,sign(Si)/v^v*,0,l,0)‘, 
r4i  =  (0, 0,  -  sign(S1*)/v/^ow’,  0,1,0, 0)*, 
r42  =  (0, 0, 0,-sign(Si*)/vW5',  0,1,0)*. 

The  left  eigenvectors  are 

*11  =  j  (Pv/Pv,  l/\/::p((,0,0,0,0,pe/pv), 
hi  =  ^(Pv/Pv,-'i-/\/z:pO,0,0,0,0,pe/pv), 
kl  =  (  PPe/Pv>  0,  0,  0,  0,  0,  Pe/Pv)i 
(6-686)  hi  =  ^  (0, 0,  y/ftov*  sign(Sj ) ,  0, 1 , 0, 0) , 

*22  =  i(0,0,0,Vw^sign(S1‘), 0,1,0), 

Ui  =  ^(0, 0,  -i/ipyu*  sign(SJ),  0, 1,0, 0), 

*42  =  I  (0, 0, 0,  -  sfpav*  sign(Sj( ) ,  0, 1 , 0) . 
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Case  B2.  Bjf  ±  0,  B|  =  =  °>  and  =  ca-  The  matrix  A  has  eigenvalues 


(6.69) 


A,  —  Cat  A2  —  0,  A3  Ca, 


with  mi  =  m3  =  3,  m2  =  1.  The  eigenvectors  are  the  same  as  in  Case  Bl. 

Case  C.  B[  ±  0,  and  Bf-  +  Bf  ±  0.  The  matrix  A  has  eigenvalues  A,  given  by  (6.62), 
with  rm  =  1,  i  =  1,  . . .  ,7.  The  right  Eigenvectors  are 

n  =  r(-Cf),  r7  =  r(cf),  r3  =  r(-c5),  r5  =  r(cs), 


r4=  1,0,0,0,55,55,-^  + 


pv  B;2+B?2n‘ 


(6.70a) 


r2=  0,0 


Pe  2p0 

BJ  sign(Bj)  Bj  sign(B 


vw 


—,BZ,  — B5,0  )  , 


„  BJ  sign(Bi)  BJsign(BJ) 
r6=(,0’0’-  27 
The  left  eigenvectors  are 

h  ~  K~ c/)>  ^7  =  Kcf)i  ^3  — 


„)• 


(6.706) 


l  -  Pe 
U  —  —  — 

Pv 

k  = 


(p, 0,0,0,-  —  , 

V  Mo 


Mo 


,1  K 


*2  ' 


*6  = 


2 (B?  +  Bs  j 
1 

2(BJ2  +  BJ2) 


(0, 0,  sign(BJ),  -y/ntiv*BZ  sign(BJ),  BJ,  -BJ ,  0), 

(0,0,  ->/Mov*BJ  sign(BI),  y/nov'Bl  sign(B*),  BJ,  -BJ,0). 


Assumption  6.1  is  implied  by  Assumption  2.1  for  the  nonlinear  system,  which  will  be 
verified  in  Section  9.  Notice  that  we  have  assumed  k  >  0.  By  the  above  expressions  for 
the  eigenvectors  of  A,  we  see  that  Assumption  6.2  is  satisfied  for  the  following  choices  of 
parameters: 

(i)  v  >  0,  [i  >  0,  and  1/cr  >  0; 

(ii)  v  =  [i  =  0,  and  1/cr  >  0,  while  B{  ^  0; 

(iii)  v  >  0,  fi  >  0,  and  1/cr  =  0,  while  B{  ^  0. 

Now  we  apply  Theorem  6.2  to  (6.50),  (6.61).  The  Green’s  function  of  the  linearized 

equations  of  magnetohydrodynamics  is 

G(x,t)=G*(r,f)  +  0(l)(l  +  ir5rTe"!i^L  +H - 

3 


(6.71) 
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where  G*  is  the  Green’s  function  of  a  diagonalizable  parabolic  system 
(6*72)  wt  +  Awx  =  B*wxx , 

the  A j  are  the  eigenvalues  of  A  given  by  (6.62),  (6.63),  and  H  =  £  e~akt6(x  -  /3kt)Vk . 

k=i 

The  matrix  B*  is  given  by  (6.5).  Notice  that  we  have  computed  all  the  eigenvectors  of  A. 
Hence  the  computation  of  G*  is  straightforward.  For  example,  in  Case  C  we  have 

(6-73) 


where 


Mi  =  Mr  =  M(c/),  M3  =  Ms  =  m(c«). 

M(c)  =  l(c)Br(c)  =  [2  +  2  B?(B?  +  Bfj/p2]-1 


(6.74)  '  y*c2Pe^eP  6v)  +  v*  +  v'p2  +  B^  +  +  ’ 

M4  =  UBr,  =  ^A, 

V*Pv 

M2  =  M6  =  I2BT2  =  leBr6  =  ~  +  - — I — - 

2^;*  2anQv*2 

The  other  cases  are  somewhat  tedious.  We  now  concentrate  on  obtaining  H.  Denote  the 
eigenvalues  of  B  as  z/ji 


(6.75)  Ui  =  0,  v2  =  ~,  j/3  =  i/4  =  ii  i/b  =  u6  =  — 1~— 

v*  v*  a  ix  0v*2 

The  corresponding  right  eigenvectors  of  B  are 


Pi  =  (1, 0, 0,  o,  £?2 ,  b3  ,  (bj2  +  b3*2)/(2Mo)  -  ov/ee)\ 

P2  =  (0)l,0,0,0,0,0)(, 
p3  =  (0,0, 1,0, 0,0,0)*, 
p4  =  (0,0,0, 1,0, 0,0)*, 
p5  =  (0,0,0,0,l,0,^//io)', 

p6  =  (o,o,o,o,oIi,B;/#i0)t, 

p7  =  (0,0,0, 0,0, 0,1)*, 


(6.76a) 
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and  the  left  eigenvectors  are 


qi  =  (1,0, 0,0, 0,0,0), 

q2  =  (0,1, 0,0, 0,0,0), 
q3  =  (0,0, 1,0, 0,0,0), 

(6.76b)  q4  =  (0,0,0, 1, 0,0,0), 

q5  =  (-B2*,0,0,0,l,0,0), 
q6  =  (-53*,0,0,0,0,l,0), 

q7  =  ((B*2  +  B*2)/(2/i0)  +  ev/ec) 0,0,0,  -B*2/ti0, -B37^0,l), 


satisfying 


c\ipj  —  6ij,  i,j  —  1, . . . ,  7. 


Case  (i).  v  >  0,  fi  >  0,  and  1/a  >  0.  In  this  case  m!  =  m  =  1  by  Theorem  6.2,  and 


V  =  piqu  =  0. 


By  Remark  6.4, 


&  ~  X]  , ..  M 

i=2  ^ 

a  =  tr(ASAV)  —  qi  AS Api  =  -v*pv/v, 


where  pt,  is  defined  in  (6.52).  Hence 


(6.77) 


H  =  ev'Mvb(x)p  iqi. 


4 

Case  (ii).  i/  =  p  =  0,  1/a  >  0,  and  B\  ^  0.  In  this  case  m  =  4,  and  Qo  =  E  PjqT  Direct 

j=i 


computation  shows  that 


(6.78) 


/° 

=  (P1,P2,P3,P4)  Pq 
\  0 


-1 

0 

°\ , 

/ql 

0 

0 

0 

q2 

0 

0 

0  1 

q3 

0 

0 

0/ 

V  q4 

=  /?ipiqi  +  02p2h, 
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where 

A, 2  =  Ty/'-Pv, 

(6'79)  fi  „  „  1  1 

Pi, 2  =  Pi  -  Pl,2p2,  qi,2  =  xqi  -  ^5 - <\2- 

1  1,2 

It  is  easy  to  see  that  the  right-hand  side  of  (6.78)  is  the  spectral  representation.  Hence 
with  the  notations  in  Remark  6.4, 


(6.80) 


Notice  that 


—  Qk  —  P/c^/e?  k  —  1,2, 

03  —  0,  Qs  =  Qo  Qi  —  0,2  —  p3^3  +  p4^4- 


5  =  E~Mr 


-  Vn 

j= 5  3 


We  then  have  for  k  =  1,  2, 


(6.81) 


Qfc  =  qkASApk  =  ~(B*22  +  Bf)  +  = 


a  >  0 


by  (6.52).  By  (6.80)  and  direct  computation, 


Therefore, 


In  conclusion, 


Q3AS^4Q3  —  — — ^-(p3q3  +  p4q4)  =  av*Bl2Q2. 

^5 


a3  =  av*B{2,  V3  =  Q3. 


(6.82)  H  =  e  °‘t[5{x  +  ty/^pv)plql+8{x-ts/^fv)p2(\2\+e  <TV‘Bi2t6(x)(p3q3  +  p4q4); 

where  a  and  pfc,  qfc,  k  =  1,2,  are  defined  in  (6.81)  and  (6.79). 

Case  (iii).  v  >  0,  p  >  0,  1/<t  =  0  and  B[  ±  0.  In  this  case  m  =  3,  and 

<2o  =  Pi9i  +  p595  +  p696- 

Direct  computation  shows  that  Q0AQ0  =  0.  Hence  f3  =  0  and  Q  =  Q0.  Notice  that 

S  =  —p2q2  +  — (p393  +  p4fl4)  H - p7<l7- 

^2  V3  V7 
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We  now  compute  otk  ■  We  have 
(6.83) 


QASAQ  =  (pi,p5,p6) 


_  BZ  Bi  \  1 

Pv  >  - *  j  - *  J  H - 

VoV*  fMQV *  )  is3 


(0 

\ 


flQV* 


If  Bf  +  Bf  =  0,  then 

QASAQ  =  -^-piqi  +  —  (p5qs  +  Peqe)' 
v 

The  right-hand  side  is  the  spectral  representation  ol\V\  +  a2B2.  Therefore, 

i>vv*t  B*2t 

(6.84)  H  =  8(x)[e  »  piqi  +  e~~^  (p5qs  +  Peqe)] . 

If  BJ2  +  B32  ^  0,  (6.83)  has  the  spectral  representation 


QASAQ  —  a\P\  +  a2B2  +  0:3^3, 


where 

(6.85) 


*2 


Oil  = 


VOV 


“2'3  =  2 


u>2  ±  +  ApvB^v* /{pav p) 


~f~  «1, 


2  pvv *  ,  B?  +  B?  ,  Bf 

UJ  — - 1 - 1 - , 

V  vvo  VVo 

V\  =  [B^2p5(\b  +  ^22p6q6  -  ^2  ^3  (p6^5  +  p5^6)]/ (^22  +  Bg2), 
.2  T-l 


P2,Z  = 

Therefore, 

(6.86) 


(a2)3-«i)2-^^(B2*: 

VO  V*Pv 


Bl 


[(<*2,3  “  <*l)pl  "  ^2,3^2P5 


,  ,  .  Ct2,3&2  ,  a2,3^3 

(a2,3  -  cti )qi  +  -  qs  +  <|6 


Vo  v  Pv 


Vo  V  Pv 


H  =  6{x)  [e~aitVi  +  e“Q2'B2  + 


OC2,zBIP§] 


where  a*  and  Vi,  i  —  1,  2,  3,  are  given  by  (6.85). 


7.  Energy  Estimate 


In  this  section  we  present  the  energy  estimate,  which  is  due  to  Kawashima,  for  the 
Cauchy  Problem  (2.1),  (2.2)  of  viscous  conservation  laws.  The  reader  is  referred  to  [Ka] 
and  the  references  therein. 

Theorem  7.1.  Suppose  that  Assumptions  2. 1-2.3  hold.  Let  m  >  2  be  an  integer  and 
assume  that  uo  £  Hm.  If  ||uo||//m  is  small,  then  the  initial  value  problem  (2.1),  (2.2)  has 
a  unique  global  solution  u(x ,  t)  satisfying 

um  £  C°([0, 00);  Hm)  n  C'ao,  oo); 

(7.1)  u(2)  6C0([0,oo);Fm)nC'1([0>oo);ffm-2), 

€  L2([ 0,  oo);  H m“1),  €  L2([ 0,  oo);  Hm), 

where  uw(x,t)  and  u^2\x,t)  are  the  orthogonal  projections  of  u{x,t)  =  on 

M  and  AfL  respectively;  f0  is  the  mapping  in  Assumption  2.2;  and  A f  is  the  null  space  of 
£?(/o(u))/o(u).  The  solution  satisfies 

(7-2)  sup  +  /  (ll«x(-,r)||2fm_1  +  ||42)(-,r)||2fm)  dr  <  C\\u0\\2Hm 

if  Q 

for  all  t  >0,  where  C  is  a  constant. 

Before  proving  Theorem  7.1,  we  restate  (2.1),  (2.2)  and  Assumptions  2.1-2.3  here  for 
convenience.  We  consider  the  Cauchy  problem 

(7-3)  ut  +  f(u)x  =  ( B(u)ux)x , 

(7*4)  u(a:,0)  =  u0(x), 

Here  (7.3)  satisfies  the  following  conditions: 

(i)  There  exist  an  entropy  pair  U(u)  and  F(u).  Set 

(7-5)  A0(u)  =  V2U(u),  A(u)  =  f(u). 

Then  (VU) A  =  VF,  Ao  is  positive  definite,  AqA  is  symmetric,  and  A0B  is  symmetric 
semi-positive  definite. 
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(ii)  For  small  u,  B(u)  /  0.  There  exists  a  smooth  one-to-one  mapping 

(7.6)  u  =  /o(«),  /o(0)  =  0, 

such  that  the  null  space  A f  of 

(7.7)  B(u)  =  B(f0(u))f0(u) 

is  independent  of  u.  Moreover,  A fL  is  invariant  under  fb(u)tAo(fo(u )),  and  B(u )  maps  R 
into  A/'1. 

(iii)  Any  eigenvector  of  A(0)  is  not  in  the  null  space  of  B( 0).  By  Lemma  6.6,  this  implies 
that  there  exists  a  constant  matrix  K  €  RnXn,  such  that  KAo(0)  is  skew  symmetric  and 
{~[KA0A  +  (KAqA)*]  +  AoB}(0)  is  positive  definite. 

Assumption  (ii)  allows  us  to  consider  such  a  system  (7.3)  whose  viscosity  matrix  B(u) 
has  a  null  space  depending  on  u,  but  after  transform  (7.6)  the  system  for  u, 

fo(u)t  +  /(/o(^))x  =  {B(u)ux)x, 

has  a  new  viscosity  matrix  B(u)  whose  null  space  is  independent  of  u.  See  examples  in 
Section  9. 

Since  the  local  existence  is  standard  [K],  to  prove  Theorem  7.1  it  is  sufficient  to  obtain 
the  following  a  priori  estimate: 

Lemma  7.2.  Assume  conditions  (i)-(iii).  Let  m  >  2  be  an  integer  and  let  T  >  0  be  a 
constant.  Let  u(x,t )  be  a  solution  to  (7.3),  (7.4)  satisfying  (7.1)  with  [0,  oo)  replaced  by 
[0, T\.  Set 

(7.8)  Nm(t)2  =  sup  \\u(-,T)\\2Hm  +  (  (\\ux{',T)\\2Hm-i  +  \\u^\',r)W2Hm)  dr. 

o <T<t  Jo 

If  Nm{T)  is  bounded  by  a  small  positive  constant  independent  ofT,  then  Nm(T)  <  C\\uq\\h™  , 
where  C  is  a  constant  independent  of  T. 

Proof.  To  simplify  the  notation  we  denote  ||  *  \\l?  by  ||  ■  ||.  Without  loss  of  generality  we 
assume  that  L/ (0)  =  F( 0)  =  0  and  /( 0)  —  0.  Set 


£{u)  =  U(u)  -  V£/(0)u. 


QUASILINEAR  HYPERBOLIC-PARABOLIC  SYSTEMS 


89 


Clearly  S  is  equivalent  to  \u\2  for  small  u  since  A0  =  V2U  is  positive  definite.  By  condition 
(i)  and  (7.3), 


S(u)t  +  (F(u)  -  VU(0)f)x  =  (VU(u)  -  VU(0))(B(u)ux)x. 

Integrate  this  equality  over  R  x  [0,  t]  for  0  <  t  <  T.  We  have 

pt  pOO 

(7.9)  c1\\u\\2(t)+  /  (VU(u))xB(u)uxdxdT<c2\\u0\\2, 

J  0  J — oo 

where  Ci  and  c2  are  positive  constants.  By  conditions  (i)  and  (ii)  the  integral  in  (7.9)  is 

pt  p  oo 

/  /  vJxAq(u)B(u)ux  dxdr 

JO  J — oo 

=  [  [  >  c3  f  ||«^||2(r) dr, 

JO  J- oo  Jo 

where  c3  >  0  is  a  constant  and  u ^  is  the  projection  of  u  on  Af1 .  Thus  (7.9)  implies  that 

(7-10)  IM|2W  +  f  \\u^W2{r)  dr  <C\\uo\\2 

Jo 

for  t  e  [0,  T],  where  C  is  a  constant  independent  of  T.  Here  and  later,  Nm(T)  is  assumed 
to  be  bounded  by  a  small  positive  constant  independent  of  T. 

We  then  estimate  the  L 2  norm  of  derivatives.  Applying  /dx^ ,  1  <  j  <  m,  to  (7.3), 

multiplying  it  by  A0(w)  from  the  left,  and  taking  the  inner  product  (■,  •)  with  dju/dxj ,  we 


1  /dju 

2  \  dxi 


.  dju\  /  dju  dj  \ 

’Mu)d^/t  +  \d^’Mu)d^{A{u)Ux)) 

/dju  ,  ,  &  .  ,  \  1  /&u  A  ,  . 

=  \m’Mu)mWu)v*)x)  +  2\d^’Mu)t 


&u\ 
dxi  /  ’ 


Integrate  this  equation  over  R  x  [0,  t]  for  0  <  t  <  T.  By  integration  by  parts  and  condition 

(i)  we  have 

(7.11) 

cilt5|  {t)+l  LS^,Mu)^{B{u)ux))dxdT 

<  cM3)f-Jo  I_oo(^Mu)[^j(A(u)ux)  -A(u)^tydxdT 

+  \L  L(S’iMu)A{u))*]S)dxdT 


-I  JS^'Mu)x^{B{u)Ux))dxdT+\J  Jl(%S’Mu)tS)dxdT' 
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where  C\  and  C2  are  again  positive  constants.  Let  P  be  the  projection  operator  on  A/"-1", 
which  is  a  constant  matrix.  By  conditions  (i)  and  (ii)  the  second  term  on  the  left-hand 
side  of  (7.11)  is 


r  f°  (£y(/o(*)fix),  A0(u)£-(PB(u)u^  dxdr 

-ii 


t  poo  /  fsj  \ 

Vr,  !'MMu)Pi-{B{u)u^) )  dxdr  + 


'CO 

=  jT  f°  dxdr  +  0(l)Nm(t) 

>-k:k x 


(r)dr-CNm(t)3, 

where  C3  >  0  is  a  constant,  and  we  have  noticed  m  >  2,  1  <  j  <  m.  Thus  (7.11)  is  reduced 
to 


(7.12) 


dJu 


dxi 


\t)+  / 

Jo 


t  m  Qjum  ..2 


dxi 


(T)dr<C\\uo\\2Hm+CNm(t): 


1  <j<m, 


for  t  €  [0,  T]. 

It  is  sufficient  to  obtain  an  estimate  for  the  L2([0,  T];  Hm~l)  norm  of  ux.  Rewrite  (7.3) 
as 


(7.13) 
where 

(7.14) 


ut  +  A(0)i*x  =  h, 


h  =  (A(0)  -  A(u))ux  +  (B(u)ui2%. 


From  (7.13)  we  have  for  0  <  j  <  m  —  1,  0  <  t  <  T, 

(7.15) 

&>+1u\  .  .  rl  r°°  /&+1 


•t  roo 


a 


/di+1u  f 

\dxi+1'KA°^>  dxidt 


dxdr  +  J*  f  KA0(0)A(0) 


d^lu 
dx i+i 


dxdr 


-Hl(^KMo)^)dxdT' 


where  K  is  as  in  condition  (iii).  For  the  first  term  on  the  left,  we  integrate  by  parts  with 
respect  to  x  and  then  with  respect  to  £,  which  gives 


pt  r  oo 

JO  J — 0< 


J\^'KMQ)l^dxdT  = 


=  -Jl(^’KM0)^)ix’t)dX 


OO 
t  roo 


+ 


/_  {^\KAo(0)ut1])dx  +  Jo  J^^,KM0)d^)dxdT. 
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Notice  that  KA0( 0)  is  skew  symmetric.  Therefore, 


dj+1u  y  di+1u\ 

dxi+i ’  KA° ^  dJdt )  dxdT 


+  l £  (u^^MoWJ+v)  dx 


-  “  ^IMI  W  “  C\\Uo\\2Hm 

since  0  <  j  <  m  -  1.  For  the  second  term  on  the  left-hand  side  of  (7.15),  we  decompose 
KA0(0)A(0)  into  symmetric  and  skew  symmetric  parts.  Clearly  the  contribution  of  the 
skew  symmetric  part  is  zero.  Notice  that  {\[KA0A  +  (KAqAY]  +  Ao5}(0)  is  positive 
definite  by  condition  (iii).  Thus 

L  LS^,KMo)m^)dxdT 

=  1  j_J^^A\^AoA  +  {KA0Af)+A0Bmd^yxdr 


dj+1u  dj+1u\ 

■qTTTT  ,  A)(0)ff(0)  s )  dxdr 


f  dj+lu 

l0  dx3  +  1 


1  f1 

(r)  dr  —  C  / 

Jo 


dxi+1  / 
&+1uW 
~dx^ 


( r)dr-CNm(t )3, 


where  c  >  0  is  a  constant,  and  we  have  used  condition  (ii)  in  the  last  step.  From  (7.14)  we 
have  for  0  <  j  <  m  —  1, 

M  ^(M(“)-^(0)K)  «+  (B(u)uW)  ( t ) 


Hence  the  right-hand  side  of  (7.15)  is  bounded  by 


(718)  i|s«?||||Sl(r)*,scj( 

a  am 


dr  +  CNm(t\ 


dT+hJo  +CNm(t)\ 


where  e  >  0  is  an  arbitrary  constant.  Choose  a  small  e.  Then  (7.15)-(7.18)  imply  that  for 
0<t<T,0<j  <  m  —  1, 

ft\\dj+1u  2  f  rf  } 

(7'19)  Jo  1^+1  (r)*-C'|||u|&..(t)+j(  ll«i2)li^(r)drj  <C|K||^+CiVm(t)3. 
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Sum  up  (7.19)  for  0  <  j  <  m  -  1  and  multiply  the  result  by  a  small  constant  e  >  0. 
Sum  up  (7.12)  for  1  <  j  <  m.  The  addition  of  (7.10)  and  these  results  gives 

IMItf-M  +  +  lka)llfr»)W*  <  C||uo||2w™  +CNm(tf,  0  <t<T. 

Jo 

Taking  the  sup-norm  fot  t  we  have  Nm(T)2  <  CUttollnm  +  CNm(T)3,  which  implies 
Nm(T)  <  C||«o||//">  for  small  Nm{T).  1=1 


8.  Systems  with  Non-Di agon aliz able  Linearizations 


We  now  prove  the  main  result  of  this  paper,  Theorem  2.3.  Follow  the  notations  in 
section  2.  Decompose  the  solution  u  to  (2.1),  (2.2)  in  the  right  eigenvector  directions 
of  /'( 0).  With  A i,  Ti  and  k  being  eigenvalues  and  matrices  consisting  of  right  and  left 
eigenvectors  of  /'( 0),  given  by  (2.5),  (2.6),  we  have 

s 

(8.1)  u(x,t)  =  '}TriUi(x,t),  Ui(x,t) —  liu(x,t),  z  =  l,...,s. 

i=l 

Rewrite  (2.1)  in  terms  of  the  U{ : 


(8.2)  Ufa  -f-  XiUix  —  TjUjxx  “t“  9ix  “1“  hixi  V  —  1,  .  .  .  ,  5. 

J= 1 

where 


(8.3)  gi  =  gi(u)  =  A iUi  -  lj(u ),  hi  =  ht{u)  =  k(B(u)  -  B(0))ux. 

Let  the  Green’s  function  associated  with  the  linearization  of  (8.2)  be  G(x,  t).  By  Duhamel’s 
principle, 

/OO  rt  pOO 

G{x-y,t)u(y,0)dy+  /  /  G(x-y,t-T)[g{y,T)  +  h{y,T)]vdydT, 

- oc  Jo  J  —  oo 

where 


The  equations  for  the  diffusion  waves  are  given  by  (2.11): 


(8.6)  eit  +  \ieix  +  -lif"(0)(riei,riei)x  =  liB(0)rieixx,  t  = 

Let  the  Green’s  function  of  the  linearization  of  (8.6)  be  G*(x,t).  We  have 

/OO  pt  poo 

G*(x  —  y,  t)0(y,0)  dy  +  /  /  G*{x  -  y,t  -  T)g*{y,r)  dydr, 

■OO  Jo  J  —  OO 


93 


94 

where 

(8.8) 

Set 

(8.9) 


TAI-PING  LIU  AND  YANNI  ZENG 


e  = 


e , 


(A 


9i=-7;hf"mri9i,ri0i)i  t  =  l,..., 


v  =  u  -  9  =  (ux  -  Oi, . . . ,us  -  0a)f  =  (vi, .  •  -  ^s)*- 


(8.10) 


+ 


+ 


Equations  (8.4)  and  (8.7)  imply  that 

/oo  poo 

G”{x-y,t)v(y,0)dy+  /  (G  -  G*)(x  -  y,  t)u(y, 0)  dy 

-oo  J  —  oo 

pt  poo 

/  /  G*(x-y,t-T){g- g* +  h)y(y,T)dydr 

[  [  (G  —  G*)(x  —  y,t  —  r)(g  +  h)y{y,r)  dydr. 

J  0  J —oo 

As  in  Section  4,  we  may  assume  by  Lemma  2.1  that 

/?B(0)ri  =  diag(/iii, . . .  ,/iimi),  i  =  1,...  ,s, 
with  fiij  >  0,  j  =  1,  ...  ,  raj.  Thus  the  componentwise  version  of  (8.10)  takes  the  form 

/°°  I  (x-y-Xjt)2  fOO 

. rj - -e  ^  vi:j(y,  0)  dy+  /  (G  -  G*)ij(x  -  y,t)u(y,0)dy 

-oo  J —oo 

^  r*  r°° _ i_ 

Jo  J  —  oo  -  r) 


(a;-i/-A£(t-r)) 

=e  (t_T)  (fly  -  9,*,-  +  hij)y{y, r)  (tydr 


/•r  poo 

+  {G-G*)ij(x  -y,t-r){g  +  h)y(y,T)dydT, 

Jo  J — oo 

J  —  15  •  '  •  7  ^  *  *  *  5  ^5 


where  Ujj,  pij,  g*j  and  hij  are,  respectively,  the  jth  components  of  Vi,  gi,  g*  and  hi, 
(G  -  G*)ij  is  the  corresponding  row  of  G  —  G* .  We  have  similar  formula  for  higher 


QUASILINEAR  HYPERBOLIC-PARABOLIC  SYSTEMS 


95 


derivatives.  In  general,  for  Z  =  0,  1,  . 
dl 


(x-y-Xjt)2  Ql 

e  44”‘  ^nvij(y<°)dy 


(8.11) 


Qyl 

f°°  r)1 

+  J_  (G -G*)ij(x-y,t)—lu(y,0)dy 

+U 


rt  poo 

+  /  /  (G  -  G*)ij(x  -  y,t  - 

J  0  J  —  oo 


dyl 

(X-T /-Xj(t-T))2  ^+1 

-oo  B,J<  ’  ^i+r (9ij-9tj+hij)(y,T)dydi 

dl+1 


dyl+1 

j  =  z  = 


(S  +  ^)(y,  T)dydr, 


Applying  Theorem  6.15  to  the  linearizations  of  (8.2)  and  (8.6),  we  obtain  for  l  =  0,  1,  ...  , 


dl 


(8-12)  g^(G  -  G*)(x,i)  =  0(l)(t  +  1  J-iH*1  £e- 


(x-Afct)J 

Ct 


*=1 


m'  f 

fc=l  fc'=0 


where  m'  is  a  nonnegative  integer  not  more  than  the  dimension  of  the  null  space  of  13(0), 

ak  >  0  and  @k  are  constants,  and  Qkk,  are  polynomial  matrices  in  t  with  degrees  not  more 
than  k',  k  =  l,  ...  ,m',  k'  =  0,  ,  l. 

We  now  perform  a  priori  estimate  through  a  sequence  of  lemmas.  Adopt  the  notations 
in  Section  2, 

lM*>  f;  A)  =  \(x  -  A (t  +  l))2  +  t  +  l]~a/2, 

(8.13)  t;  A)  =  [|x  -  A(t  +  1)|  +  t  +  l]-“/2, 

S 

*)  =  y>3/2(*,  t;  Xi)  +  (j>i (x,  t ;  Ai)V>3/4(a:,  A*)]. 

i=l 

Set 


(8.14)  M(t)=oSup(mgj|||ui(-,r)$(.,r)-1||L„  +  ||,;lx(.,r)(r  +  l)^(.,r)-i|| 


+ 


2  / 

E 

i=i  v 


<92^,- 


&r2* 


(■^) 


L°° 


<92i+1v?- 


<9x2*+1 


(■^) 


j(r  +  l)32'|. 
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Then  for  — oo  <  x  <  oo,  t  >  0,  i  =  1,  . . .  ,  s, 

|t>i(a:,t)|  <  M(t)4>(x,f),  \vix(x,t)\  <M{t)(t  +  l)'t*{x,t), 

<  M(t)(t  +  I)-1, 

<M(t)(t  +  1)"L 


(8.15) 

(Pvt.  J 
te* (Xlt) 

j 

d3v, .  . 

94V{  . 

5 

d^Vi .  , 

1  =  0,  1, 
(8.16) 


r°° _ i_ 

J-oo  \/4tt/ 


(g-y-A,t)2  Ql 


, _ 'T^i~lvij(y,0)dy  =  0(1)6*  (t  +  l)-^(x,t), 

s/topPji  dyl 

where  8*  is  the  norm  of  the  initial  data  defined  in  (2.24)- 

Proof.  By  (8.9),  (8.1),  (2.2),  (2.12),  (2.20)-(2.24)  and  (3.1),  we  have  for  1  <  i  <  s, 

Vi(x,0)  =  MiOr,0)  -  6i(x, 0)  =  kuo(x)  -  6,(x, 0),  ||vi(-,0)||w11  =  0(1)5*, 

/OO 

^«i(x,0)dx  =  0, 

IMU»  =  0(1)5*,  M*)|  =  0(1)5*(1  +  |*|)-1, 

where  x 

rn{x)=  /  Vi(y,0  )dy. 

J  —  OO 

For  l  =  0,  denote  the  left-hand  side  of  (8.16)  as  I.  For  \x  —  A;£|  <  (t  +  l)1^2, 

/oo  -I  (x-y-Ait)2  f00  „  (x-y-Xjt)2 

—===e  ri'ij(y)dy  =  0(  1)  /  t  le  4(#i«+;,t  \riij(y)\dy 

-oo  y/^TTVijt  J-oo 

=  0(1)6* (t  +  l)-1  -  0(l)6^3/2(*,  t]  Ai), 

where  is  the  jth  component  of  77* ,  and  e  a  small  positive  constant.  We  have  assumed 
t  >  1  in  the  above  estimate.  The  case  i  <  1  is  clear.  For  (t  +  l)1/2  <  \x  —  A*t|  <  t  +  1,  say, 
x  -  \{t  >  0,  then 

/OO  _ (x  — y  — A^t)2 

t-1e  |»jy(y)|dy 

-OO 

*(x-Ait)/2  ^  _(a:-Ait)5 


(8.18) 


=  0(1)  /  t  *e  C‘  |%-(l/)|dy 

J  —CO 

+  0(1)  [°°  t-1e-(*~*^t‘>a'5*(l  +  y)-1dy 

=  0(l)5*r1e-i£z^i2i  +  0(l)5*ri(l  +  1*  -  M)“l 
=  0(l)5*(t  +  l)_2|x  -  Aii|— 1  =  0(l)5*^/a(a:,<iAi). 
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The  case  that  x  -  Xrf  <  0  is  similar.  For  \x  -  A<*|  >t  + 1,  again  consider  x-\it>  0.  Then 


mj_ 


(x-\it)/2 


mL 


t  \e  cl*  \vij{y,0)\dy 


1  (as  —  y  —  A.i)2  c 

r^e - + 


-/(*-A4t)/2 

=  0(l)^t-2e-(  ct°  +  0(l)tf*|a:-  A*t|~5 

-  0(l)**|s  -  \it\~i  =  0(l)tf*^1(arl  t;  A,)^/^, i;  A*). 

Similarly  the  estimate  holds  for  x  -  \t  <  0.  We  thus  have  proved  (8.16)  for  l  =  0.  For 
/  =  1,  the  case  that  t  >  1  is  clear  by  integrating  the  left-hand  side  by  parts;  while  the  case 
that  t  <  1  is  similar  to  the  case  when  /  =  0.  n 


Lemma  8.2.  Under  the  assumptions  of  Theorem  2.3,  we  have  for  -oo  <  x  <  oo;  t  >  0, 

/°°  1  (x-y-Xjt)2  Ql 


OO  v4  KIHjt 

where  2  <  l  <  5  is  an  integer. 


Proof.  Denote  the  left-hand  side  of  (8.20)  as  I.  For  t  >  1,  we  perform  integration  by  parts, 


-°mj_ 


l  d~  1  (a;  Ajt)  ill 

t  *e  o,  |»y(y)0)|dy  =  O(l)(l+l)-T-||Vy(.,0)||z,i. 


Equation  (8.20)  follows  by  (8.17).  For  t  <  1 
/  =  0(1)  ^(-'0)  =0(1) 


1  =  ioo  =  0(1)IM,>°)llff<+>  =  0(1)<5‘(*  +  1)— = 


by  (2.24). 


Lemma  8.3.  Under  the  assumptions  of  Theorem  2.3,  we  have  for  — oo  <  x  <  oo,  t  >  0, 

1  =  0,1, 

/OO 

(G  ~  G*)ij(x  -  y,t)—lu(y,Q)dy  =  0(1)6*  (t  + 

Proof.  By  (8.12),  (8.5),  (8.1)  and  (2.24), 

(8.22) 

/OO 

(G  -  G*)ij(x  -  y,  t)u(y,  0)  dy 

■OO 

/OO  5  2  m' 

0(l)(t  +  l)~*t~2  ]Pe  cr*  |u0(2/)|  dy  +  y^O{l)e~t/c\uQ{x  -  f3kt)\ 

00  *=i  k=  i 

S -\  POO  tx_y_x  t\2 

=  0(1)  /  (t  +  l)“irie-—  toh*  \u0(y)\dy  +  0(l)<5*e-f/c(l  +  \x\)~i. 

t.  1  J  —  oo 
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Set 

/OO  ( 2C  "i  ^  ^2 

( t  +  l)~h~^e~  ~  |Uo(y)|  dy. 

-OO 

By  (2.24),  for  \x  -  \kt\  <  (t  +  l)i, 

4  =  0(l),5*(f  +  l)-1  =  0(l)6*(t  +  l)-h3/2(x,t-,  AO- 
Similar  to  (8.18),  for  (t  +  l)1/2  <\x-  Xkt\  <  t  +  1,  we  have  by  (2.24) 

Ik  =  0(l)5‘(t  +  1  +  0(l)fi*(f  +  1)^(1  +  \x  -  Afcf|)_‘ 

=  0(l)6’(t  +  l)-*1p3/2(x,t',  Afc). 

And  similar  to  (8.19),  for  \x  —  A^t|  >  t  +  1,  we  have 

Ik  =  0(l)S*(t  +  \k)i>3/4(xi  t]  Afc). 

Together  with  (8.22),  (8.23)  and  (8.13), 

7  =  0(l)6*(i  +  l)-i^(x,t). 

Therefore,  (8.21)  is  proved  for  1  =  0.  For  l  =  1,  (8.21)  follows  by  integration  by  parts. 

□ 

Lemma  8.4.  Under  the  assumptions  of  Theorem  2.3,  we  have  for  — oo  <  x  <  oo,  t  >  0, 
1  =  2,  ...  ,5, 

/oo  cd  t 

(G  -  G*)ij(x  -y,t)—lu(y,0)dy  =  0(1)6* (t  +  l)-1"*. 

Proof.  Similar  to  the  proof  of  (8.20),  (8.24)  is  clear  by  integration  by  parts  and  (8.12). 

□ 


From  (8.3),  (8.8)  and  (8.9),  we  have  the  following  expressions: 


(8.25)  Qi  —  9i  —  ^ iui  —  hf(u)  +  ri@i) 

=  -hi'ZrmrjOj^e,)  -  li'Ef'WMiM) 

s  s 

+  0(1)  J2  |0jl(kl  + 1^1)  +  0(1)  E  N2  +  E  ^-l3 

j/fc  j- 1  3= 1 
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{9i~9t)x  =  +  Y  (NM 

j,k= 1 

+  IMM  +  |0,flM  +  Mk*|)  +  o(i)  Y \0jx\\0k\, 

&k 

cf  r  &L. 

=0(1)  |K,t)||L^  —(.,«)  +\\ux{;t)\\2W^.~ 

L°°  L  Udj  1,00 


Let  Oa  be  defined  by  (3.8): 

(8-28)  0<y(x,t;  A,/x)  =  (t  +  1)~§  e~(  mu++i)))  ? 

where  a  >  0,  (x  >  0  and  A  are  constants. 

Lemma  8.5.  Under  the  assumptions  of  Theorem  2.3,  we  have  for  -oo  <  x  <  oc,  t  >  0, 

1  =  0 ,1, 


(8-29»  u 


f°°  1  (x-y-AjCt-r))2  Ql- j-1 

T>  dy^('9ij~9^^y,T^dydT 


~  0(1)(<5*  +  M{t))2(t  + 


Proof.  Denote  the  left-hand  side  of  (8.29)  as  /.  By  (8.25),  (8.26),  (3.1),  (8.15),  (8.13)  and 
integration  by  parts,  we  have  for  l  =  0,  1, 


a. 


(x  —  y-Xiit-T))- 


OC  y/AnqMYt  -  r) 

dl+1  f  1 


Y.  f"(0){rk0k,  rk0k)(y, t)]  dydr 


+  0(1) 


dyl+1  I  2 


(t  -  [*.  (tf.  +  M(t))  05/2(y,  r;  Xk,y) 

k 


+  M(r)2  y^(r  +  1)  4  ^3/2(y,  r;  A*)]  dydr 


+  0(1)  jt  j  (f-r)-V  ^)(->  [tf*(6*  +  M(r))^^5/2+I(y,r;Afc,/i) 

2  °°  fc 

+  M(t)2  ^(t  +  1)“ 4  “ 2^372(2/, T;  A*.)]  dydr 


==/i  +  /2  +  /3, 
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where  e,  p  >  0  are  constants.  We  see  that  the  first  integral  h  in  (8.30)  is  II  in  (4.14) 
when  l  =  0,  and  II  in  (4.16)  when  l  =  1.  Hence  by  (4.27)  and  (4.33), 


(8.31) 


o 

h  =  0{l)6*2(t  4-  l)-*  '^2‘ip3/2(x,t-,\k). 

k—1 


To  estimate  the  second  integral  in  (8.30),  we  use  Lemmas  3.2,  3.3,  3.5  and  3.6: 

ft/2  /»oo 


ft/ 2  fOO  (x-V-A,;(t-r))2 

=  0(l)(<5*+M(f))2rJ  /  /  (t -r)_1e 

JO  «/ — oo 

■  ^[^5/2 (?/i  r5  A*,  a)  +  (r  +  1)  A'$z/‘2.{y,'T\\k)\  dydr 

k 

=  0(1)(«*  +  M(t))2H  YJ&3/ At)  +  ^i(*,t;A*)V>3/4(*.t;^fc)]. 


(8.32) 


where  we  have  used  (3.24)  and  the  fact  that  for  >  |A i  —  Afc|, 

\x  -  A i(t  +  1)|"S |rc  -  A k(t  +  1)|“*  char { min (Ai,  A *)(f  +  1)  +  -K'x/t  +  1  <  £ 
(8.33)  <  max(Ai?  Afc)(t  +  1)  -  K\Jt  +  1} 

—  0(1)  [0i  {x,  t;  Ai)^3/4(^>  A*)  +  0i(x,  ^  Afc)V'3/4(^5  At)]. 

The  third  integral  in  (8.30)  is  estimated  in  the  same  way, 


(8.34) 


J3  =  0(1)(<5*  +  M(f))2(t  +  l)“4$(x,  *). 


Hence  (8.29)  is  proved  by  (8.30)-(8.32)  and  (8.34).  D 

Lemma  8.6.  Under  the  assumptions  of  Theorem  2.3,  we  have  for  -00  <  x  <  00,  t  >  0, 

1 


t  roo 


<“»  LI 


e  —rvr(gij  ~  9y)(v, 0  dydr 


dyl+1 
=  0(1)(6*  +  M(t))2 


+  if  l  =  2,3 

(t  +  1)"2,  if  1  =  4,  b 


Proof.  Denote  the  left-hand  side  of  (8.35)  as  7.  By  integration  by  parts, 
■t/2 


7  = 


/•t/2 

0(1)  /  (t-T)-5_1H(sy -fly)(-,T)IU‘^ 

Jo 

+  0(1)  f  ( t 

Jt/2 


dyl 


(dij  9ij)(’i  T . 


dr. 


L°° 
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By  (8.25),  (3.1),  (8.15)  and  (8.13)  it  is  clear  that 


/oo 

[  +  M(t))(t+ 1)_  I  \6k ( y ,  r)|  +  (M(r)$(j/,  r))2]  dy 

-OO  7. 


0(1)(6*  +  M(t))2(t  +  1)-*. 


Also  by  (8.27),  (8.1),  (8.9),  (8.8)  and  Lemma  2.2, 


(9ij  -  9ij)(-,T~)  —0(1)  (y,  II Ok  +  (y "2  fyjl(0k+vk)  ^ 

(\  ^  ol 

ll®*»  +  "fcj/llw1-2.”)  +  Q^if"(ty(riOi’ri0i) 


=  0(l)(S* +M(r)) 


2f(r  +  l)-3/2,  if  1  =  2,  3 
1  (r  +  l)-1,  if  l  —  4,  5 


Thus  (8.35)  is  straightforward. 
By  (8.3)  we  have 


hij  =  lij(B(u)  -  B(0))ux  =  bij(u,ux)  +  0(l)|w|2|ux| 

S 

—  ,  bij(rk6k,  Tk^kx)  T  0(1)  y  ^  |$A;  |  life's  | 

k= 1  fc#fc' 

+  O(i)  5^[|Vfc||^fe/*|  +  |0fe||v*'*|  4-  \vk\\vk'x\  +  |<9fc|2|^x|], 


where  bij  is  a  bilinear  form.  Similarly, 

(8.37)  hijx  =  0(l)(\u\\uxx\  +  \ux\2) 

=  0(  l)fe(|ff*||fl*.«|  +  |  l^fc'xrr  |  H"  |^/c  | \^k'xx  |  "h  \Vk  \  l^k': 


+ ym*\2 + h;fcxp) , 

JU 


dl  f  dl+1 

loo  =  o(1)[N,»0IU“ 


~Hi^c(*J  £)  Hl00  o  /W(*>^)  d"  [|ttx(*j 


,  /  >  2. 
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Lemma  8.7.  Under  the  assumptions  of  Theorem  2.3,  we  have  for  —oo  <  x  <  oo,  t  >  07 


(8.39) 


(x-y-iMt-r))2  Ql  +  l 


y/4irfJLij(t  -  r) 


iViiO-T) 


Qyl+l 


hij(y,r)  dydr 


=  0(1)(<5*  +  M(t))2< 


(t  +  l)"1, 
(t+  l)“i, 


if  1  =  0,  1 
ifl  =  2,  3, 
i/I  =  4,  5 


where  a  >  0  is  arbitrarily  small.  Moreover,  if  all  the  A k,  1  <  k  <  s,  are  simple,  a  =  0. 


Proof.  Denote  the  left-hand  side  of  (8.39)  as  I.  For  l  =  0,  1,  we  have  by  (8.36),  (8.37), 
(3.1),  (8.15)  and  (8.13), 


t  r  OO 


(x- j/-Aj(t-r))< 
(t-r) 


OO  y/4irfiij(t  -  r) 
t—  1  fOO 


(8.40) 


-a 

♦<**> n 

+  m(t)2(t  + 1)~*  Y!  ^3/2(1/.  t;  a*) 


dl+- 


I  _  (ac-y-^t(t-r)) 


/>t  />00 

0(1)/  /  (I-r)-1. 

Jt—lJ  —OO 


+ 


(x  —  —  r))' 

4(Mil+£)(t-r) 


^  jj^f[bij(rk8k,rkdkx)(y,T)dydT 

«*(«*+M(r))^flj(»,T;Afc,|i) 

k 

dydr 

6'(6*  +  M(r))J2e3(y,T;  Xk,it) 

dydr 


+  M(r)2(r  +  1)  1  y^(^3/2(y^;^fc)  +  Afe)^3/4(3/,r;  Afe)) 

k 

Ii+  h  +  h, 


where  e,  p  >  0  are  constants.  To  estimate  Ji,  we  use  Lemma  3.2  and  integrate  by  parts 
for  k  —  i,  and  use  Lemma  3.4  for  k  ^  i.  Thus 


h  ~  0(1)6*2  \$2 +i(x,t\  A j,  p)  log(t  +  2)  +  +  1)  2  1  +0f V’3/2 ^0] 

k 

=  0(l)6*2(t  +  l)-*-*+a^2*p3/2(x,t;\k). 

k 


QUASILINEAR  HYPERBOLIC-PARABOLIC  SYSTEMS 


103 


Then  we  apply  Lemmas  3.2,  3.3,  3.5  and  3.6  to  I2 : 


h  =  0(1)(£*  +  M(t))2<  (t  +  1)  2  f  f  ( t-T)~le 

l  J0  oo 

'  +  (r  +  l)_^^3/2 (*/,£;  Afe)]dj/dr 

k 

+  (t  +  l)-1/*  f  (t  -  r)_1(£  +  1  -  r)“ie_ ' 

it/2  i-oo 

•  X)[03/2(y,  t;  Afc,  M)  +  (r  +  1)_*  ^3/2(y,  t;  A*)]  dydr  j 
=  OC1)^*  +  M(t))2(t  +  l)-*-*  ^2ip3/2(x,t;Xk). 

k 

To  estimate  /3,  notice  that  (4  +  1  -  r)  =  0(1)  for  t  -  1  <  r  <  t.  Thus 
h  =  0(1)(6*  +  M(t))2(t  +  1)-?  I*  ^  (t  -  T)~\t  +  1  -  T)-V^I+«2 

it-1  i-oo 

•  5^[^3/2(j/,  t;  A*,  At)  +  (r  +  l)"i  ty-s/afe,  i;  A*)  +  <Ai  (2/,  r;  \k)i>3/4(y,  r,  Afc))]  dydr 

k 

=  0(1)(<5’  +  M(t))2(t  +  l)-t  yyVsMM;  Afc)  +  ^(x,  t;  \k)4>3/4 (y,T;  Afc)) 

k 

by  Lemmas  3.2,  3.3,  3.5-3.S.  We  have  proved  (8.39)  for  l  =  0,1.  We  see  that  a  >  0 
appears  only  in  Ix.  If  all  the  \k  are  simple,  then  the  0k  are  scalars,  and 


bij{rkdk,rkekx)  =  bij(rk,rk)^Y^  ■ 


We  can  avoid  a  in  this  case.  For  2  <  /  <  5,  we  have 

ft_1  l  (tc  —  y  —  Aj(t  —  t))2 

1  =  0(1)  /  (t-r)-5-ie  4().ij+«„t-r)  hijiy^dydT 

•JO  J  —  oo 


+  0(1) 


a  00  !  (s-T/-A,-(t-T))2 

J‘-r)  e  M„,r)(WT 


rt-l 

0(1)  /  (t  - 
J  0 


S/4  +  /5. 

Prom  (8.36),  (8.15), 


T)  2  2ll%(-!'f)||L~dr  +  0(l)  f  (t-r)  2  S-ihij(-,T)  dr 

Jt- 1  0JT 


ft- 1 

h=0(  1)  /  (t  —  r)_2_i((5*  +  M(t))2(t  +  l)_i  dr 
0 

J1-  /'t/ 2  _ 2 

I  (r  +  l)-§dr  +  (t  +  l)-f  /  (t-r)-i-idr 
0  it/2 

=  0(1)(«*  +  M(t))2[(t  +  l)-i-l  +  (t  +  i)-l], 
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which  is  bounded  by  the  right-hand  side  of  (8.39).  By  a  standard  interpolation  formula 
and  Theorem  7.1,  we  have  for  t  >  0, 

(8.41) 


Together  with  (8.38)  and  (8.15), 


qxM-’V 


=  0(1)(«*  + 


{t  +  1)' 
(t  +  1)' 


if  l  =  2,  3 
if  l  =  4,  5  * 


Thus 


h  =  0(l)(6*  +  M(t))2 


(t  +  1)  \ 

(«  +  !)"** 


if  l  =  2,  3 
if  i  =  4,  5  ’ 


which  is  also  bounded  by  the  right-hand  side  of  (8.39). 


□ 


Lemma  8.8.  Under  the  assumptions  of  Theorem  2.3,  we  have  for  -oo  <  x  <  oo,  t>0, 


pt  pOO  +  l 

(8.42)  /  /  (G -G*)«(*-»,t- +  /»)(», r)d»dr 

f  (t  +  l)“5^4>(a:,t),  ifl  =  0,1 
=  0(1)(£*+M(f))2i  (t  +  l)-1,  if  1  =  2,  3. 

[(t  +  1)-*,  if  1  =  4,  5 


Proof.  Denote  the  left-hand  side  of  (8.42)  as  I.  Consider  l  =  0,  1.  We  have 

I  =  h  +  h  + 13, 


where 

pt  poo  _ _ .  1 

h=Jo  j  (G  -  G*)ij(x  -  y,t  -  t)q^Th  yy--L/"(O)(rfc0fc,rfc#fc) 

+  b{rkek,rk0kxj) (y, t)  dydr, 

I2  =  J*  J°°  ( G  -  G*)y  ( x  -y,t-  (5  +  \Lf"(0){rk9k,rk8k)j  (2/,  r)  dydr, 

pt  poo  +  l  ^ ^  \ 

/3  =  y  J  (G-G*)ij{x-y,t-T)-^i^[h-22b{rk9k,rk6kx)j(y,T)dydT, 

and  b(w\,W2)  is  the  n-vector  valued  bilinear  form  of  rui,  W2  £  Rn  which  has  bij(w  1,102)  85 
the  jth  component  of  its  ith  block.  To  estimate  h,  we  integrate  by  parts  with  respect  to 
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x  for  0  <  r  <  t/ 2  and  use  (8.12): 


ft/2  POO  (x-iz-AfeCt-x))2 

=  /  /  0(l)(*-r)  ^  6*202(y,T-,\k,,n)dydT 

Jo  J~°°  k,k>= i 

p£/2  m'  s 

+  /  0(l)5*2e-(t-r)/c  EE0  2 {x  -  fait  -  t),  r;  A*. ,  /u)  dr 

1,0  *=ifc'=i 

pt  POO  3  2 

+  /  /  0(l)((-r)-i(t  +  l-r)-lj]  e-— 0(*-‘rT)>  6t2e3+l(y,T;\k,,n)dydT 

pt  m!  s 

+  /  0(l)6*2e_(t_T)/c  O3+i{x-0k(t-T),r;Xkl,fi)dT 

Jt'2  k=ik'=i 

=OW2(t+l)-*\  ±  /*  r(t-r)-Vls=m^Wyfr;W)^ 

m'  s 

+  EE  /  e~(t-T)/c02(x  -  hit  -T),r\\k',n)dT  , 

k=i  k'= i  -70 


where  M  >  0  is  a  constant.  Applying  Lemmas  3.2,  3.3  to  the  first  term,  Lemma  3.9  to  the 
second  term,  and  using  (8.33),  (3.24),  we  see  that  A  is  bounded  by  the  right-hand  side  of 
(8.42).  By  (8.3),  (8.9),  (3.1),  (8.15),  we  have 

d2n 

0^2  =  0(l)(|u||u«,|  +  W2) 

=  0(l)(tf*  +  M(f))2  [£  93(x, t;  A fc,  A*)  +  (t  +  l)"1*^,  t)l , 

Hfe=l 

O^h 

(M  l^arxx |  d*  |wx||'W;E;E|  + 


0(1)(6*  +  M(t))2  ^03(M; Xk,fi)  +  {t  +  l)_13>(x,t) 


Together  with  (8.26),  (8.8),  (8.36),  (8.37),  and  making  use  of  (8.12),  we  estimate  I2  and 
Is  in  a  similar  way  as  I\ : 

ft  POO  3  2 

h  =  /  /  0(l)(t  -  T)-i-i(«  +  1  -  r)-i  Y  e-l'-rvXT”  (6*  +  M(r))2 

Jo  J~°°  kpi i 

'  [07/2(2/, t;  Afc/,^)  +  (r  +  l)~*ip3/2(y, r;  A*,-))  di/dr 

ft  m'  p  S 

+  /  0(l)e-(t-T)/cy](A*  +  M(T))2  j]fls(ic-/?t(«-r),T;Alk>,|*) 

-70  *=i  U'=i 
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+  (r  +  l)  X^(x-  Pk{t-r),T) 


dr , 


(x  —  y  —  A l.  ( t  —  t ) _ 

e  (<5*  +  M(r))2 


/*t-l  /*oo  JL 

h=  [  /  ofiKt-r)-1--  ^ 

Jo  J  oo  fc,fc'=l 

■  [67/2{y,T-,  \k',v)  +  (r  +  1)- hs/2 (y,T;Xk')]dyd.T 

rt  m'  r  5 

+  /  0( l)e-(t“T)/c,]T(l5* +  M(r))2  ^03(x-/?fc(t-T),T;Aifc.,/*) 

■'O  fe=l  U'=l 


+  (t  +  1) 

t  r  oo 


dr 


+ 


t—  1  J  —  oo 


1  _  (x-y-Afc(t-T-))^ 

0(l)(i  -  t)~"2  (t  +  1  -  r)“2  c(*-r)  (6*+M(r))J 


fc=i 


-  s 

D  d3(,y,T-,\k',fJ.)  +  (r  +  l)_1$(j/,r) 

-fc'  =  1 


dydr. 


The  first  term  of  I2  is  estimated  in  the  same  way  as  in  the  proof  of  Lemma  8.7,  see  /2  and 
/3  in  (8.40).  The  second  term  of  I2  can  be  written  as 


0(1)(<5*  +  A/(t))2  Y  /  e~{t~T)lc{r  +  1)  4$(x-/3jc(t-T),T)dT 

frt  do 


m' 

=  0(l)(tf*  +M(t))2(t  +  iyi  Y"  /  e-(t"T)/c$(a:-/3fc(t-T),r)dr: 

fc=iJo 


which  is  bounded  by  the  right-hand  side  of  (8.42),  using  Lemmas  3.9  and  3.10.  Clearly  J3 
is  estimated  in  the  same  way. 

For  2  <  l  <  5,  we  have  by  (8.12) 


/  = 


t/  2  /*oo  ^J+l 


—  oo 

t  rOO 


5xz+1 


(G  -  G*)ij(x  -  y,t-  r)g(y,  t)  dydr 


rl  r°°  ft  dl 

U_ o^G^G^x-y’t~T)w9{y'T)dydT 

+  J  J  T (G-G*)ij(x-y,t-T)hy(y,T)dydT 

rt  POO  ft 

+  J  J  ^(G-G*)ij(x-y,t-T)—lh(y,T)dydT 

rt/l  /*t/2 

=  J  0(l)(i-r)-£^||fl(-,r)||Lidr  +  ^  0(l)e-(‘-T)/c||s(^T)||w.+i.~  dr 
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+  /  0(l)(t-r)  "(t  +  l-r)  2  dT 

Jt/2  oyl  Loc 

+  f  0(l)e-^)/C  JLg(.,T)  dr 

Jt/2  oyl  ^1,00 

ft-  1 

+  /  0(l)(t-r)-*(t  +  l-r)“*||V-,T)|Uoodr 

./o 

+  [  0{\)e-^yc\\hy{;r)\\w^dr 

Jo 

+  [  0(l)(t-r)-i(*+l-r)-i  dr 

./t-i  oy  Loo 

+  f  0(l)c-<‘-T>/°  dr. 

it-l  ^2/  M/l.oo 


,-(*-r)/C 


dr. 

1^1,00 


By  (8.25),  (8.8),  (8.27),  (8.41),  (8.38),  we  have 

||<?(-,r)||Ll  =  0(1)(6*  +  M(r))2(r  +  l)-*, 

di  (r  +  l)“3,  1  —  2,3 

Loo=°m*  +  M(T))>}  (r  +  l)-1,  f  =  4,5, 

,  (r  +  l)-^,  /  =  6 

di  [  ('7_  +  l)“i)  *  =  1 

^(‘,r)  L»=0(1)(^  +  M(r))  1  (T+1)-1’  *  =  W 

l  (r  +  1)"*,  /  =  5,6 

Thus  I  is  bounded  by  the  right-hand  side  of  (8.42). 

From  (8.11)  and  (8.14),  applying  Lemmas  8. 1-8.8,  we  conclude  that 


M{t)  <C6*  +C{6*  +M(t))2. 


If  M(t)  and  6*  are  small,  we  have  M(t)  <  C6*.  By  continuity,  M(t)  <  CS*  if  <5*  is  small. 
Together  with  (8.15),  (2.25)  follows. 

We  have  proved  Theorem  2.3.  From  Theorem  6.15  and  Remark  6.4,  we  can  find  explicitly 
Qko  in  (8.12).  It  is  easy  to  see  that  if  the  left  eigenspace  of  B  associated  with  the  zero 
eigenvalue  is  independent  of  u ,  then  QkoLB(u)  =  0  for  all  it,  where  L  is  the  matrix 
consisting  of  the  /*.  Thus  when  we  perform  a  priori  estimate  for  dlVij/dxl  using  (8.11),  the 
highest  derivative  of  v  involved  in  the  source  terms  is  dl+1v/dxl+1  instead  of  dl+2v/dxl+2. 
To  obtain  the  pointwise  estimate  for  v(x,  t),  it  is  sufficient  to  have  an  L°°  bound  for  vx(x,  t) 
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rather  than  a  pointwise  estimate.  In  fact,  it  is  sufficient  to  require  that  ||u*(-,£)||l~  decays 
as  (t  +  I)"1,  K*(-,t)|U“  as  (t+l)-i,  and  ||vXM(-,t)||L~  is  bounded.  Therefore,  in 
Theorem  2.3  it  is  sufficient  to  require  uo  €  Hi{ R)  rather  than  Uo  £  Hs( R).  And  this  is 
Remark  2.5. 

We  now  outline  the  proof  of  Theorem  2.6.  On  the  right-hand  side  of  (8.11),  the  first 
two  terms  give  the  contribution  of  the  initial  data.  In  (2.27)  we  have  improved  the  decay 
rate  of  w0  and  u'0  to  \x\~i  as  x  oo.  These  two  terms  then  become 

(8.43)  0(l)£*(f+  l)-* 

k=  1 


for  Z  =  0,  1,  where  i/S  is  defined  in  (2.19),  c.f.  (4.14),  (4.22).  Define 


(8.44)  M(t)  ~  sup  max 

V  ;  0<r<t 


(||ut(-,r)$i(-,r)-1||i,t»  +  |M-,r)(r +  !)*$<(•,  r) 


+ 


d 2 


3  / 

s( 

1=2  X 


dx 2 
d21-1 


(r  +  1) 


dx21' 


L°° 

lvi('>T) 


9 21 


dx21 


(T  + 


l)^ 


where 


(8.45)  =  ip3/2(x,t;Xi)  +  '^2'4>(x,t-tXj). 

jV* 

Similar  to  Lemma  8.8,  we  can  show  that  the  fourth  term  on  the  right-hand  side  of  (8.11) 
is 


(8.46)  0(1)(^+M(t))2(t  +  1)  ’$i(M),  *  =  0,1. 


The  third  term  has  the  same  bound  (8.46).  The  proof  is  similar  to  the  proofs  of  Lemmas 
8.5  and  8.7  except  that  we  need  to  apply  Lemma  3.4  to 


<8  i7)  17 


(X-V-Xjjt-T))4 
4  Hij(t-T) 


0  j  —  OO  v4 7TMij(t  ~  r) 


dl+1 

dyl+l 


-lij  ^  f"{0){rk8k,rkVk){y,  r) 

k^i 


dydr ,  l  =  0, 1. 
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in  (8.29),  and  to  I\  in  (8.40).  To  estimate  (8.47),  we  require  ||vio:x(*>  £)IIl°°  =  0(l)M(t)(t  + 
l)~i.  In  turn  we  need  uq  G  H9( E)  to  close  the  a  priori  estimate. 

We  see  that  in  (8.30)  the  contribution  of  the  nonlinear  source  term 
r  1 

frJ+1  ~2 lijJ^f"(°)(r^k,rk0k)(t/,T)  ,  1  =  0,1, 

^  k^i 

gives  the  rates  (t  -f  l)-3/2  and  (t  +  l)-2  for  Vi  and  ViX,  respectively,  away  from  the  char¬ 
acteristic  directions.  These  two  rates  can  not  be  improved. 

We  now  outline  the  proof  of  Theorem  2.9  when  the  linearization  of  (2.1)  is  not  diagonal- 
izable.  By  the  assumption  B  is  a  nonsingular  constant  matrix.  Thus  in  (8.11)  h  =  0,  and 
in  (8.12)  the  <5-functions  disappear.  For  l  =  0,  1,  the  first  and  third  terms  in  (8.11)  were 
estimated  in  Section  4,  see  (4.13)-(4.16),  (4.31)  and  (4.35).  The  second  term  is  estimated 
in  a  similar  way  as  the  first  term,  making  use  of  (8.12),  and  the  bounds  are 

s 

0(1)6*  and  0(l)£*(t  +  1)“*, 

fc= 1 

respectively,  for  l  =  0,  1.  The  fourth  term  is  estimated  in  the  proof  of  Lemma  8.8,  with 
h  =  0,  6  =  0,  and  the  bounds  are 

s 

0(l)(r  +  M(t))2y^(x,f;Afe)  and  0(1)  (6*  +  M(t)f(t  +  1)"*, 

fc  =  l 

respectively,  for  l  =  0,  1.  Here  we  notice  that  we  only  need  the  L°°  norm  of  the  Vjx  instead 
of  the  pointwise  estimate,  since  B  is  a  constant  matrix,  and  that  the  VjXX  do  not  appear 
since  we  no  longer  have  ^-functions.  The  property  (2.32)  then  follows  from  (4.31)  and 
(4.35). 

At  last,  we  point  out  that  the  solution  u  to  (2.1),  (2.2)  is  approximated  by  the  solution 


ut  +  f{u)x  =  B*uxxt  -oo  <  x  <  oo,  t  >  0, 
u(x ,  0)  =  wo(^)?  —  oo  <  x  <  oo 


at  a  higher  rate  than  by  the  diffusion  waves,  where 


B*  =  (r,,  ■■  ■  ,rs) 


hB(0)n 


lsB(0)rs )  \ls 


110 


TAI-PING  LIU  AND  YANNI  ZENG 


In  fact,  let 


(8.50)  u* (x,t)  =  riU*(x,t),  v‘(x,t)  =  Ui{x,t)  -  u*(x,t). 

i- 1 


Then  similar  to  (8.11)  we  have 
(8.51) 

/OO 

(G  -  G*)ij(x  -  y,t)u(y,0)  dy 

-OO 


t  r  oo 


+ 


a 


oo 
t  poo 


1  _  (s-y-AjU-r))2  f) 

—7====e  4*“i(t'T)  ^-{9ij(u)-g*j{u*)  +  hij(u))(y,T)dydT 

V4t Tinjit-T)  oy 


r  r°°  d 

+  J  J  (G  -  G*)ij(x  -  y,t  -  t)  —  {g{u)  +  h(u)){y,r)dydr, 
j  =  1,  •  •  •  ,nrii,  i  -  1, . . .  ,s, 


where  g  and  h  are  given  by  (8.3).  Except  the  contribution  from  gij(u)  -  gij{u*)  in  the 
second  term,  the  right-hand  side  of  (8.51)  has  been  estimated  in  Lemmas  8.3,  8.4,  8.7  and 
8.8;  it  decays  faster  than  Vij(x ,  t)  by  a  factor  (t  +  l)“*+e,  where  £  >  0  is  arbitrarily  small. 
Then  by  a  priori  estimate  we  have  the  final  remark  of  this  section: 

Remark  8.9.  Under  the  assumptions  of  Theorem  2.3,  the  solution  u  to  (2.1),  (2.2)  is 
approximated  by  the  solution  u*  to  (8.48),  with 


(8.52)  Ui(x ,  i)  -  u*(z,  t)  =  l%u(x,t)  -  ku*(x ,  t) 

=  0(1)6 *  [(*  +  l)-i+£$(x,  t)  +  (t  +  l)~i\x  -  A i(t  +  l)|“1+£ 

■  char{(i  +  l)2/£  <  lx  “  A*(t  +  1)|  <  e(t  +  1)}], 

where  e  >  0  is  arbitrarily  small.  From  (8.52)  we  see  that  u{x,t)  -  u*(x,t)  decays  as 
( t  +  l)~1+£  along  the  characteristic  directions,  and  as  ( t  +  l)“3/2+£  away  from  those 


directions. 


9.  Applications  to  Continuum  Mechanics 


We  now  apply  our  general  result  to  two  important  examples:  the  Navier-Stokes  equa¬ 
tions  and  the  equations  of  magnetohydrodynamics. 


Example  9.1.  The  Navier-Stokes  equations  describe  the  classical  1-D  compressible,  vis¬ 
cous  flow  or  1-D  thermoelastic  flow.  In  Lagrangian  coordinates  the  equations  read: 


(9.1) 


vt-ux  =  0, 

<  Ut  +  Px  —  (^~WX)X, 

>  (e+h2)t  +  (pu)x  =  (*ex  +  *uux)xi 


where  v>  u,  p,  e  and  0  are,  respectively,  the  specific  volume,  velocity,  pressure,  internal 
energy,  and  temperature  of  the  fluid.  For  simplicity  we  assume  that  the  viscosity  p  is  a 
nonnegative  constant  and  the  heat  conductivity  k  is  a  positive  constant.  Besides,  we  have 
the  equation  of  state.  Denote  the  entropy  by  s.  Then  among  the  thermodynamic  variables 
v,  p,  e,  0  and  s,  we  can  choose  two  of  them  as  independent  variables  and  regard  the  rest 
as  smooth  functions  of  them.  We  write 


p  =  P(v> e )  =  P(v,  s)  =  p(v,  0),  0  =  0(u,  e)  =  0(v,  s), 

(9.2) 

s  =  s{v,  e)  =  g(v,  0),  e  =  e(u,  s)  =  e(v,  0), 

and  assume 


(9.3) 


Pv  =  foP(V’  d)  < 


0e  -  Q-0(v,e)  >  0, 


Pe  =  fcP(V<  e)  ^  0. 


The  thermodynamic  law  de  =  6ds  -  pdv  implies 


(9.4) 


ev  =  - p ,  es  =  6,  eg  =  6se , 

ev  =  8pe  -p,  se  =  1/6,  sv  =  p/6. 


From  (9.2)  and  (9.4)  we  have 


Pv=Pv-  PPe,  6v  =  -Ps  =  =  -6^-, 

se  ee 

(9.5)  pv  =  pv  +pe0v  =  pv  -  pv  =  Pv  -f-  evpe  =pv- 

ee  0e 

0  1 

Pe  =  PsSe  =  -j  =  -g(6v-pee),  6v=6v+p6e  =  —ev6e. 


Ill 
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To  examine  that  Assumptions  2. 1-2.3  are  satisfied,  we  rewrite  (9.1)  in  the  form 


(9.6) 

Here 

(9.7) 

and  we  have 

(9.8) 


w 


Wt  +  f(w)x  =  (B(w)wx)x. 


=  (v,u,e+  -uy,  f(w)  =  {-u,p,pu)\ 


0 

B(w)  =  I  0 


-1 


A(w)  =  f'(w)  =  I  Pv  -UPe  Pe 


\Upv  p  U  pe  Upe  J 

We  consider  w  around  a  constant  state  w*  —  (v*,  0,  e*)  with  v*t  e*  >  0.  Thus  u,  /'( 0)  and 
B(0)  is  Section  2  correspond  to  w  —  w*,  f(w*)  and  B(w*)  here.  For  Assumption  2.1,  we 
take  U(w)  =  -s  and  F(w)  =  0  as  the  entropy  pair.  Thus  by  (9.4), 


VC/  —  (  sv,'iise,  se )  ^(  Pj w,  1), 


(9.9) 


(svv 

USye. 

~Sye 


Together  with  (9.5)  and  (9.3), 


det  Aq  —  se(svvsee  sve)  — 


USye  Sve 

Sg  XL  See  USee 
XlSee  See 


-pvde  +  (p0e  -  Qv)~J 


—  04  (  +  Pe@ v)  —  ^4  @ePv  ^  9, 


(9.10) 


— svv  =  — det  A0  +  szvc  )  =  i-O'^e det  A>  +  sje)  >  0, 

See  \  Se 


K 


det  I 

usve  S e  XL  See 


)- 


U * 


=  —svvse  H - det  Aq  >  0. 

Se 


Therefore  Aq  is  positive  definite  and  U  is  strictly  convex  with  respect  to  w.  It  is  clear  that 
(VC/)/'  =  0  by  (9.8)  and  (9.9).  From  (9.7)  and  (9.9), 


j  f  k>Q%)< 


,Sve  KXL0eSve  KiO  e  Sye 

AqB  =  -  I  Ku0vSee  /iSc  “  KU20eSee  Ku9eSee 
\  — hzQySee  fcXlOpSee  f^OpSee 
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which  is  symmetric  since 


(9.11) 


OeSVe  —  @v$ee  —  — . 


It  is  straightforward  to  varify  that  A0B  is  semi-positive  definite  by  looking  at  the  minors. 
For  Assumption  2.2,  we  let 


(9.12) 


w  =  (v,u,6)t,  w  =  f0(w). 


Hence 

/I  0  0\  /0  0  0  \ 

f'0{w)  =  0  1  0  ,  B=  0  £  0  , 

\e„  u  eg)  \0 

where  we  have  used  (9.5).  Clearly  the  null  space  Af  of  B  is  independent  of  w  or  w. 
(However,  the  null  space  of  B  depends  on  w.)  We  see  that  Af1-  is  spanned  by  {773}  or 
{'12,  r] 3},  depending  on  n  =  0  or  n  >  0,  where 


%  =  (0,1,0)*,  %  =  (0,0,1)*. 

Thus  is  invariant  under 

Svv  $ve^v  0  0  \ 

o  Jo). 

SveC$  USeeCQ  — SeeCQ  J 

It  is  clear  that  B  maps  R"  into  Af1 .  Assumption  2.3  has  already  been  varified  in  Sec¬ 
tion  6.2.  We  notice  that  the  left  eigenspace  of  B(w)  associated  with  the  zero  eigenvalue  is 
independent  of  w.  By  Theorem  2.3  and  Remark  2.5,  if  the  initial  data  w(x,  0)-w*  e  H4( R) 
and  satisfies  the  smallness  assumptions  in  Theorem  2.3,  then  the  solution  to  the  Navier- 
Stokes  equations  has  a  time  asymptotic  state  as  three  diffusion  waves  along  the  directions 
dx/ dt  =  Ai,  1  <  i  <  3,  where 


Xl  =  -V-Pt»  >2=0,  A3  =  v^* 

are  the  eigenvalues  of  A(w*),  p*v  denotes  the  value  of  pv  at  the  constant  state.  The  first 
and  the  third  diffusion  waves  are  nonlinear,  corresponding  to  genuinely  nonlinear  fields  of 
the  inviscid  system,  while  the  second  one  is  linear,  corresponding  to  a  linearly  degenerate 
field.  The  difference  between  the  solution  and  the  diffusion  waves  is  then  given  by  (2.25). 
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Example  9.2.  [LLP]  In  Lagrangian  coordinates,  the  equations  of  magnetohydrodynamics 
characterizing  the  flow  of  a  conducting  fluid  in  the  presence  of  magnetic  field  are 


(9.13) 


'  vt  -  uix  =  0, 

Ult  +  [P  +  2f^(B2  +  Bi)]x  “  (%ulx)x, 

U2t  -  =  (%U2x)x, 

u3t  ~  (~j^BlBz)x  ”  (^3x)x> 

<  {vB2 )t  ~  (BiU2)x  ~  {^^B2x)xi 
(vBs)t  -  ( Blu3)x  =  (~^B3x)x, 

(e+l(u2  +  u2  +  u2)  +  ^v(B2  +  B2))t 

+  [(P  +  2^(B2  +  B3))Ul  ~  mB'(B2U2  +  B3U3^x 

k  =  ["tXi1iix  +  ^(U2U2X  +  ^3W3x)  +  ^x  +  +  B3B3x)]x, 


where  v,  u  =  (tii,ii2j  «3),  P,  B  =  (BuB2 ,B3),  e  and  6  represent  the  specific  volume, 
velocity,  pressure,  magnetic  induction,  internal  energy  and  temperature  respectively,  p  = 
p(v,e),  0  =  0(u,e),  B*  is  a  constant,  /x,  i/  >  0  are  the  two  coefficients  of  viscosity,  p0  >  0 
is  the  magnetic  permeability,  k  >  0  the  heat  conductivity  and  1/cr  >  0  the  electrical 
resistivity.  Again,  we  have  the  equation  of  state. 

Denote  the  entropy  by  s.  Introduce  the  notations  in  (9.2)  and  assume  (9.3).  We  have 
the  thermodynamics  relations  in  (9.4)  and  (9.5).  Denote  the  total  energy  and  full  pressure, 
respectively,  as 


(9.14) 


E  —  e  4-  ^(ti?  4-  u\  +  u\)  +  —  (B2  +  BVh 

P  =  P+^-(Bl  +  Bj). 

ZfiQ 


We  write  (9.13)  in  the  form  of  (9.6), 


Wt  +  f(w)x  —  (B(w)wx)x-> 


with 


W  =  (v,Ui,U2,U3,vB2,  vBz,  E)1 , 

f(w)  =  {-ul,P,-—B[B2,-—BlB3,-B\u2,-B*luz,Pu1  -  j-B\{B2u2  +  B3u3))' 
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(9.15a) 


where 


(9.156) 


Clearly, 


(9.16a) 


where 


B(w)  = 


0 

0 

0 

0 

b2 

_ 


afiov^ 


anov< 

0 


a/iov* 


\  ai  a2ui  a3u2  a3u3  a462  04.63 


°\ 

0 

0 

0 

0 

0 

6/ 


b  = 


_L  (T&  _L  U^\(  b  1  \ 


«2 


1 

Mo 


c7  =  -peU\U3 


M  , 

b 

1 

=  -  -  o,  a3 

=  ~  “  6, 

a4 

+  09 

V 

V 

Mo 

crpQV2 

° 

-1 

0 

0 

0 

0 

0 

Cl 

-Pc«l  - 

-Pe02 

-PeU3 

C2B2 

C2  B3 

Me 

C362 

0 

0 

0 

-c3 

0 

0 

'M  = 

C3#3 

0 

0 

0 

0 

“Cs 

0 

0 

0 

-*r 

0 

0 

0 

0 

0 

0 

0 

-*r 

0 

0 

0 

V  c4 

C5 

C6 

C7 

cs 

eg 

PeUi 

(B22  +  -Bi)(% 

-J). 

C2  = 

±(i- 

Me  ), 

c3  = 

K 

\  * 

v) 

Mo 

/ 

Pqv' 

6? 

\B2u2  +  B3u3),  C5 

=  P- 

~Peu\, 

c6  = 

~PeU\U2  -  - 

b\b2 


B{B: 

Mo 


c8  =  c2u\B2  - 


■#1*02 

Mo^ 


Mo 


c9  =  c2u\Bz  - 


B{uz 

(10V 


We  consider  w  around  a  constant  state 
(9.17) 


w*  =  (u*,  0,0,0  ,v*B*lV*B^E 


*\t 


where  E*  =  e*  +  (B2*2  +  Bf)v*/(2^0),  v*  >  0,  e*  >  0,  B2*  and  B*3  are  constants. 

In  Assumption  2.1  again  we  take  t/(w)  =  -s  and  F(w)  =  0  as  the  entropy  pair. 


Therefore, 


(9-18a)  VI/  =  -s„£l  -  seVe, 

(9.186)  A0  e=  V2t/  =  -svvV\Vl  -  see(Ve)4Ve  -  seV2e  -  s„e((Ve)‘r/i  +  v\' Ve), 
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where 


m  =  o), 


Ve 


=  (— 
V  2^o 


(9.19) 


V2e  = 


■{bI  +  bD 

,-u  i 

,-^2 

J. 

Mo 

B]+Bl 

0 

0 

0 

Bz_ 

Mov 

0 

-1 

0 

0 

Mov 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

-1 

0 

B  2 

o 

0 

0 

_ 1_ 

Mov 

b3 

0 

0 

0 

Mov 

0 

Mov 

0 

0 

0 

0 

0 

MoV 

0 

0 

0 

0 

_ 1_ 

MOV 

0 


o\ 

0 

0 

0 

0 

0 

0  ) 


Clearly  the  sum  of  the  first  three  terms  on  the  right-hand  side  of  (9.18b)  is  positive  definite 
by  (9.10),  (9.4)  and  (9.3).  Notice  the  special  form  of  the  last  term.  To  show  (9.18b)  is 
positive  definite,  we  only  need  to  show  det  A0  >  0.  By  eliminating  the  first  to  sixth  rows 
of  the  second  term  on  the  right,  we  have 

„5 


det  Aq  —  ^  2  ( svvsee  Sve)  '> 

M 0V 


cl.  (9.10).  Hence  U  is  strictly  convex.  From  (9.16),  (9.18a),  (9.19)  and  (9.4)  we  can  vanfy 
{VU)f  =  0.  From  (9.19)  and  (9.15)  we  have  by  direct  computation  that 

(Ve)B  =  6Ve  -1 — 0vt]i. 


Thus  (9.18b),  (9.11)  and  (9.15b)  imply  that 

A0B  =  [-s«(Ve)‘  -  sver)\](Ve)B  -  se{V'2e)B 


(9.20) 


=  -1  S'MVe  +  Si)f  (Ve  +  %i)  -  se(V2e)B. 

V  Ve  ve 


Clearly  the  first  term  on  the  right  is  symmetric  semi-positive  definite.  By  direct  compu- 
tation,  the  second  term  is 


- s , 


VMoV3 

0 

0 

0 


Vr  U 

0 


0 

0 

0 

V 

0 

0 

0 


b2 


0 

0 

0 

0 

0 


S3 

o-MoV3 

0 

0 

0 

0 


0\ 

0 

0 

0 

0 

0 

0J 
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which  is  also  symmetric  semi-positive  definite.  Hence  AqB  is  symmetric  semi-positive 
definite. 

For  Assumption  2.2  we  let 


(9.21) 
Then 

(9.22) 


W  =  w  =  fo(w). 


/o(«>)  = 


1 

0 

0 

0 

B2 

Bs 


(9.23)  B(w)  =  B(w)f^w) 


/0 

0 

0 

0 

0 

0 

0 


i  So  -l-Bi  .  ~ 

V  2  HQ  +  ev  ul  u2 


0 

0 

0 

1 

0 

0 

U3 


0 

0 

0 

0 

V 

0 

B2V 

Mo 


0 

0 

0 

0 

0 

V 

B3v 

Mo 


0 

0 

a 

v 

0 

0 

0 

7W2 


0 

0 

0 

i± 

V 

0 

0 

>3 


0 

0 

0 

0 

1 

<7  fJ>oV 

0 

<7Vov 


0\ 

0 

0 

0 

0 

0 

eej 

0 

0 

0 

0 

0 

1 

afj,Qv 

B, 


0\ 

0 

0 

0 

0 

0 

V 


Clearly  the  null  space  M  of  B(w)  is  independent  of  w. 

Prom  (9.22),  (9.18b), 

(9.4),  (9.5) 

and  (9.11)  we  have 

/ 

e 

0 

0 

0 

0 
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Notice  that  we  have  assumed  k  >  0.  Hence  Af1  is  invariant  under  /o(iu)Mo(w).  It  is  also 
clear  that  the  range  of  B(w)  is  in  A fL. 

Assumption  2.3  has  been  varified  in  Section  6.2  for  the  following  cases: 


(i)  v  >  0,  fi  >  0,  and  1/a  >  0  (finitely  conducting); 

(ii)  v  =  n  =  0,  and  1/a  >  0,  while  ^  0; 

(iii)  v  >  0,  /i  >  0,  and  1/cr  =  0  (perfectly  conducting),  while  B{  ^  0. 
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For  these  cases  our  main  result  applys.  Again  we  notice  that  the  left  eigenspace  of  B(w) 
associated  with  the  zero  eigenvalue  is  independent  of  w.  Thus  if  w(:r,0)  —  w*  G  H4(R) 
and  satisfies  the  smallness  assumptions  in  Theorem  2.3,  the  initial  value  problem  of  the 
equations  of  magnetohydrodynamics  has  an  asymptotic  solution  which  is  a  superposition 
of  diffusion  waves.  From  Section  6.2  A(w*)  has  eigenvalues  in  nondecreasing  order  as 

(9.24)  Cf ,  ca,  cs ,  0,  c$ ,  ca ,  cj , 

where  c/,  ca,  cs  are  the  fast,  the  Alfven,  and  the  slow  wave  speeds,  respectively,  given 
by  (6.63).  If  B{  =  0,  cs  =  ca  =  0.  In  this  case  we  have  a  multiple-mode  diffusion  wave 
along  the  particle  path  dx/dt  =  0,  and  two  single- mode  diffusion  waves  along  the  fast  wave 
directions.  If  B{  ^  0  and  BJ  =  £3  =  0,  either  c/  =  ca  or  cs  =  ca.  In  this  case  we  have  a 
single-mode  diffusion  wave  along  the  particle  path,  and  two  multiple-mode  diffusion  waves 
along  the  Alfven  wave  directions.  In  addition  we  have  two  single-mode  diffusion  waves 
along  the  fast  (slow)  wave  directions  if  —pv  >  <?a  (~pv  <  c^),  where  pv  takes  value  at  the 
constant  state.  If  B{  7^  0  and  B%2  +  B 3 2  ^  0,  all  the  eigenvalues  in  (9.24)  are  distinct, 
and  we  have  single-mode  diffusion  waves  in  all  the  characteristic  directions.  The  difference 
between  the  solution  and  the  asymptotic  solution  is  then  given  by  (2.25). 
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Large  Time  Behavior  of  Solutions  for  General  Quasilinear 
Hyperbolic-Parabolic  Systems  of  Conservation  Laws 
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Physical  models  of  conservation  form— such  as  the  compressible  Navier-Stokes 
equations,  and  MHD  and  full  nonlinear  elasticity  equations — are  not  uniformly 
parabolic,  but  rather  hyperbolic-parabolic.  This  memoir  gives  a  self-contained  anal¬ 
ysis  of  nonlinear  interactions  of  dissipation  waves  as  well  as  the  hyperbolic  aspects 
of  general  systems.  It  introduces  the  new  pointwise  estimates  of  Green  functions 
and  coupling  of  nonlinear  waves. 
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